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ABSTRACT 

The  optical  rotatory  dispersion  of  lactoperoxidase 
and  its  cyanide  and  fluoride  complexes  was  studied  over 
the  spectral  region  210-500  nm,  and  that  of  the  azide 
complex  from  350-450  nm.  Results  of  the  measurements  of 
the  reduced  mean  residue  rotation  at  233  nm  lead  to  the 
conclusion  that  there  are  no  significant  changes  in  protein 
conformation  upon  the  binding  of  these  ligands  to  lacto¬ 
peroxidase.  Lactoperoxidase  was  estimated  to  have  17% 
a-helical  character  at  pH  7.0.  Results  of  studies  in  the 
Soret  region  indicate  that  lactoperoxidase,  unlike  horse¬ 
radish  peroxidase,  hemoglobin,  and  myoglobin,  exhibits  a 
negative  Cotton  effect  as  do  its  cyanide,  azide,  and  fluoride 
complexes.  The  binding  of  these  ligands  apparently  causes 
an  alteration  of  the  geometry  of  the  heme  group  with  respect 
to  the  protein  moiety. 

The  kinetics  of  the  formation  of  the  primary  lactoper- 
oxidase-hydrogen  peroxide  compound  (compound  I)  at  25 
have  been  studied  by  steady  state  methods.  The  second  order 
rate  constant  ^  is  pH-independent  over  the  pH  region 
investigated,  having  a  value  of  (9.2  ±  0.9)  x  10  M  sec 
An  anomalous  effect  of  formate  buffer  on  the  kinetics  of 
the  formation  of  compound  I  is  demonstrated. 

The  kinetics  of  the  oxidation  of  iodide  ion  by  lacto¬ 
peroxidase  compound  II  have  been  studied  as  a  function  of  pH 
at  25°  and  an  ionic  strength  of  0.05.  Both  a  first  order 
and  a  second  order  dependence  on  the  concentration  of  iodide 
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ion  were  detected.  The  second  order  rate  constant  for  the 
reaction  of  lactoperoxidase  compound  II  with  iodide  ion 
decreases  from  4.2  x  10  to  2.6  x  10  M  sec  and  the  third 
order  rate  constant  decreases  from  2.7  x  lO'*'0  to  1.4  x  10^ 

M  sec  with  increasing  pH  over  the  pH  range  2.9  -  10.1. 

The  pH  dependence  of  the  reaction  is  explained  in  terms  of 
an  acid  dissociation  outside  the  pH  range  of  this  study. 

The  kinetics  of  the  oxidation  of  p-cresol  by  compound 
II  of  lactoperoxidase  have  been  studied  over  the  pH  range 
2.1  -  11.2  by  the  stopped-flow  technique.  The  reaction  is 
kinetically  first  order  in  p-cresol  over  the  entire  pH 
range.  Use  is  made  of  the  diffusion  controlled  limit  to 
show  that  p-cresol  reacts  in  the  unionized  form  over  the 
pH  region  of  the  study.  The  complexity  of  the  pH-rate 
profile  is  discussed  in  terms  of  acid  dissociation  constants 
of  groups  in  the  enzyme,  and  the  ionization  of  the  substrate 

The  rate  of  exchange  of  water  with  the  iron  atom  of  the 

heme  group  of  horseradish  peroxidase  was  studied  at  several 
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pH  values  using  a  line  broadening  technique  with  O 
nuclear  magnetic  resonance.  The  results,  while  complicated 
by  experimental  error,  indicate  that  the  true  values  for  the 
rate  of  exchange  are  not  obtainable  in  the  experimental 
temperature  region  5°  —  50°  by  the  nuclear  magnetic  resonance 
line  broadening  technique. 
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Chapter  1 
Introduction 

The  hemoproteins  are  a  distinct  class  of  enzymes  and 
proteins  which  include  the  hemoglobins  and  myoglobins 
concerned  in  oxygen  transport  and  storage,  the  cytochromes, 
which  are  electron  transfer  agents,  and  several  enzymes, 
the  peroxidases  and  catalase.  While  the  physiological  role 
of  the  oxygen  transport,  storage,  and  utilization  pigments 
are  fairly  well  defined,  the  importance  of  peroxidases  and 
catalase  is  not  clear.  The  peroxidases  can  be  classified 
as  enzymes  that  catalyze  the  oxidation  of  certain  substrates 
by  such  peroxides  as  hydrogen  peroxide,  methyl  hydrogen 
peroxide,  and  ethyl  hydrogen  peroxide.  Peroxidases  are 
widespread  in  plant  materials  and  have  been  implicated  in 
the  control  of  plant  growth  (Galston,  1956) ;  they  are  also 
found  in  certain  animal  tissues.  Catalase  is  generally 
present  in  aerobic  cells,  frequently  in  sizeable  amounts. 

The  interest  in  these  enzymes  is  largely  centered  about 
their  mechanisms  of  action,  but  many  reasons  for  the  presence 
of  these  vigorous  catalysts  remain  unsolved. 

Catalase  is  an  enzyme,  the  primary  function  of  which 
appears  to  be  the  destruction  of  hydrogen  peroxide  according 
to  the  equation 

2H202  - >  2H20  +  02 

Catalases  from  different  sources  are  alike,  while  peroxidases 
are  generally  different  from  each  other  and  from  catalase. 
Catalase  from  all  sources  has  a  molecular  weight  of 
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250,000;  the  molecule  is  composed  of  four  subunits,  each 
of  which  has  as  prosthetic  group  the  heme  group  known  as 
f erriprotoporphyrin  IX,  shown  in  Fig.  1-la.  Under  certain 
conditions  catalase  can  act  as  a  peroxidase  by  using 
hydrogen  peroxide  to  oxidize  certain  compounds;  however, 
the  more  generally  observed  reaction  is  the  destruction  of 
hydrogen  peroxide  in  the  so-called  "catalatic"  reaction. 

Peroxidases  possess  some  catalatic  ability  in  the  absence 

-4 

of  reducing  agents,  but  they  are  only  10  as  effective  as 
catalase . 

Peroxidases  are  widely  distributed  among  plants;  the 
richest  sources  are  the  sap  of  the  fig  tree  and  the  root 
of  the  horseradish.  Saunders  et_  al_.  (1964)  give  an 
extensive  list  of  sources  of  peroxidase  activity  in  both 
plants  and  animals.  Some  of  the  peroxidases  which  have  had 
a  considerable  amount  of  work  performed  on  them  in  recent 
years  are  horseradish  peroxidase,  lactoperoxidase ,  thyroid 
peroxidase,  cytochrome  c  peroxidase,  and  chloroperoxidase . 

In  addition,  there  is  another  peroxidase  which  occurs  in 
leucocytes  called  myeloperoxidase  (verdoperoxidase) ;  how¬ 
ever,  relatively  little  work  has  been  done  on  this  enzyme. 

Horseradish  peroxidase  (EC  1.11.1.7;  donor  H2O2  oxido- 
reductase)  is  the  peroxidase  which  has  been  most  extensively 
investigated.  This  enzyme,  which  is  brown,  has  a  molecular 
weight  of  about  40,000  (Keilin  and  Hartree,  1951)  and 
contains  about  18%  carbohydrate  material  by  weight  (Shannon 
et  al.,  1966).  It  contains  as  prosthetic  group  ferriproto- 


3 


porphyrin  IX  (Fig.  1-la) .  Horseradish  peroxidase  was  first 
observed  to  form  spectroscopically  distinct  compounds  with 
hydrogen  peroxide  over  30  years  ago  (Keilin  and  Mann,  1937) . 
Much  work  has  been  done  concerning  the  nature  of  these 
horseradish  peroxidase-hydrogen  peroxide  compounds  by  study¬ 
ing  their  fast  reaction  kinetics,  electron  spin  resonance 
spectroscopy,  and  magnetic  susceptibilities,  and  the 
conclusions  will  be  discussed  in  some  detail  later  in 
connection  with  the  enzyme  lactoperoxidase .  Saunders  et  al . 
(1964) ,  and,  more  recently,  Hasinoff  (1970)  provide  good 
reviews  of  horseradish  peroxidase. 

Cytochrome  c  peroxidase  (EC  1.11.1.5;  cytochrome  c  : 

^2^2  oxa<^orec^uc^::ase )  '  which  catalyzes  the  oxidation  of 
ferrocytochrome  c  to  f erricytochrome  c  in  the  presence  of 
hydrogen  peroxide , 

H2C>2  +  2  ferrocytochrome  c  - >  2  f erricytochrome  c  +  2  OH 

was  discovered  in  baker's  yeast  by  Altschul  et  al .  (1940)  . 

The  enzyme  is  found  exclusively  in  aerobically  grown  yeasts. 
It  has  a  number  of  convenient  properties  that  are  difficult 
to  find  in  other  heme-containing  enzymes.  Large  quantities 
of  cytochrome  c  peroxidase  can  be  prepared  from  commercially 
available  yeasts  by  simple  chromatographic  techniques 
(Yonetani  and  Ray,  1965;  Ellfolk,  1967;  Yonetani ,  1968).  It 
is  readily  crystallized  from  dilute  salt  solutions 
(Yonetani  et  al.,  1966;  Ellfolk,  1967;  Yonetani ,  1968).  The 
crystalline  preparation  of  this  enzyme  is  free  of  isozymes 
and  has  a  consistent  purity.  The  molecular  weight  was 
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determined  to  be  34,100  (Ellfolk,  1967a)  and  the  prosthetic 
group  to  be  f err iprotoporphyrin  IX,  shown  in  Fig.  1-la 
(Altschul  et_  al_.  ,  1940  ;  Abrams  et_  al^.  ,  1942  ;  Yonetani  and 
Ray,  1965).  Unlike  other  peroxidases,  it  has  been  found  to 
contain  no  carbohydrate  material  (Ellfolk,  1967).  As 
reported  originally  by  Altschul  et  al_.  (1940)  and  Abrams 
et  al.  (1942),  the  brown  coloured  cytochrome  c  peroxidase 
reacts  with  a  stoichiometric  amount  of  hydroperoxide  to 
form  a  "red  peroxide  compound",  which  has  two  more  oxidizing 
equivalents  than  the  native  enzyme.  This  red  peroxide 
compound  is  highly  stable  in  the  absence  of  reducing  agents, 
in  contrast  to  the  peroxide  compounds  of  other  peroxidases. 
The  half-life  of  its  decomposition  is  of  the  order  of 
several  hours  at  room  temperature  (Yonetani  et  al . ,  1966). 
This  stability  has  lent  itself  to  electron  spin  resonance 
studies  of  the  peroxide  compound  and  comparisons  with  the 
native  enzyme  (Yonetani  et_  al .  ,  1966a)  .  In  addition, 
Yonetani  (1970)  reports  that  X-ray  studies  of  cytochrome  c 
peroxidase  are  planned;  this  should  give  a  clear  picture  of 
the  heme  environment  of  this  enzyme,  and  will  be  indispen¬ 
sable  for  the  understanding  of  its  mode  of  action. 

In  the  thyroid  gland  iodide  ion  is  oxidized  and 
iodinates  tyrosine  residues  of  the  protein  thyroglobulin . 

Two  iodinated  tyrosine  residues  presumably  couple  to  form 
the  hormone  thyroxine  which  is  secreted  into  the  circulatory 
system.  It  is  then  carried  to  other  tissues  upon  which 
the  hormone  exerts  its  influence.  It  has  long  been  thought 
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that  a  peroxidase  in  the  thyroid  gland  was  involved  in  the 
oxidation  of  iodide  ion.  Klebanoff  et_  ad.  (1962)  showed 
that  beef  thyroid  preparations  were  capable  of  the  iodina- 
tion  of  tyrosine.  Yip  (1966)  studied  the  properties  of  a 
peroxidase  purified  from  beef  thyroid  tissues.  The  puri¬ 
fication  was  largely  effected  by  the  use  of  ion-exchange 
chromatogrphy .  The  molecular  weight  was  determined  to  be 
50,000.  The  prosthetic  group  was  identified  as  ferriproto- 
porphyrin  IX.  Hosoya  and  Morrison  (1967a)  have  described 
the  isolation  and  purification  of  hog  thyroid  peroxidase. 

Gel  filtration  experiments  have  shown  it  to  have  a  molecular 
weight  of  104 , 000 .  The  difference  between  this  figure  and 
that  reported  for  beef  thyroid  peroxidase  may  be  due  to 
species  difference  or  aggregation.  Hosoya  and  Morrison 
(1967a)  also  found  that  hog  thyroid  peroxidase  reacts  with 

hydrogen  peroxide  about  as  quickly  as  horseradish  peroxidase 

7-1-1 

and  lactoperoxidase  (k  =  1.1  x 10  M  sec  at  pH  7.4)  and 
that  this  enzyme  can  iodinate  tyrosine  and  thyroglobulin. 
Taurog  et  al.  (1970)  have  described  an  improved  procedure 
for  the  solubilization  and  purification  of  hog  thyroid 
peroxidase;  they  obtained  29  milligrams  of  purified  enzyme 
from  23  kg  of  frozen  hog  thyroid  glands.  Their  spectral 
data  suggest  that  the  heme  in  thyroid  peroxidase  may  not 
be  f err iprotoporphyrin  IX.  Taurog  (1970)  has  shown  that 
the  thyroid  peroxidase-catalyzed  iodination  of  thyroglobulin 
is  inhibited  by  excess  iodide.  These  findings  offer  a 
possible  explanation  for  inhibitory  effects  of  large  doses 
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of  iodide  ion  on  thyroid  hormone  formation  in  humans 
(iodide-induced  myxedema) .  Most  of  the  problems  involved 
in  work  with  thyroid  peroxidase  may  be  the  result  of  poor 
yields  of  impure  preparations  ,  and  much  work  remains  to  be 
done  on  this  important  enzyme. 

Morris  and  Hager  (1966)  have  described  the  isolation 
and  properties  of  chloroperoxidase .  This  enzyme  is 
present  in  the  mold  Caldar iomyces  f umago ,  which  synthesizes 
the  chlorine-containing  fungal  metabolite  called  caldario- 
mycin.  They  determined  that  the  prosthetic  group  of 
chloroperoxidase  is  f erriprotoporphyrin  IX,  and  that  the 
molecular  weight  is  42,000.  They  also  determined  the  amino 
acid  composition  of  the  enzyme  and  noted  that  25-30%  of 
the  molecule  is  composed  of  carbohydrate  material.  Chloro¬ 
peroxidase,  unlike  horseradish  peroxidase,  can  oxidize 
chloride  and  bromide  ions;  Hager  el:  al .  (1966)  have  fol¬ 

lowed  the  chloroperoxidase-catalyzed  halogenation  of  tyrosine 
and  have  shown  that  the  relative  activity  among  the  halogen 
anions  for  the  halogenation  of  tyrosine  is  approximately 
5.1  :  4.8  :  1  for  iodide,  bromide,  and  chloride  ions,  respect¬ 
ively.  Thomas  and  Hager  (1968)  presented  evidence  that 
chloroperoxidase  can  catalyze  the  oxidation  of  iodide  ion 
beyond  the  molecular  halogen  level  to  iodate ,  i.e.,  to  an 
oxidation  state  of  +5.  More  recently,  Thomas  et  al.  (1970) 
have  reported  on  the  formation  of  peroxide  and  halide 
complexes  of  chloroperoxidase,  and  their  relation  to  the 
mechanism  of  the  halogenation.  In  addition,  Thomas  e t  al « 
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(1970a)  have  shown  that  chloroperoxidase  can  release 
oxygen  from  peroxy-acids ,  a  finding  which  has  some  bearing 
on  the  nature  of  peroxidase-peroxide  compounds.  This 
point  will  be  discussed  later  in  connection  with  lacto- 
peroxidase . 

The  subject  of  most  of  this  thesis  is  lactoperoxidase , 
isolated  from  bovine  milk,  and  which  has  also  been  detected 
in  bovine  salivary  glands  (Morrison  et  al_.  ,  1965)  and 
bovine  harderian  and  lacrimal  glands  (Morrison  and  Allen, 
1966) . 

Historical  Introduction  of  Lactoperoxidase 

While  searching  for  a  method  by  which  fresh  milk  could 
be  distinguished  from  cooked  milk,  Arnold  (1881)  discovered 
that  fresh  bovine  milk  contained  a  substance  which  rapidly 
oxidized  guaiacol  when  hydrogen  peroxide  was  added. 

Previous  studies  had  shown  that  oxidation  of  dyes  could  be 
carried  out  by  various  plant  materials  if  hydrogen  peroxide 
were  added  to  the  reaction  mixture  (Schonbein,  1863) . 
Linnosier  (1881)  demonstrated  the  presence  of  a  substance 
with  similar  activity  in  leucocytes  and  gave  the  name 
"peroxidase"  to  all  of  those  materials  which  catalyzed  the 
oxidation  of  dyes  in  the  presence  of  hydrogen  peroxide. 

For  some  time  the  peroxidase  of  milk  was  assumed  to 
be  identical  to  the  leucocyte  peroxidase;  its  occurence 
in  milk  was  believed  to  result  from  the  disintegration  of 
leucocytes.  However,  Theorell  and  Akeson  (1943)  demon¬ 
strated  that  the  two  peroxidases  had  different  properties. 


8 


These  authors  suggested  that  the  two  green  coloured 
animal  peroxidases  be  classified  as  "verdoperoxidase" 

(or  "myeloperoxidase")  for  the  leucocyte  peroxidase,  and 
that  the  milk  enzyme  be  named  "lactoperoxidase " . 

Isolation  and  Purification 

The  first  procedure  for  the  isolation  of  lactoper¬ 
oxidase  was  developed  by  Thurlow  (1925).  The  peroxidase 
activity  was  obtained  in  a  crude  fraction  by  the  use  of 
ammonium  sulfate  fractionation.  Elliot  (1932)  elaborated 
on  this  procedure  and  obtained  the  activity  in  a  brown 
coloured  fraction.  Yakushiji  (1939)  further  purified  the 
peroxidase  by  subjection  of  the  fraction  which  had  been 
obtained  by  ammonium  sulfate  precipitation  to  precipitation 
by  acetone . 

Lactoperoxidase  was  finally  obtained  in  pure  form  and 
crystallized  by  Theorell  and  co-workers  (1943,  1944). 

Starting  with  skim  milk  these  workers  obtained  a  crude 
fraction  of  lactoperoxidase  by  the  ammonium  sulfate  fraction¬ 
ation  procedure  of  Elliot  (1932)  .  Inert  proteins  were 
removed  from  this  crude  fraction  first  by  heat  denaturation 
at  70°,  followed  by  precipitation  with  basic  lead  acetate, 
and  then  by  acetone  precipitation.  Repeated  electrophoresis 
at  pH  5.9  was  required  to  remove  a  contaminating  red 
protein.  The  lactoperoxidase,  thus  purified,  was  crystal¬ 
lized  at  its  isoelectric  point  from  ammonium  sulfate  solution. 
The  crystals  were  in  the  form  of  thin  leaves. 
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The  isolation  and  purification  procedures  of  early 
workers  were  never  fully  detailed,  however ,  and  the 
reported  yields  we re  low,  although  they  did  obtain  crystals 
of  the  enzyme.  Polis  and  Shmukler  (1953)  developed  a 
more  reproducible  procedure  in  which  the  final  purification 
was  achieved  by  chromatography  on  calcium  phosphate  columns. 
Morrison  et  al .  (1957)  developed  a  method  employing  ion 

exchange  resins  for  the  isolation  of  the  enzyme.  Morrison 
and  Hultqui st  (1963)  later  simplified  and  improved  their 
previous  procedure,  so  that  yields  of  a  more  highly  purified 
preparation  were  obtained.  This  improved  procedure  is 
described  in  Appendix  1,  along  with  criteria  used  for  the 
purity  and  activity  of  samples  of  lactoperoxidase .  Although 
several  criteria  of  purity  are  used,  the  ratio  of  the 
absorbance  at  412  nm  to  that  at  280  nm  (called  P.N. ,  the 
"purity  number")  has  been  the  most  common  criterion  for 
comparing  the  purity  of  various  preparations  of  lacto¬ 
peroxidase.  A  P.N.  value  of  0.9  denotes  a  quite  pure  prep¬ 
aration  of  lactoperoxidase. 

Rombauts  e^  al.  (1967)  have  determined  the  composition 
of  lactoperoxidase.  The  molecular  weight  is  about  77,500 
and  there  is  one  iron  atom  per  enzyme  molecule.  A  quantita¬ 
tive  determination  showed  that  all  the  common  amino  acids 
were  present.  Carbohydrate  material  constituted  8%  by 
weight  of  the  enzyme  molecule.  Table  1-1  shows  the  amino 
acid  and  carbohydrate  composition  of  lactoperoxidase. 


Amino  Acid  Composition  of  Lactoperoxidase 
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Carlstrom  (1965,1966,1969)  has  claimed  that  lactoperoxidase 
exists  as  a  number  of  closely  related  isozymes;  however, 
Rombauts  et  al.  (1967)  have  argued  that  Carlstrom' s 
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findings  are  the  result  of  enzyme  degredation  during  his 
extensive  purification  and  separation  procedures. 

The  nature  of  the  prosthetic  group  (heme  group)  of 
lactoperoxidase  has  not  been  established.  Hultquist  and 
Morrison  (1963)  have  indicated  that  the  heme  prosthetic 
group  is  probably  a  derivative  of  mesoheme  9,  in  which 
two  double  bonds  are  in  conjunction  with  the  tetrapyrrol 
nucleus  and  one  or  more  hydroxyl  groups  are  attached  to 
the  side  chains.  The  linkage  between  hemin  and  protein  was 
identified  as  an  ester  bond.  Figure  1-la  shows  the  heme 
prosthetic  group  of  horseradish  peroxidase,  hemoglobin, 
and  myoglobin;  f erriprotoporphyrin  IX,  which  is  closely 
related  to  mesoheme  9  (Falk,  1964) .  Figure  1-lb  shows  the 
structure  of  the  heme  prosthetic  group  of  lactoperoxidase 
suggested  by  Hultquist  (1962) .  The  difference  between  the 
two  heme  groups  is  manifested  by  different  absorbance 
characteristics  in  the  native  enzymes  and  the  complexes 
they  form. 

Reaction  of  Lactoperoxidase  with  Oxidizing  Agents 

Spec trophotome trie  studies  have  demonstrated  the 
existence  of  enzyme-substrate  compounds  between  lactoper¬ 
oxidase  and  peroxides,  similar  to  the  complexes  which  had  been 
previously  observed  for  horseradish  peroxidase  (Keilin  and 
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Fig.  1-la.  Ferriprotoporphyrin  IX. 


Fig.  1-lb. 


A  proposed  structure  for  the  prosthetic 
of  lactoperoxidase .  (Hultquist,  1962). 


group 
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Mann,  1937  ;  Theorell ,  1942).  The  designations  compound  I 
and  compound  II  refer  to  those  spectroscopically  distinct 
species  that  form, in  order  of  appearance,  after  the 
addition  of  hydrogen  peroxide  to  peroxidases.  Using  a 
rapid  flow  spectrophotometer,  Chance  (1949,  1950)  was  able 
to  demonstrate  the  presence  of  the  labile  compound  I  of 
lactoperoxidase .  He  also  demonstrated  that  compounds  I  and 
II  of  lactoperoxidase  were  formed  with  methyl  hydrogen 
peroxide  and  ethyl  hydrogen  peroxide.  The  spectra  of  these 
compounds  could  not  be  distinguished  from  those  of  the 
hydrogen  peroxide  compounds.  Figure  1-2  shows  the  spectra 
in  the  400  nm  region  (the  Soret  region)  of  the  primary  (I) 
and  the  secondary  (II)  lactoperoxidase-hydrogen  peroxide 
compounds,  and  the  Soret  band  of  lactoperoxidase  (A). 

Chance  (1950)  also  titrated  lactoperoxidase  with  hydrogen 
peroxide,  and  this  experiment  revealed  that  one  mole  of 
peroxide  had  combined  with  one  mole  of  enzyme. 

Kinetic  studies  with  the  rapid  flow  spectrophotometer 
demonstrated  that  compounds  I  and  II  were  involved  in 
enzyme  action  (Chance,  1949,1950) .  From  the  results  of 
these  studies  Chance  concluded  that  the  reaction  mechanism 
for  lactoperoxidase  was  identical  to  that  which  had  been 
postulated  for  horseradish  peroxidase  (Chance,  1949a),  and 
which  may  be  represented  by  the  following  scheme : 


Fig.  1-2. 


The  Soret  bands  of  Lactoperoxidase  (A)  , 
compound  I,  and  compound  II. 


afrom  Chance  (1950) . 
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LP 

+  H2°2 

LP-I 

+  AH  £ 

LP-II 

+  ah2 

->  LP-I 


LP-II  +  AH 


->  LP  +  AH  * 


2  AH 


->  A  +  AH  (or  AH  - 


(1) 

(2) 

(3) 

AH)  (4) 


where  LP  represents  the  native  enzyme,  LP-I  and  LP-II 
are  the  oxidized  forms  of  the  enzyme  referred  to  as 
compounds  I  and  II,  respectively,  and  AH2  is  an  oxidizable 
substrate,  with  A  the  product  of  oxidation.  Various  kinds 
of  substrates  have  been  used  in  lactoperoxidase-catalyzed 
oxidations,  among  them,  phenols,  amines,  and  inorganic 
ions  such  as  iodide  ion  and  ferrocyanide  ion. 

The  reaction  of  lactoperoxidase  with  hydrogen  peroxide 
to  give  compound  I  was  shown  to  be  very  rapid.  The  somewhat 
slower  rate  of  reaction  with  methyl  hydrogen  peroxide  and 
the  still  slower  rate  of  reaction  with  ethyl  hydrogen  peroxide 
suggested  that  there  was  a  correlation  between  the  size  of 
the  substrate  molecule  and  the  rate  of  reaction.  Thus  it 
appeared  that  the  iron  of  lactoperoxidase  was  not  completely 
accessible  to  substrate. 

The  nature  of  the  intermediates  involved  in  oxidations 
by  horseradish  peroxidase  and  lactoperoxidase  is  still  a 
subject  of  some  uncertainty.  Virtually  no  work  has  been 
performed  with  lactoperoxidase  in  this  regard,  but  since 
there  is  believed  to  be  a  general  mechanism  for  the  action 
of  peroxidases,  it  may  be  considered  valid  to  extend  the 
conclusions  obtained  from  horseradish  peroxidase  to  lacto- 
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peroxidase . 

Compound  I  of  horseradish  peroxidase  (HRP-I)  retains 
the  two  oxidizing  equivalents  of  the  hydroperoxide  used 
in  its  formation,  and  is  in  an  oxidation  state  of  V 
(Chance,  1952).  Compound  II  of  horseradish  peroxidase 
(HRP-II )  can  be  formed  from  compound  I  and  is  in  an 
oxidation  state  of  IV  (George,  1953) .  A  number  of  struc¬ 
tures  have  been  proposed  for  HRP-I  and  HRP-II.  For  HRP-II 
these  include  the  Fe(IV)  and  the  Fe(IV)0  (the  ferryl  ion) 
structures  proposed  by  George  (1953a)  .  Brill  and  Williams 
(1961)  agreed  that  the  ferryl  structure  describes  HRP-II, 
but  proposed  for  HRP-I  the  existence  of  two  components  (1) 
a  simple  complex  with  hydroperoxide  in  the  sixth  coordina¬ 
tion  position  and  (2)  an  oxidized  porphyrin  ring  compound. 
Brill  and  Sandberg  (1968)  later  put  forward  proposals  for 
the  structure  of  HRP-I  including  (1)  iron  in  a  quadrivalent 
state  plus  a  free  radical  or  (2)  iron  in  the  ferric  state 
plus  a  biradical  on  the  porphyrin  or  protein.  Moss  et  al . 
(1969)  stated  that  the  results  of  their  Mossbauer  study 
on  HRP-I  and  HRP-II  were  compatible  with  a  Fe(IV)  type 
structure  for  HRP-II,  but  they  also  concluded  that  there 
was  no  difference  in  the  configuration  of  the  iron  between 
HRP-I  and  HRP-II,  indicating  that  the  second  equivalent 
on  HRP-I  is  not  stored  on  the  iron.  Peisach  et  a_l.  (1968) 
and  Blumberg  et  al_.  (1968)  proposed  structures  of  HRP-I 
and  HRP-II  consistent  with  low  temperature  optical  spectra, 
electron  paramagnetic  resonance  spectra,  and  previously 
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reported  magnetic  susceptibility  data.  For  HRP-II  they 
postulate  that  two  oxygen  atoms  reside  in  the  sixth 
ligand  position  of  the  iron  atom: 

[HRP  (heme  dQ6  )”00H]IV 

For  the  structure  of  HRP-I  it  was  suggested  that  the  Z_ 
axial  ligands  are  not  present  or  that  they  provide  only 
a  weak  component  to  the  ligand  field  consistent  with: 

[HRP  (heme35/2)  ] 111 ,  [X] 11 

The  two  oxidizing  equivalents  are  stored  not  on  a  ligand 
but  either  on  the  protein  or  the  porphyrin  represented  as 
group  X. 

Dolphin  et  aJL.  (1971)  contend  that  the  optical  spectrum 
of  HRP-I  characterizes  it  as  a  porphyrin  7T-cation  radical 
because  of  its  similarity  to  the  optical  spectrum  of  a 
cobalt  porphyrin  Tr-cation  radical.  Thus  one  of  the  two 
oxidizing  equivalents  of  HRP-I  is  accounted  for  by  the  loss 
of  an  electron  from  the  porphyrin  r-orbitals ,  and  the 
second  oxidation  may  then  occur  from  the  porphyrin  ring,  the 
protein,  or  the  metal.  Thomas  et_  a_l.  (197  0a)  have  shown 
that  a  unique  reaction  catalyzed  by  a  related  enzyme, 
chloroperoxidase ,  is  the  evolution  of  oxygen  from  substituted 
peroxides  such  as  ethyl  hydrogen  peroxide  and  peroxy  acids . 
Although  horseradish  peroxidase  reacts  with  the  latter  two 
compounds  to  form  compound  I ,  it  can  not  complete  a  cycle 
of  oxygen  formation  and  regeneration  of  the  native  enzyme 


19 


(Chance,  1948,1949b).  Provided  that  this  oxygen  evolution 
reaction  proceeds  through  a  compound  I  intermediate, 
def ini te  conclusions  regarding  the  chemical  nature  of 
compound  I  might  be  drawn  once  the  origin  of  the  oxygen 
atoms  are  known.  Along  this  line  of  approach,  Hager  et  al . 
(1971)  have  shown  that  the  oxygen  atoms  in  molecular  oxygen 
evolved  from  peroxy  acids  arise  from  different  substrate 
molecules.  This  leads  to  the  conclusion  that  compound  I 
must  retain  at  least  one  oxygen  atom  after  reaction  with 
the  first  peroxy  acid  molecule.  In  a  subsequent  step, 
compound  I  reacts  with  a  second  peroxy  acid  molecule  to 
form  molecular  oxygen. 

Further  work  is  necessary  to  define  the  chemical 
nature  of  peroxidase-peroxide  compounds;  however,  it  should 
be  emphasized  that  any  conclusions  regarding  the  nature  of 
HRP-I  and  HRP-II  and  compound  I  of  chloroperoxidase  may 
not  necessarily  apply  to  the  structures  of  LP-I  and  LP-II. 

The  characteristic  property  and  function  of  enzymes 
is  the  catalysis  of  chemical  reactions.  Any  fundamental 
study  of  this  catalytic  function  must  be  based  on  quantita¬ 
tive  measurements  of  the  rate  of  the  catalyzed  reaction. 

From  the  effect  of  varying  the  reaction  conditions  on  the 
rate,  inferences  may  be  made  about  the  mechanism  of  enzyme 
action.  Ideally  such  kinetic  studies  should  be  brought 
into  relation  with  chemical  and  structural  studies  on  the 
enzyme  in  order  to  obtain  a  definite  picture  of  the  process, 
but  this  is  only  possible  if  the  enzyme  has  been  obtained 
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in  a  high  degree  of  purity.  Many  enzymes  have  not  been  so 
purified  and  kinetic  studies  are  the  only  approach  possible 
at  present.  It  is  well  known  that  biological  systems, 
particularly  living  cells,  are  more  sensitive  to  changes 
of  temperature  and  pH  than  most  non-biological  chemical 
reactions,  and  this  is  due  largely  to  the  properties  of 
the  enzymes  on  which  these  systems  depend.  A  knowledge 
of  enzyme  kinetics  thus  helps  in  the  understanding  of 
biological  phenomena.  A  particularly  important  tool  in  the 
investigation  of  the  kinetics  of  enzyme  catalyzed  reactions 
is  the  variation  of  the  pH  of  the  reaction  medium,  and  it 
is  this  approach  that  is  used  in  the  kinetic  investigations 
described  in  this  thesis. 
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Chapter  2 

The  Effect  of  Ligand  Binding  on  the  Optical 
Rotatory  Dispersion  of  Lactoperoxidase 

Introduction 

A  beam  of  plane  polarized  light  may  be  considered  to 
be  composed  of  two  circularly  polarized  rays,  one  whose 
vector  rotates  clockwise  and  one  counterclockwise  as  the 
beam  advances.  A  medium  containing  asymmetric  molecules 
transmits  the  two  components  with  unequal  velocity,  that 
is,  it  has  different  refractive  indices  for  the  two 
circularly  polarized  rays.  If  neither  ray  is  absorbed, 
the  two  rays  have  a  resultant  which  is  a  plane-polarized 
ray  whose  vibration  is  in  a  plane  which  has  been  rotated 
with  respect  to  that  of  the  incident  ray.  The  plane  of 
polarization  moves  through  a  steadily  increasing  angle 
as  the  beam  passes  through  the  medium.  The  amount  of 
rotation  is  measured  as  the  angle  between  the  incident 
and  emergent  beams,  and  is  normally  expressed  as  specific 
rotation,  [a]^,  which  is  given  by 


where  X  and  t  specify  the  wavelength  at  which  measurements 
are  made,  and  the  temperature,  respectively,  a  is  the 
observed  rotation  in  degrees ,  1  is  the  path  length  of  the 
sample  in  decimeters,  and  c  is  the  number  of  grams  of 
solute  in  one  ml  of  solution.  Since  the  indices  of 
refraction  are  dependent  on  wavelength,  different  wave- 


lengths  of  light  are  rotated  by  different  amounts, 
giving  rise  to  the  phenomenon  known  as  optical  rotatory 
dispersion . 

If  the  medium  also  absorbs  the  two  circularly 
polarized  components  unequally,  the  emergent  light  beam 
is  elliptically  polarized  and  the  medium  exhibits 
circular  dichroism.  The  combination  of  unequal  absorp¬ 
tion  (circular  dichroism)  and  unequal  velocity  of  trans¬ 
mission  (optical  rotation)  of  left  and  right  circularly 
polarized  light  in  the  region  in  which  optically  active 
absorption  bands  are  observed  is  a  phenomenon  called 
the  Cotton  effect.  The  O.R.D.  curve  of  a  Cotton  effect 
is  a  symmetrical  curve  shaped  like  a  derivative  of  an 
absorption  peak  and  is  centered  at  a  wavelength  corre¬ 
sponding  to  the  maximum  of  the  absorption  peak.  The 
maximum  of  a  Cotton  effect  curve  is  known  as  the  peak, 
and  the  minimum  as  the  trough.  When  the  peak  occurs 
at  a  wavelength  longer  than  that  of  the  trough,  the 
Cotton  effect  is  termed  "positive";  in  the  reverse 


"^Abbreviations  used  in  this  chapter:  O.R.D.  ,  optical 
rotatory  dispersion:  LP ,  LP-CN,  LP-N^ ,  and  LP-F ,  ferric 
lactoperoxidase  and  its  cyanide,  azide,  and  fluoride 
complexes,  respectively;  PGA,  poly-a-L-glutamic  acid; 
P.N.,  (purity  number),  ratio  of  the  absorbance  of  a 
solution  of  LP  at  412  nm  to  the  absorbance  at  280  nm. 
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situation,  the  Cotton  effect  is  "negative". 

Cotton  effects  have  been  classified  into  two  groups. 
Using  the  terminology  suggested  by  Blout  (1964), 

"intrinsic  Cotton  effects"  reflect  the  optical  rotatory 
power  of  the  protein  itself,  and  "extrinsic  Cotton 
effects"  are  generated  by  a  chromophoric  site  on  the 
protein  or  by  a  chromophoric  molecule  interacting  with 
an  asymmetric  site  on  the  molecule. 

Optical  rotation,  generated  by  a  metal-protein 
complex,  when  manifested  as  a  specific  Cotton  effect, 
may  reflect  simultaneously  the  chemical  composition  of 
the  metal-protein  ligand  site,  its  configuration,  and 
also  the  spatial  disposition  of  the  chromophoric  site 
with  respect  to  neighboring  groups.  Hence,  studies  of 
optical  rotatory  dispersion  of  metal-protein  systems 
disclose  features  in  addition  to  those  discernible  from 
absorption  spectra.  The  heme  proteins  serve  to  illus¬ 
trate  this  point. 

The  conformation-dependent  Cotton  effects  have 
provided  a  means  for  characterizing  secondary  structures. 
The  depth  of  the  233  nm  trough  has  been  suggested  as  a 
measure  of  the  helical  content  in  proteins  (Simmons  et.  al . 
1961).  However,  its  quantitative  significance  is  still 
not  certain  because  of  the  lack  of  agreement  in  reference 
scales . 

The  233  nm  trough  method  has  been  used  to  study  the 
effect  of  ligand  binding  on  the  conformation  of  heme 
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proteins.  It  has  been  found  that  while  myoglobin  shows 
no  change  in  helical  content  upon  the  binding  of  ligands 
(Samejima  and  Yang,  1964;  Breslow  et  al . ,  1965),  and 
while  hemoglobin  shows  no  change  upon  the  binding  of 
most  ligands  (Beychok,  1964) ,  hemoglobin  seems  to  undergo 
an  8.5%  decrease  in  helicity  upon  the  binding  of  oxygen 
(Brunori  et  al. ,  1967) .  A  similar  study  on  horseradish 
peroxidase  (Ellis  and  Dunford,  1968)  showed  no  change  in 
helical  content  upon  the  binding  of  the  ligands  fluoride, 
cyanide,  and  hydroxide  ions.  The  present  study  was 
undertaken  with  a  view  to  determining  whether  a  gross 
conformational  change  takes  place  upon  the  binding  of 
ligands  to  lactoperoxidase .  This  chapter  describes  in 
greater  detail  the  work  already  published  (Maguire  and 
Dunford ,  1971) . 

Experimental  methods 

Lactoperoxidase  was  isolated  by  the  method  of 
Morrison  and  Hultquist  (1963) .  Further  purifications 
were  carried  out  by  column  chromatography  with  Sephadex 
G-100  and  G-200  gels.  Activity  tests  were  performed  on 
lactoperoxidase  using  the  guaiacol  method  of  Maehly  and 
Chance  (1954) .  LP  solutions  were  dialyzed  against  the 
appropriate  buffer  and  passed  through  a  Millipore  filter 
before  use.  The  P.N.  of  all  solutions  was  at  least  0.82. 
The  maximum  P.N.  attained  was  0.93.  All  buffers  were  of 
an  ionic  strength  of  0.05.  Buffers  used  in  the  order  of 
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increasing  pH  were:  glycine-HCl,  formate,  acetate, 
phosphate,  tris-HCl ,  glycinate,  sodium  hydroxide-disodium 
hydrogen  phosphate,  and  sodium  hydroxide-potassium  chloride. 
Overlapping  buffer  systems  were  used  to  check  for  possible 
buffer  effects.  Inorganic  chemicals  were  of  reagent 
grade  and  were  used  without  further  purification. 

Lactoperoxidase  concentrations  were  determined  from 
absorbance  measurements  at  412  nm  using  a  Cary  14  spectro¬ 
photometer.  The  molar  absorptivity  used  was  = 

1.14  x  10^  M  ^cm  (Morrison  et  al_.  ,  1957)  .  Typical 
concentrations  of  lactoperoxidase  used  were  7  x  10  6  M 
for  O.R.D.  measurements  in  the  ultraviolet  region  and 
5  x 10  M  for  measurements  in  the  visible  region.  All  pH 
measurements  were  made  with  an  Orion  model  801  digital  pH 
meter  equipped  with  a  Fisher  combination  electrode. 

Optical  rotatory  dispersion  spectra  were  obtained 
using  a  Jasco  model  O.R.D./UV5  recording  spectrometer. 

Both  1  mm  and  5  mm  cells  were  used.  At  least  five  O.R.D. 
determinations  were  carried  out  on  the  LP  samples  at  each 
pH  value,  except  for  pH  7.0  where  14  determinations  were 
made.  Multiple  determinations  were  carried  out  to 
determine  the  reproducibility  of  the  O.R.D.  spectrum  for 
a  given  sample  (maximum  observed  deviation  of  1%)  and  to 
check  the  variation  in  spectrum  with  change  in  concentra¬ 
tion  of  lactoperoxidase. 
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The  LP-F  complex  was  studied  at  pH  4.5  in  a  solution 
which  was  0.025  M  with  respect  to  NaF ,  LP-CN  at  pH  7.0 
with  0.01  M  KCN ,  and  LP-N3  at  pH  3.9  with  0.5  M  NaN3 .  All 
spectra  of  the  complexes  were  taken  six  times  and  were 
reproducible  to  within  2%. 

For  all  solutions  used  in  this  study,  the  maximum 
absorption  was  less  than  two  in  the  wavelength  regions  of 
interest.  This  procedure  was  followed  to  help  ensure  that 
the  O.R.D.  spectra  obtained  were  not  caused  by  rotatory 
artifacts  produced  by  regions  of  high  absorbance  (Urnes 
and  Doty,  1961).  The  values  of  the  specific  rotation  [a] 
were  calculated  using  the  molar  concentrations  determined 
spectrophotometrically  and  a  value  of  78,000  for  the 
molecular  weight  of  lactoperoxidase .  All  measurements 
were  made  at  25°. 

Results 

Figure  2-1  shows  the  O.R.D.  curve  for  lactoperoxidase 
over  the  v/avelength  range  210-350  nm,  plotted  as  corrected 
specific  rotation,  [a'],  vs.  wavelength.  Within  experimen¬ 
tal  error,  the  O.R.D.  curves  of  the  fluoride  and  cyanide 
complexes  of  LP  are  the  same  as  that  of  pure  LP .  The 
O.R.D.  spectrum  of  LP-N3  could  not  be  obtained  in  this 
region  because  of  the  high  absorbance  of  the  azide  ion.  The 
curve  in  Fig.  2-1  shows  a  minimum  at  233  nm  which  is 
characteristic  of  the  conformation-dependent  trough  of  the 
negative  Cotton  effect  centered  at  224  nm,  which  is 
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exhibited  by  poly-a-amino  acids  in  the  a-helical  conforma¬ 
tion.  For  a  quantitative  measure  of  the  amplitude  of 
this  trough,  we  have  used  the  reduced  mean  residue  rotation 
[r']  233  /  as  suggested  by  Simmons  et  cLL.  (1961).  The  values 
of  [R  ] 2 3 3  f°r  LP  and  its  cyanide  and  fluoride  complexes 

were  calculated  from  the  relation 

M 

[ R  ^  233  =  100  2  ~Z  ^  233 

n2  3  3  +  2 

where  the  refractive  index  of  water,  n22?>,  was  taken  as 
1.39,  and  the  mean  residue  molecular  weight,  M  ,  as  115. 
Lactoperoxidase  contains  about  8%  carbohydrate  material 
by  weight  (Rombauts  et  al_.,  1967),  and  thus  a  corrected 
specific  rotation  t 06 '  ]  2  3  3  Was  usec*  basecl  on  an  effective 
protein  molecular  weight  of  72,000.  Table  2-1  contains 
the  values  obtained  for  C R *  1 2 3 3  and  percent  helix  content 
for  pure  LP  over  the  pH  range  3.4-10.3.  Table  2-2  contains 
similar  O.R.D.  results  for  LP-F  and  LP-CN  at  233  nm. 

The  values  of  percent  helix  show  no  trend  with  pH. 

In  general,  an  LP  solution  with  a  higher  P.N.  exhibited 
a  slightly  lower  helix  content  (about  2%  in  going  from  P.N. 
0.83  to  0.88)  .  This  may  be  due  to  degradation  as  a  result 
of  extensive  purification  procedures.  The  use  of  [ R *  1 2 3 3 
to  calculate  the  helical  content  of  proteins  seems  to 
result  in  an  underestimate  in  many  cases  (Urnes  and  Doty, 
1961) ,  which  may  be  due  to  the  fact  that  poly-a-L-glutamic 
acid  has  been  used  as  the  standard  for  100%  helicity. 
Although  there  seems  to  be  agreement  that  [ R '  1333  ecIual 
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Table  2-1.  [R  ]233 

a  function  of  pH. 

and  percent  helix 

for  LP  solutions 

pH 

"~^R  ^233 

%  helix 

3.4 

4480 

19 

4.0 

3490 

11 

4.6 

4190 

17 

o 

• 

LO 

3930 

15 

5.6 

5100 

24 

6.1 

3940 

15 

7.0 

4270 

17 

7.6 

4360 

18 

8.0 

4040 

16 

8 . 2 

4490 

19 

8 . 8 

3710 

13 

9.1 

3360 

10 

9.5 

4650 

20 

9.9 

4030 

16 

10.3 

4510 

19 

30 


Table  2-2.  Reduced  mean  residue  rotation  at  233  nm  and 
helix  content  for  lactoperoxidase  and  its  fluoride  and 
cyanide  complexes. 


Species 

PH 

“ [R  J  233 

%  helix 

LP 

7.0 

4270 

17 

LP-CN 

7.0 

4190 

17 

LP 

4.5 

4080 

16 

LP-F 

4.5 

4220 

17 
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to  -2000°  for  random  PGA,  there  is  a  wide  disparity  in  the 
reported  values  for  completely  helical  PGA  (Simmons  et  al . 
1961;  Jirgensons,  1965  ;  Yang  and  Samejima  ,  1963)  .  After 


Yang  (1967),  we  will  assume  that  [R  ] 2  3  3  (helix)  =  -15,000 
and  [ R *  ]  333  (c°il)  =  -2000°  and  obtain  a  rough  estimate  of 
the  a-helical  content  by  simple  interpolation: 


/ 


O 


fraction  of  a-helix 


( [ R ' ]233  +  2,000) 

13,000 


(3) 


The  randomness  of  the  [ R * 1333  results  would  probably 
have  obscured  any  concentration  or  buffer  effect,  if  these 
were  present.  At  either  extreme  of  pH,  denaturation 
occurred.  O.R.D.  determinations  were  not  carried  out  below 
pH  3.4  or  above  10.4.  Figure  2-2  shows  the  extrinsic 
Cotton  effects  associated  with  the  heme  Soret  bands  of 
lactoperoxidase  and  its  fluoride,  cyanide,  and  azide 
complexes . 

The  activity  of  lactoperoxidase  was  determined  by  the 
oxidation  of  guaiacol  (Maehly  and  Chance,  1954),  which  was 
monitored  spectrophotometrically  at  470  nm.  It  was 
observed  that  the  activity  of  lactoperoxidase  was  independ¬ 
ent  of  the  P.N.  exhibited  (from  P.N.  0.50  to  0.90), 
indicating  that  protein  or  peptide  contaminants  do  not 
interfere  with  the  oxidation  reaction. 


Discussion 

The  values  of  [R  ]  222  f°r  LP-ligand  complexes 

were  about  the  same  as  for  pure  lactoperoxidase,  the 
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Fig.  2-2.  Optical  rotatory  dispersion  in  the  Soret  region  of  lactoperoxidase 
and  its  fluoride  (—  —  —),  cyanide  ( - )  ,  and  azide  ( - )  complexes. 


difference  being  within  experimental  error.  It  is 
concluded  that  the  binding  of  these  ligands  to  lactoperox- 
idase  has  no  significant  effect  on  the  amount  of  a-helical 
conformation  of  the  enzyme.  Lactoperoxidase  is  calculated 
to  contain  17%  a-helix  at  pH  7. 

The  data  in  Fig.  2-1  indicate  the  presence  of  Cotton 
effects  due  to  side-chain  chromophores  in  the  260-290  nm 
region.  They  were  studied  in  some  detail  but  were 
difficult  to  resolve.  However,  it  is  possible  that  they 
could  interfere  with  the  measurement  of  peptide  optical 
activity . 

The  numerical  values  of  the  a-helical  content  of 
lactoperoxidase  and  its  derivatives  should  be  viewed  with 
reservation.  The  trough  method  is  attractive  in  that  it 
provides  a  direct  detection  of  the  helical  and  random 
coiled  conformations,  but  it  neglects  the  complications 
due  to  the  possible  existence  of  6-forms  of  poly-peptides 
with  different  rotatory  characteristics  (Iizuka  and  Yang, 
1966;  Sarkar  and  Doty,  1966;  Davidson  et  al.,  1966)  and 
also  to  the  Cotton  effects  of  side-chain  chromophores  on 
the  enzyme  itself  (Vournakis  et_  al_.  ,  1968).  In  addition, 
there  may  be  forms  of  helices  other  than  the  a-helix 
which  show  different  rotatory  characteristics.  The 
rotatory  strength  of  the  3^q  helix  has  been  shown  to  differ 
from  that  of  the  a-helix  and  to  display  a  chain-length 
dependence  (Woody  and  Tinoco ,  1967).  Urry  (1967)  has 
shown  that  the  heme  chromophore  may  be  another  possible 
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source  of  optical  activity  in  the  peptide  region. 

Finally,  the  difference  in  reference  scales  for  100% 
helical  PGA  should  be  recalled. 

The  O.R.D.  curves  of  lactoperoxidase  and  its  deriva¬ 
tives  in  the  Soret  region  are  interesting  in  that 
negative  Cotton  effects  are  displayed.  It  has  been 
reported  that  cytochrome  c  peroxidase  also  shows  a 
negative  Cotton  effect  in  the  Soret  region  (Willick  et  al., 
1969).  The  Soret  O.R.D.  curves  of  LP  at  pH  3.9  and 
pH  4.5  are  identical  to  that  at  pH  7.0.  The  amplitudes 
of  the  Soret  Cotton  effects  are  as  follows:  LP ,  71°; 

LP-F ,  40°;  LP-N^  ,  56°;  and  LP-CN,  26°;  the  enzyme  concentra- 

-5 

tion  was  4.96  x 10  M  in  all  solutions.  The  reduction  of 
the  amplitude  of  the  Soret  Cotton  effect  for  the  complexes 
of  lactoperoxidase  may  indicate  some  conformational 
change  in  the  immediate  environment  of  the  heme  group 
induced  by  the  bound  ligand. 

While  the  absorption  spectrum  of  catalase  is  similar 
to  that  of  horseradish  peroxidase  and  to  spectra  of  ferri- 
hemoglobin  and  f errimyoglobin ,  the  O.R.D.  spectrum  of 
catalase  is  remarkably  different.  Horseradish  peroxidase, 
f errihemoglobin ,  and  f errimyoglobin  each  exhibit  a  positive 
Cotton  effect  in  the  Soret  region  while  in  catalase  the 
effect  is  negative.  Ullmer  and  Vallee  (1963)  have 
correlated  this  observation  with  an  orientation  of  the 
heme  group  in  catalase  different  from  the  other  heme 
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proteins  mentioned.  Marked  differences  in  the  access¬ 
ibility  of  the  heme  group  in  catalase  and  horseradish 
peroxidase  are  suggested  by  the  fact  that  the  rates  of 
catalysis  by  catalase  decrease  sharply  as  the  size  of  the 
substrate  increases,  a  phenomenon  not  observed  with  horse¬ 
radish  peroxidase.  The  O.R.D.  data  for  catalase  and 
horseradish  peroxidase  seem  pertinent  to  these  observations, 
since  they  demonstrate  that  catalase  and  horseradish 
peroxidase  provide  different  asymmetric  environments  for 
their  respective  heme  moieties.  In  this  regard,  it  is 
indicated  that  lactoperoxidase  also  has  a  heme  group  with 
a  different  accessibility  to  substrate  than  the  heme 
group  of  horseradish  peroxidase,  because  of  its  negative 
Cotton  effect.  A  steric  effect  is  also  exhibited  by  lacto¬ 
peroxidase  similar  to  catalase.  The  rate  of  reaction 
decreases  with  increasing  substrate  size  in  the  reaction 
between  lactoperoxidase  and  hydrogen  peroxide,  methyl 
hydrogen  peroxide,  and  ethyl  hydrogen  peroxide  (Chance, 


1949)  . 


f 
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Chapter  3 


The  Kinetics  of  the  Formation  of  the  Primary  Lactoperoxidase- 

Hydrogen  Peroxide  Compound 


Introduction 

Lactoperoxidase  (LP)  is  known  to  combine  with  hydrogen 
peroxide  to  form  certain  compounds  capable  of  oxidizing 
various  inorganic  and  organic  reducing  agents.  George 
(1953)  and  Yamazaki  and  Souzu  (1960)  have  suggested  that 
free  radicals  should  be  formed  from  substrates  as  inter¬ 
mediates  in  peroxidatic  oxidation.  The  proposed  mechanism 
as  applied  to  LP  is 


LP  +  H2C>2 
LP-I  +  AH2 
LP-II  +  AH2 
2  AH  • 


LP-I 

LP-II  +  AH- 

LP  +  AH- 

A  +  AH2  (or  AH-AH) 


(1) 

(2) 

(3) 

(4) 


where  LP-I  and  LP-II  are  the  oxidized  forms  of  the 
enzyme  referred  to  as  compounds  I  and  II  respectively, 

AH2  is  the  substrate  reducing  agent,  AH-  the  free  radical 
intermediate  (such  as  semiquinone  in  the  oxidation  of 
hydroquinone) ,  and  A  is  the  oxidized  product.  Chance 

(1949)  determined  the  value  of  k^  at  pH  7.0  to  be 

7-1-1 

2  x 10  M  sec  .  Our  interest  lay  in  determining  the  pH 
dependence  of  k^ .  This  chapter  describes  in  greater  detail 
the  work  already  published  (Maguire  et  al . ,  1971). 
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Application  of  the  steady-state  approximation  to  the 
reaction  scheme  leads  to  the  following  relations: 

=-k1  [LP]  [H202]  +  k3[LP-IIJ[AH2]  =  0  (5) 


d  1  =  kn  [LP  ]  [H,0,  ]  -  k,  [LP-I  ]  [AH,  ]  =  0 


2  2 


(6) 


"[LdtII]  =  k2[LP-I] [AH2]  -  k3[LP-II] [AH2]  =  0 


(7) 


=  k2  [LP-I]  [AH2]  +  k3[LP-II]  [AH2] 

The  enzyme  conservation  relation  is 

[LP]  +  [LP-I]  +  [LP-II] 

From  Equations  7  and  8, 


2k4[AH-]z  =  0  (8) 


=  [LP] 


0 


(9) 


Also 


Therefore 


[AH  *] 


d [Products ] 
dt 


pT  [LP-II] [AH  ] 


k4[AH*] 


(10) 


(ID 


d [Products ] 
dt 


k3 [LP-II] [AH2] 


(12) 


An  expression  for  [LP-II]  can  be  derived  in  terms  of 

[LP ] q  by  using  the  enzyme  conservation  relation  Eq .  9  and 

the  steady-state  Eqs .  5-7: 

kik2[LP]0[H2O2] 


(k1k2  +k3k3)  [H202]  +k2k3[AH2] 


[LP-II] 


(13) 
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Therefore 


d [H202] 

dt 


d  [Products  ] 
dt 


klk2k3  0  ^-H202^  ^  AH2  ^ 

(klk2  +klk3}  [H2°2]  +k2k3tAH2] 


(14) 


The  term  [LP] ^  refers  to  the  total  concentration  of  LP  in 
solution.  If  the  concentration  of  the  hydrogen  donor, 

AH2 ,  is  sufficiently  large  that  the  last  term  in  the 
denominator  of  Eg.  14  is  dominant,  then 


-d[H202] 

dt 


kl 0  tH2°2 1 


=  kn  ,  [H_0„] 

lobs  2  2 


(15) 


If  the  initial  concentration  of  lactoperoxidase  ([LP]q) 

is  known,  k^  may  be  obtained  from  the  relation 

k. 


ki  = 


lobs 


[LP] 


(16) 


0 


It  should  be  noted  that  Eqs .  15  and  16  are  independent  of 
any  radical  recombination  steps.  For  example,  if  reaction 
step  (4)  were 

AH*  +  AH2  - >  Products  (17) 

the  expression  for  -d[H2C>2]/dt  would  be  the  same.  There 
are  three  necessary  and  sufficient  conditions  that  allow 
the  derivation  of  our  kinetic  results.  The  first  condition 
is  that  there  should  be  one  obligatory  step  in  the  reaction 
scheme  as  in  step  (1) .  The  second  condition  is  that  there 
should  be  no  other  step  in  which  LP  or  H202  appear  as 


reactants,  and  the  third  condition  is  that  the  substrate 
concentration  should  be  in  excess  of  the  concentration  of 
hydrogen  peroxide,  and  sufficient  to  make  step  (1)  rate¬ 
determining.  The  finding  that  the  rate  is  first  order  in 
[ LP ] q  and  and  independent  of  [AH^]  or  the  nature 

of  AH^  constitutes  a  check  that  the  above  conditions  exist 
in  our  experiments. 

Experimental  Section 

Lactoperoxidase  with  a  P.N.  of  0.65  -  0.90  was  isolated 
from  cow's  milk  by  the  procedure  of  Morrison  and  Hultquist 
(1963)  and  lyophilized.  Stock  solutions  of  lactoperoxidase 
for  the  experiments  were  prepared  by  dialyzing  the  lyophil¬ 
ized  material  against  a  pH  7.0  phosphate  buffer  of  ionic 
strength  0.05.  The  stock  solutions  were  stored  at  5°  and 
were  used  within  three  days.  Lactoperoxidase  concentrations 

were  determined  from  absorbance  measurements  at  412  nm, 

5-1-1 

using  a  molar  absorptivity  of  1.14  x io  M  cm  (Morrison 
et  al . ,  1957) .  All  kinetic  measurements  were  performed  at 
25°  using  a  Cary  model  14  spectrophotometer.  Slidewires 
for  the  absorbance  ranges  0-0.2  and  0-2  were  used.  Doubly 
distilled  water  was  used  to  prepare  all  solutions.  Inorgani 
chemicals  were  of  reagent  grade  and  were  used  without 
further  purification.  Hydroquinone  was  obtained  from  the 
Baker  and  Adamson  Chemical  Company,  guaiacol  from  the  Sigma 
Chemical  Company.  Hydrogen  peroxide  concentration  determina 
tions  were  performed  according  to  the  method  of  Ovenston 
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and  Rees  (1950) .  An  Orion  model  801  digital  pH  meter  in 
conjunction  with  a  Fisher  combination  electrode  was  used 
for  all  pH  measurements. 

Steady  state  methods  were  employed  in  the  determina¬ 
tion  of  k^.  A  very  large  excess  of  reducing  agent  was 
present  in  solution,  sufficient  to  make  reactions  (2)  and 
(3)  much  faster  than  reaction  (1) .  Reaction  (4)  is  also 
presumed  much  faster  than  (1) .  Therefore  reaction  (1)  is 
rate-limiting  in  the  oxidation-reduction  cycle. 

The  rate  of  reaction  (1)  was  followed  spectrophoto- 
metrically  by  observation  of  the  rate  of  oxidation  of  an 
appropriate  hydrogen  donor.  The  scheme  depicted  in  Eq .  (1) 

to  (4)  may  be  simplified,  but  since  the  rate-limiting  step 
is  reaction  (1) ,  neglect  of  other  free  radical  intermediates 
in  the  oxidation  of  AH^  to  A  does  not  complicate  the 
predicted,  and  observed,  first-order  kinetics  for  reaction 
(1) .  It  should  be  noted  that  if  the  concentration  of 
hydrogen  donor  is  not  sufficiently  large,  the  rate  of 
appearance  of  product  will  be  governed  by  rather  complicated 
kinetics . 

Hydroquinone  was  used  as  the  hydrogen  donor  in  the  pH 
region  3.0  to  7.0;  guaiacol  from  pH  7.0  to  10.8.  Guaiacol 
was  used  at  high  pH  because  hydroquinone  apparently  under¬ 
goes  a  rapid  rate  of  oxidation  in  basic  solution  to  benzo- 
quinone  and  polymerization  products. 

When  hydroquinone  was  used  as  the  hydrogen  donor,  the 
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usual  concentrations  were:  [LP] ,  2  x 10  M;  [E^O^] , 

-5  -3 

2  x 10  M;  and  [hydroquinone] ,  1  x 10  M.  The  reaction  was 

followed  spectrophotometrically  at  290  nm.  Experimentally 

3  -1 

determined  molar  absorptivities  at  290  nm  are,  3.16  x  10  M  cm 

2  -1  -1 

for  hydroquinone  and  3  x 10  M  cm  for  benzoquinone . 

Within  experimental  error  the  observed  absorbance  decreases 
were  the  same  as  those  calculated  from  the  molar  absorp¬ 
tivities  and  the  concentration  of  hydrogen  peroxide. 

When  guaiacol  was  used,  the  usual  concentrations  were: 
[LP]  ,  2  x  10  ^  M;  f  2  x  10  ^  M ;  and  [guaiacol]  ,  10  ^  M. 

The  reaction  was  followed  spectrophotometrically  at  470  nm. 

The  product  of  the  oxidation  of  guaiacol  is  not  known  with 
certainty.  There  are  absorption  maxima  at  470  and  420  nm 
which  decay  after  some  time,  indicating  further  reaction. 

Booth  and  Saunders  (1956)  claim  that  the  color  production 
at  470  nm  is  due,  at  least  in  part,  to  the  production  of 
3 ,  3 ' -dimethoxy- [bi-2 , 5-cyclohexadienyl-l-ylidene ] -4 , 4  '  - 
dione .  The  molar  absorptivity  of  the  product  of  the  oxida¬ 
tion  is  unknown;  however,  the  absorbance  changes  during  the 
course  of  the  reaction  were  proportional  to  the  concentra¬ 
tions  of  hydrogen  peroxide  used,  in  the  range  5  x  10  ^  M  to 
4  x 10"5  M. 

The  rate  of  autoxidation  of  hydroquinone  or  guaiacol , 
negligible  in  their  respective  pH  regions,  was  compensated 
for  by  the  presence  of  the  same  hydrogen  donor  in  the 
reference  and  reaction  cuvettes.  The  hydrogen  donor 
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was  added  to  the  reaction  cuvette  first,  then  usually 
hydrogen  peroxide,  followed  by  lactoperoxidase .  The  order 
of  addition  of  hydrogen  peroxide  and  enzyme  made  no  differ¬ 
ence  to  the  total  absorbance  change  or  rate  constant  observed. 

The  problem  of  the  kinetic  results  being  vulnerable  to 
contributions  in  absorbance  (either  positive  or  negative) 
from  the  slower  secondary  reactions  (for  example,  further 
reactions  of  benzoquinone ,  perhaps  polymerization)  was 
minimal  on  the  time  scale  of  the  primary  oxidation.  For 
example,  using  the  concentrations  specified  above,  we  did 
not  notice  any  further  reaction  of  benzoquinone  for  at 
least  an  hour,  whereas  the  primary  oxidation  of  hydroquinone 
to  benzoquinone  was  over  in  less  than  half  that  time.  In 
other  words,  a  good  "infinity"  value  of  absorbance  was 
obtained.  These  remarks  also  apply  to  the  use  of  guaiacol 
as  the  reducing  agent. 

Buffers  used  in  the  order  of  increasing  pH  were :  citric 
acid-sodium  citrate,  formic  acid-sodium  formate,  acetic  acid- 
sodium  acetate,  potassium  dihydrogen  phosphate-disodium 
hydrogen  phosphate,  tris (hydroxymethyl ) aminomethane 
hydrochloride-tris (hydroxymethyl) aminomethane ,  glycine- 
sodium  glycinate,  and  sodium  bicarbonate-sodium  carbonate. 
Overlapping  buffer  systems  were  used  to  check  for  possible 
buffer  effects.  All  buffers  were  of  an  ionic  strength  of 
0.05. 

The  values  of  k,  .  were  obtained  by  analyzing  all  first 

lobs 

order  curves  with  a  non-linear  least  squares  program  on  an 
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IBM  360/67  computer.  About  fifty  points  were  taken  per 
curve  over  about  six  half-lives.  Values  of  were  then 
obtained  using  Equation  (16) .  Order  studies  (plots  of 
klobs  — *  [LP]0>  were  performed  at  pH's  3.01,  4.95,  and 

7.01  with  hydroquinone  as  the  hydrogen  donor,  and  at  pH's 
7.00  and  9.07  with  guaiacol  as  the  hydrogen  donor. 

Results 

Figure  3-1  shows  a  typical  experimental  first  order 
curve  for  the  oxidation  of  hydroquinone.  The  oxidation 
of  guaiacol  also  followed  first  order  kinetics.  Values  of 
k^  were  found  to  be  independent  of  the  P.N.  of  the  lacto- 
peroxidase.  Table  3-1  contains  values  of  k^  (with  standard 
deviations)  as  a  function  of  pH  with  hydroquinone  as  the 
reducing  agent;  Table  3-2  contains  values  of  k^  (with 
standard  deviations)  as  a  function  of  pH  with  guaiacol  as 
the  reducing  agent;  and  Table  3-3  contains  values  of  k^ 
(with  standard  deviations)  as  a  function  of  the  pH  of 
formate  buffer,  with  hydroquinone  as  the  reducing  agent. 

It  is  seen  that  guaiacol  gives  the  same  results  as  hydro¬ 
quinone  for  k^ ,  within  experimental  error,  in  the  pH  region 
at  which  kinetic  experiments  overlapped,  e.g.,  at  pH  5.0 

and  7.0.  A  concentration  effect  of  guaiacol  on  the 

-2 

observed  rate  constant  was  found  in  the  10  M  region  as 
illustrated  in  Fig.  3-2.  To  facilitate  discussion,  k| 
will  be  defined  as  that  second  order  rate  constant  which 
is  obtained  by  dividing  k^^g  by  [LP]q.  Only  in  the  limit 
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Table  3-1 

as  a  Function  of  pH  with  Hydroquinone  as  Reducing  Agent. 

Reaction  conditions  are  described  in  the  Experimental  section, 
and  is  calculated  using  Equation  (16) .  Values  of  k^ 
obtained  at  pH  4.95  and  7.01  are  results  of  order  plots 

(klobs^-  [“V* 

Number  of 


pH 

kl  ' 

(M 

-1 

sec 

k) 

★ 

Buffer 

determinations 

3.10 

(8.4 

+ 

0.7) 

X 

106 

CIT 

6 

3.37 

(7.6 

+ 

0.6) 

X 

106 

CIT 

6 

3.62 

(9.7 

+ 

0.9) 

X 

106 

CIT 

6 

3.83 

(9.9 

+ 

0 .9) 

X 

106 

CIT 

6 

3.96 

(8.7 

+ 

0.8) 

X 

io6 

CIT 

6 

4.14 

(7 . 0 

+ 

0.4) 

X 

io6 

CIT 

6 

4 . 38 

(7.6 

+ 

0.9) 

X 

io6 

CIT 

6 

4.59 

(8.8 

+ 

0.5) 

X 

io6 

CIT 

6 

3.81 

(5.7 

+ 

0.3) 

X 

io6 

AC 

6 

3.94 

(9.6 

+ 

0.5) 

X 

io6 

AC 

6 

4.14 

(1.1 

+ 

0.1) 

X 

io7 

AC 

6 

4.37 

(1.2 

± 

0.1) 

X 

io7 

AC 

6 

4 . 53 

(1.1 

+ 

0.1) 

X 

io7 

AC 

6 

4 . 77 

(9 . 9 

+ 

1.3) 

X 

io6 

AC 

6 

4.95 

(1.1 

+ 

0.1) 

X 

io7 

AC 

18 

5.16 

(1.4 

+ 

0.1) 

X 

io7 

AC 

6 

5 .36 

(1.1 

+ 

0.1) 

X 

io7 

AC 

6 

5.59 

(1.5 

+ 

0.1) 

X 

io7 

AC 

6 

5.78 

(1.4 

+ 

0.1) 

X 

io7 

AC 

6 

5.96 

(1.1 

+ 

0.2) 

X 

io7 

P 

10 

6.15 

(1.2 

+ 

0.1) 

X 

io7 

P 

10 

6.38 

(1.0 

+ 

0.1) 

X 

io7 

P 

10 

6.59 

(1.2 

+ 

0.1) 

X 

io7 

P 

10 

6.80 

(1.1 

+ 

0.1) 

X 

io7 

P 

10 

continued 


< 


Table  3-1  continued 


Number  of 


,  ,  -1 

■  1 

* 

PH 

k^  (M  sec 

) 

Buffer 

determinations 

7.01 

(9.4  ±0.8)  x 

io6 

P 

52 

7.17 

(1.2  ±  0.1)  x 

io7 

P 

10 

7.37 

(1.0  ±  0.1)  x 

io7 

P 

10 

Buffer  key:  CIT,  citrate;  AC,  acetate;  P,  phosphate. 
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Table  3-2 

k^  as  a  Function  of  pH  with  Guaiacol  as  Reducing  Agent. 

Reaction  conditions  are  described  in  the  Experimental  section, 
and  is  calculated  using  Equation  (16) .  Values  of  k^ 


obtained  at 

pH  7  . 

0  0  and 

9 . 07 

are  results 

of  order  plots 

(klobsVS-  tLP]0>- 

pH 

ki 

,  -1  -1, 

(M  sec  ) 

* 

Buffer 

Number  of 

determinations 

4 .97 

(1.0 

±0.1) 

x  107 

AC 

6 

7.00 

(1.1 

±0.1) 

x  107 

P 

15 

8.07 

(1.0 

±0.1) 

X  107 

T 

6 

9.04 

(1.0 

±  0.1) 

x  107 

G 

6 

9 .07 

(1.0 

±0.1) 

X  107 

T 

15 

9.93 

(9.9 

±0.2) 

x  106 

G 

6 

10.02 

(9.8 

±0.2) 

x  106 

CAR 

6 

10.76 

(9.7 

±0.3) 

x  106 

CAR 

12 

* 

Buffer  key 

:  AC, 

acetate;  P 

,  phosphate; 

T,  tris; 

G,  glycinate;  CAR,  carbonate. 
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Table  3-3 

as  a  Function  of  pH  of  Formate  Buffer  with  Hydroquinone 
as  Reducing  Agent. 

Reaction  conditions  are  described  in  the  Experimental 
section,  and  is  calculated  using  Equation  (16) .  The 
value  of  k^  at  pH  3.01  was  determined  by  an  order  plot 

(klobs^-  [LP  V  • 

Number  of 


PH 

kl 

,  -1 
(M  sec 

h 

determinations 

3.01 

(7.7 

±1.0) 

X 

104 

14 

3.18 

(1.2 

±0.1) 

X 

io5 

6 

3.27 

(1.9 

±0.1) 

X 

105 

6 

3.40 

(1-9 

±0.1) 

X 

105 

6 

3.50 

(3.0 

±0.1) 

X 

105 

6 

3.60 

(2.9 

±0.1) 

X 

105 

6 

3.74 

(4.7 

±0.3) 

X 

io5 

6 

3.79 

(4.2 

±0.1) 

X 

io5 

6 

4 .03 

(6.4 

±0.4) 

X 

io5 

6 

4.25 

(9.9 

±0.3) 

X 

105 

6 

4.39 

(1.2 

±0.1) 

X 

io6 

6 

4.67 

(1.7 

±0.1) 

X 

io6 

3 

4.96 

(3.4 

±0.1) 

X 

io6 

6 

% 
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of  sufficiently  high  concentration  of  reducing  agent  is 

I 

k,  ,  a  true  first  order  rate  constant,  and  k,  equal  to 

k-^ ,  the  true  second  order  rate  constant  for  compound  I 

formation  (Eq.  16) .  Figure  3-2  shows  that  at  pH  7.00,  the 

value  of  k^  increases  with  increasing  guaiacol  concentra- 

-2 

tion  until  the  latter  reaches  about  5  x  10  M;  the  value 

of  k|  then  becomes  independent  of  guaiacol  concentration 

to  10  "*■  M  and  hence  equal  to  k^  .  At  pH  10.02,  k^  becomes 

independent  of  guaiacol  concentration  (hence  equal  to  k-^) 

-2 

when  the  latter  reaches  about  7  x 10  M  as  shown  in 
Fig.  3-3.  At  pH  10.76,  where  guaiacol  is  ionized,  a 
number  of  experiments  were  performed  in  which  the  ionic 
strength  was  varied  from  0.1  to  0.3  and  guaiacol  concentra¬ 
tions  were  about  0.1  M.  No  ionic  strength  effect  on  the 
rate  constant  k^  was  observed. 

Figure  3-4  is  a  plot  of  vs. [LP]q  at  pH  7.01. 

From  the  slope,  k^  is  calculated  by  least  squares  analysis 

6  “1 

to  be  (9.4  ±  0.8)  x  10  M  sec  .  The  intercept  is  zero 
within  experimental  error,  showing  first  order  behaviour; 
this  was  the  case  for  all  such  plots  of  k^^g  vs.  [ LP ] q 
at  other  pH  values. 

A  plot  of  k^  vs.  pH  over  the  pH  range  3.0  to  10.8  is 
shown  in  Fig.  3-5.  The  filled-in  circles  represent  values 
of  k^  obtained  when  experiments  were  performed  in  formate 
buffer.  The  line  through  the  rest  of  the  points  was 
obtained  by  least-squares  analysis.  The  slope  of  the  line 
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[lp]  *io9m 


Fig.  3-4.  A  plot  of  k^Q,  vs .  [LP]  at  pH  7.01  using  hydro- 
quinone  as  reducing  agent.  This  plot  is  typical  of  the 
five  order  plots  performed  at  the  pH's  described  in  the 
Experimental  section.  The  slope,  equal  to  k, ,  has  a  value 

r  —  i  _ -|  -L 

of  (9.4  ±0.8)  x 10  M  s  ,  calculated  by  least-squares 
analysis . 
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Fig.  3.5.  A  plot  of  k-^  vs.  pH  from  pH  3.0-10.8.  The  filled-in  circles  represent  values 
of  obtained  in  formate  buffer.  The  values  of  k^  are  those  in  Tables  3-1 , 3-2  and  3-3.  The 
slope  of  the  line  through  the  open  circles  is  zero  within  experimental  error,  (1.1  ±1.2)  xl 
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is  zero  within  experimental  error,  indicating  a  constant 

6  —  1  —  1 

value  of  of  (9.2  ± 0.9)  x  10  M  sec  .  In  formate  buffer 
however,  the  value  of  decreases  from  3.4  x  10^  M  ^sec  ^ 
at  pH  4.96  to  7.7  x  10^  M  ^sec  ^  at  pH  3.01. 

It  was  necessary  to  store  the  lactoperoxidase  in  stock 
solutions  at  pH  7.0.  It  was  noticed  that  if  the  LP  v/ere 
stored  at  pH  3.0,  for  example,  the  values  of  k1  progressively 
decreased  with  time,  indicating  deactivation  of  the  lacto¬ 
peroxidase  at  this  extreme  of  pH. 

Discussion 

It  is  seen  from  Fig.  3-5  that  the  rate  constant  for 
the  formation  of  lactoperoxidase  compound  I  from  lacto¬ 
peroxidase  and  hydrogen  peroxide  is  essentially  pH-independ- 
ent  over  the  pH  region  investigated.  The  value  of 

(9.2  ±0.9)  x 10°  M  sec  is  to  be  compared  with  Chance's 

7  -1  -1 

(1949,  1950)  value  for  ^  of  2  x  10  M  sec  ,  slightly  more 
than  twice  the  value  reported  here.  One  possible  reason  for 
the  discrepancy  may  be  that  the  value  for  the  extinction 
coefficient  at  412  nm  is  incorrect  by  a  factor  of  two.  This 
would  cause  the  concentrations  of  LP  to  be  incorrect  by  a 
factor  of  two,  giving  rise  to  a  rate  constant  k^  that  is 
different  by  a  factor  of  two.  However,  although  the 
extinction  coefficient  may  be  in  error,  it  would  appear 
unlikely  that  it  would  be  in  error  by  this  amount. 

Another  possibility,  that  the  enzyme  cycle  could  be 
affected  by  the  attack  of  a  free  radical  on,  for  example, 

LP ,  thereby  reducing  the  amount  of  LP  in  the  cycle,  could 
be  discounted  on  the  grounds  that  first  order  kinetics 


are  observed  throughout  the  course  of  a  particular  reaction. 

There  could  be  some  unforeseen  source  of  error  in  our 

results,  but  to  the  best  of  our  knowledge,  the  results  yield 

accurate  values  of  k^.  On  the  other  hand,  Chance  (1949, 

1950)  gives  little  experimental  detail  concerning  the  condi- 

7-1-1 

tions  under  which  his  value  of  2  x 10  M  sec  was  obtained. 
No  mention  is  made  of  the  number  of  determinations  or  the 
reproducibility  of  the  value  of  k1  obtained;  Chance  (1950) 
does  mention  obtaining  values  of  k^  as  low  as  9  x 10^  M  ^sec 
and  as  high  as  2.7  x 107  M"1sec_1.  There  is  also  no  mention 
made  of  ionic  strength  or  temperature,  although  Morell 
(1954)  claims  that  Chance's  experiments  were  performed  at 
temperatures  varying  from  25°  to  30°.  For  these  reasons,  it 
is  concluded  that  the  value  of  ^  reported  here  is  more 
reliable . 

There  is  some  evidence  that  k^  and/or  k^  of  the  reac¬ 
tion  scheme  Eqs.  1-4  are  not  pH-independent ,  as  is  k][. 

This  is  implied  by  the  k^  vs_.  [guaiacol]  plots  at  pH  7.00 
(Fig.  3-2)  and  pH  10.02  (Fig.  3-3).  The  concentration 
region  in  which  k^'  becomes  independent  of  guaiacol  concen¬ 
tration  (hence  equal  to  k^  occurs  at  higher  guaiacol 
concentration  at  pH  10.02  than  at  pH  7.00,  as  shown  in  the 
Results  section.  This  indicates  that  it  is  more  difficult 
to  attain  the  conditions  under  which  reaction  (1)  is  rate- 
limiting  at  this  higher  pH,  i.e.,  the  concentration  of 
guaiacol  must  be  increased  in  order  to  maintain  pseudo- 
first  order  conditions.  This  could  be  explained  by  a 
decrease  in  some  other  rate  constant  (s)  in  the  cycle  with 
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increasing  pH.  In  fact,  it  will  be  shown  in  the  next  two 
chapters  that  does  decrease  with  increasing  pH  for  the 
reactions  between  LP-II  and  either  iodide  ion  or  p-cresol. 

Figure  3-5  also  shows  the  effect  of  formate  buffer 
on  the  kinetics  of  the  formation  of  compound  I.  It  was 
determined  that  in  pH  3.60  formate  solution  there  is  no 
shift  in  the  lactoperoxidase  Soret  band  maximum,  but  there 
is  a  7%  decrease  in  the  molar  absorptivity  compared  to 
standards  either  in  pure  water  or  pH  7.0  phosphate  solution. 
Fig.  3-6  and  3-7  are  plots  of  k^  vs .  [hydroquinone]  obtained 
from  experiments  performed  in  pH  3.6  acetate  and  pH  3.6 
formate  solutions.  In  pH  3.6  acetate  buffer,  k^  becomes 

independent  of  hydroquinone  concentration  (hence  equal  to 

-4 

k^)  at  about  4  x 10  M,  while  in  pH  3.6  formate  buffer  this 
occurs  at  a  higher  hydroquinone  concentration  of  about 
5  x  10-4  M. 

More  striking  is  the  vertical  difference  between  the 

two  plateaus  of  Fig.  3-6  and  Fig.  3-7.  The  value  of  k^ 

(i.e.,  the  limiting  value  of  k^)  in  pH  3.6  formate  buffer 
5-1-1 

is  3.5  x 10  M  sec  while  the  value  of  k^  in  pH  3.6 
acetate  buffer  is  7.5  x 106  M  ^sec  1,  a  factor  of  20 
higher.  This  difference  can  be  thought  of  as  a  measure  of 
the  extent  of  binding  of  formate  ion  to  lactoperoxidase, 
and  can  be  used  to  calculate  a  dissociation  constant  for 
the  lactoperoxidase-formate  ion  complex.  In  the  equation 
K  =  [LP] [HCO~]/ [LP-HCO“] , [HCO“]  =  0.05  (the  ionic  strength 
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of  all  buffers  in  this  study) ,  and  the  ratio  [LP ] / [LP-HCO" ] 
was  calculated  from  the  ratio  of  rate  constants 
obtained  in  formate  buffer  and  another  buffer  at  the  same 
pH  (such  as  acetate  or  citrate,  again  at  the  same  ionic 

strength  of  0.05).  The  dissociation  constant  thus  calcula- 

-4  -2 

ted  varie s  from  4.2x10  M  at  pH  3.01  to  1.85x10  Mat 

pH  4.96.  Fig.  3-8  is  a  plot  of  the  negative  logarithm  of 

the  dissociation  constant  (pK)  vs_.  pH.  The  slope,  obtained 

by  least  squares  analysis,  is  -(0.77  ±0.03).  A  value  of 

less  than  unity  for  the  slope  may  indicate  that  there  is 

a  pK  of  the  free  enzyme  in  this  region.  Chance  (1952a)  has 

demonstrated  that  the  binding  of  formate  buffer  to  catalase 

involves  only  a  small  perturbation  of  the  Soret  spectrum, 

and  also  that  the  dissociation  constant  of  the  catalase- 

-4 

formate  ion  complex  varies  from  5  x 10  M  at  pH  5.6  to 

-  5 

1.6  x 10  M  at  pH  4.1.  Chance  (1952a)  further  demonstrates 
that  a  plot  of  the  negative  logarithm  of  the  dissociation 
constant  (for  the  catalase  formate-ion  complex)  vs.  pH 
is  linear  with  a  slope  of  -1,  indicating  that  for  catalase 
there  is  no  ionizing  group  on  the  enzyme  which  affects  the 
binding  reaction. 
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Chapter  4 

The  Kinetics  of  the  Oxidation  of  Iodide  Ion 
By  Lactoperoxidase  Compound  II 

Introduction 

The  horseradish  peroxidase-catalyzed  oxidation  of 
iodide  ion  was  first  observed  by  Bach  (1904,  1907)  at  the 
turn  of  the  century.  In  the  thyroid  gland  a  peroxidase 
with  physical  properties  similar  to  those  of  the  horse¬ 
radish  enzyme  oxidizes  iodide  ion  as  a  step  in  thyroid 
hormone  biosynthesis  (Hosoya  and  Morrison,  1967a). 

Hosoya  and  Morrison  (1967b)  have  demonstrated  that  lacto¬ 
peroxidase  is  more  active  in  catalyzing  the  oxidation  of 
iodide  ion  than  is  thyroid  peroxidase,  myeloperoxidase, 
or  horseradish  peroxidase.  The  chloroperoxidase  of 
Caldariomyces  fumago  plays  a  similar  role  in  the  synthesis 
of  caldariomycin ,  which  contains  chlorine  (Hager  et  al . , 
1966) .  Hence  the  peroxidase-catalyzed  oxidation  of  halide 
ions  is  of  particular  interest.  The  present  study  was 
undertaken  for  two  reasons.  One  reason  was  that  a  study  of 
this  sort  would  be  preliminary  to  a  study  of  a  system  in 
which  the  kinetics  of  the  iodination  of  tyrosine  would  be 
studied,  thereby  providing  a  close  analogy  to  the  action  of 
thyroid  peroxidase  in  the  thyroid  gland.  Another  reason 
this  study  was  undertaken  was  to  probe  the  active  site  of 
compound  II  of  lactoperoxidase  in  the  hope  that  the  log 
k-pH  profile  might  provide  clues  as  to  the  identity  of 
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amino  acid  groups  on  the  enzyme  which  participate  in  the 
reaction.  This  chapter  describes  in  greater  detail  the 
work  in  press  (Maguire  and  Dunford,  1972)  . 

Experimental  Section 

Lactoperoxidase  was  obtained  as  a  lyophilized  powder 

from  Calbiochem.  In  this  form,  the  samples  exhibited  a 

P.N.1  of  about  0.6.  The  enzyme  was  purified  by  gel 

filtration  on  Sephadex  G-200  Superfine  at  4  °  using  a 

phosphate  buffer  of  pH  7  and  ionic  strength  0.05  as  eluant. 

Enzyme  fractions  obtained  in  this  way  which  exhibited 

P.N.'s  greater  than  0.8  were  used  in  this  study.  The 

concentration  of  lactoperoxidase  was  determined  spectro- 

photometrically  at  412  nm  using  a  molar  absorptivity  of 
5  —  1  —  1 

1.14  xlO  M  cm  (Morrison  et  al . ,  1957). 


Abbreviations  used  are:  LP ,  lactoperoxidase;  LP-I ,  the 
primary  lactoperoxidase-hydrogen  peroxide  compound;  LP-II, 
the  secondary  lactoperoxidase-hydrogen  peroxide  compound; 
HRP-II,  the  secondary  horseradish  peroxidase-hydrogen 
peroxide  compound;  P.N.,  purity  number,  the  ratio  of  the 
absorbance  of  a  solution  of  lactoperoxidase  at  412  nm  to 
that  at  280  nm;  y,  ionic  strength;  kQbs ,  the  pseudo-first 
order  rate  constant  for  the  reaction  of  LP-II  with  iodide 
ion;  kn  ,  k„ :  second  and  third  order  rate  constants,  respect- 
ively,  for  the  reaction  of  LP-II  with  iodide  ion;  [  ]  ,  molar 
concentration;  A  V,  A  A:  changes  in  voltage  and  absorbance, 
respectively . 
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The  buffers  used  in  this  study,  prepared  from  reagent 
grade  chemicals  and  of  an  ionic  strength  of  0.05,  are 
listed  in  Table  4-1.  An  Orion  model  801  digital  pH  meter 
in  conjunction  with  a  Fisher  combination  electrode  was 
used  for  all  pH  measurements. 

Potassium  iodide  was  obtained  from  Alfa  Inorganics 
and  the  McArthur  Chemical  Company,  and  sodium  iodide  from 
Orion  Research,  Inc.  The  three  iodide  samples  were  found 
to  exhibit  identical  chemical  behavior  within  experimental 
error . 

Solutions  of  hydrogen  peroxide  were  prepared  by 
diluting  a  30%  solution  of  hydrogen  peroxide  obtained  from 
the  Fisher  Scientific  Company.  Concentrations  of  hydrogen 
peroxide  were  determined  according  to  the  method  of  Ovenston 
and  Rees  (1950)  .  The  water  used  in  all  solutions  was 
distilled  once  from  alkaline  permanganate  and  once  from 
glass . 

The  rate  of  reaction  of  LP-II  with  iodide  ion  was 
followed  spectrophotometrically  at  412  nm,  the  wave-length 
of  maximum  absorbance  by  native  lactoperoxidase .  The 
concentration  of  iodide  ion  was  at  least  ten  times  that  of 
LP ,  maintaining  pseudo-first  order  conditions.  LP-II  was 
prepared  in  most  cases  by  the  addition  of  one  equivalent 
of  hydrogen  peroxide  to  a  LP  solution,  and  had  a  half-life 
of  over  20  minutes  at  most  pH's.  In  tris-HCl  buffer, 
about  1.2  to  1.5  equivalents  of  hydrogen  peroxide  were 
required,  which  may  be  due  to  the  presence  of  some  reducing 


' 
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agent  in  this  buffer.  LP  and  LP-II  were  found  to  exhibit 
an  isosbestic  point  at  422  nm,  and  LP-II  an  absorption 
maximum  at  432  nm,  in  agreement  with  the  spectra  obtained 
by  Chance  (1950)  . 

The  reactions  between  LP-II  and  iodide  ion  at  25°  were 
studied  using  a  Cary  14  recording  spectrophotometer  (equip¬ 
ped  with  0-0.2  and  0-2  absorbance  slidewires) ,  and  a 
stopped-flow  apparatus  constructed  in  this  laboratory 
which  has  a  dead-time  of  6  milliseconds.  This  stopped- 
flow  apparatus  is  nearly  identical  to  that  described  by 
Hasinoff  (1970)  except  that  the  cell  compartment  is 
constructed  of  polypropylene  instead  of  plexiglass.  The 
stopped-flow  data  were  obtained  as  an  amplified  photo¬ 
multiplier  voltage  print-out  digitalized  at  30  equally 
spaced  intervals  of  time.  The  optical  detection  system 
has  been  described  elsewhere  (Ellis  and  Dunford,  1968)  . 

In  stopped-flow  experiments,  one  driving  syringe 
contained  the  LP-II  in  unbuffered  aqueous  solution,  and 
the  other  syringe  contained  iodide  ion  in  the  appropriate 
buffer.  This  procedure  minimized  the  possibility  of 
enzyme  denaturation  at  the  low  pH’s  at  which  some  experi¬ 
ments  were  performed.  In  addition  the  total  absorbance 
change  was  kept  small  (A A  <  0.04)  so  that  the  relative 
voltage  changes  observed,  AV,  were  proportional  to  A  A. 

The  kinetic  data  were  analyzed  as  described  by 
Roman  et  al .  (1971) .  At  each  pH  and  for  a  given 
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concentration  of  iodide  ion,  3-5  experiments  were  performed 
on  the  Cary  spectrophotometer,  or  about  10  experiments 
were  performed  if  the  stopped-flow  apparatus  was  used. 

When  it  was  possible,  concentrations  of  iodide  ion  were 
used  which  would  yield  reaction  rates  which  overlapped  the 
time  scales  pertaining  to  the  Cary  spectrophotometer  and 
the  stopped-flow  apparatus;  using  both  instruments, 
experiments  at  the  same  concentration  of  iodide  ion  (and 
the  same  pH)  yielded  values  for  a  pseudo-first  order  rate 
constant,  kQ^s ,  which  were  equal  within  experimental  error. 

Blank  reaction  rate  experiments  were  conducted  at 
various  pH's  in  the  absence  of  LP ,  and  in  no  case  was  an 
appreciable  rate  measured  (by  observing  absorbance  changes 
at  353  nm,  an  absorbance  maximum  of  tri-iodide  ion)  for 
the  uncatalyzed  reaction  between  iodide  ion  and  hydrogen 
peroxide . 

The  products  produced  in  the  reaction  between  LP-II 
and  iodide  ion  (iodine  atoms  or  molecular  iodine)  may 
react  further  with  lactoperoxidase ,  perhaps  with  aromatic 
amino  acid  residues.  It  was  noticed  in  experiments 
performed  using  the  stopped-flow  apparatus  that  there  exists 
a  second,  much  slower,  reaction.  At  pH  5,  this  reaction 
is  50  times  slower  than  the  reaction  between  LP-II  and 
iodide  ion,  and  so  caused  no  detectable  interference  in 
measurements  of  the  rate  of  the  reaction  of  LP-II  with 
iodide.  Calculation  of  rate  constants  using  the  initial 
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time  course  of  the  LP-II — iodide  reaction  led  to  values 
which  were  slightly  higher  (ca.  3%)  than  those  rate 
constants  obtained  by  computer  analysis  of  the  whole  first 
order  curve.  Because  this  value  is  within  experimental 
error,  it  was  not  necessary  to  apply  the  correction,  and 
it  was  concluded  that  this  second,  unidentified,  reaction 
did  not  interfere  with  observation  of  the  reaction  between 
LP-II  and  iodide  ion. 

Results 

The  kinetics  of  the  reaction  between  LP-II  and  iodide 
under  conditions  of  high  iodide  ion  concentration  in 
relation  to  the  concentration  of  LP-II  are  consistent  with 
the  differential  rate  expression 

— d  [ LP.~Hl  =  k  ,  [LP-II]  (1) 

obs L  J 

dt 

Linear  semilogarithmic  plots  (over  four  or  more  half-lives) 
of  AV  (or  AA)  vs.  time  proved  the  validity  of  Equation  1. 

An  example  of  an  experimental  curve  is  shown  in  Fig.  4-1; 
Fig.  4-2  shows  a  semilogarithmic  plot  of  the  curve  in 
Fig.  4-1,  and  also  plots  of  curves  obtained  at  two  other 
pH  values. 

Fig.  4-3  is  a  plot  of  kQbs  vs.  the  concentration  of 
iodide  ion,  using  data  obtained  at  pH  6.95.  The  non-linear¬ 
ity  is  evidence  that  a  higher  order  in  concentration  of 
iodide  ion  is  involved.  Values  of  k^^g  as  a  function  of 
the  concentration  of  iodide  ion  were  found  to  fit  the 


relation 
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Fig.  4-3. 
|i  =  0.05. 
simply  fi 


A  plot  of  kQbs  vs.  [I  J  at  pH  6.95  phosphate. 
The  curvature  shows  that  the  reaction  is  not 
st  order  in  the  concentration  of  iodide  ion. 
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kobs  =  kitI_]  +  k2[1-]2  <2) 

where  and  are  second  and  third  order  rate  constants 
respectively.  A  plot  of  kobs/[I“]  vs.  [I"]  yields  k2 
from  the  slope  and  k^  from  the  intercept.  A  plot  of 
^obs^1-]  vising  data  obtained  at  pH  6.95  is  shown 

in  Fig.  4-4.  The  plot  curves  up  at  low  concentrations  of 
iodide  ion  because  of  the  spontaneous  decay  of  LP-II.  The 
values  of  the  rate  constants  k^  and  k2  were  determined  by 
least  squares  analysis.  The  rate  constant  k^  varies  from 
(2.6  ±0.4)  xio"1  M_1sec_1  to  (4.20  +0.44)  x  106  M_1sec-1 
over  the  pH  range  10.07  -  2.86,  while  the  rate  constant  k2 
varies  from  (1.41  ±0.17)  x 103  M-2sec_1  to  (2.72  ±0.24)  x 
10^^  M  2sec  ^  over  the  same  pH  interval. 

The  rate  constants  k^  and  k2  are  plotted  logarithmic¬ 
ally  vs.  pH  in  Fig.  4-5  and  4-6,  and  listed  in  Table  4-1. 

The  slope  (obtained  by  least  squares  analysis)  of  the  log 
k2  vs.  pH  plot  is  -1.03  ±0.10;  the  slope  of  the  log  k^  vs. 
pH  plot  is  -0.98  ±0.04.  The  error  in  k^  is  greater  than 
the  error  in  k2  because  of  the  interference  of  the  spontan¬ 
eous  rate  of  decay  of  LP-II  at  the  low  concentrations  of 
iodide  ion  in  plots  such  as  Fig.  4-4. 

A  slight  perturbation  of  the  Soret  spectrum  of  lacto- 
peroxidase  was  observed  when  the  concentration  of  iodide  ion 
was  in  very  large  excess  over  that  of  lactoperoxidase . 

Fig.  4-7  shows  a  difference  spectrum  at  pH  7  of  LP  +  I" 
vs.  LP.  The  same  difference  spectrum  is  obtained  at  pH  4.6. 
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Fig.  4-5.  A  semilogarithmic  plot  of  vs.  pH  for  the 
reaction  of  LP— II  with  iodide  ion .  The  error  bars  on  the 
points  are  standard  deviation s  obtained  by  analysis  of 
plots  such  as  Fig.  4—4  at  each  pH.  The  solid  line  is  a 
least  squares  fit,  with  a  slope  of  -  (0 . 98  ±  0 . 04 )  . 
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Fig.  4-6.  A  semilogarithmic  plot  of  k2  vs.  pH  for  the 
reaction  of  LP-II  with  iodide  ion.  The  error  bars  on  the 
points  are  standard  deviations  obtained  by  analysis  of 
plots  such  as  Fig.  4-4  at  each  pH.  The  solid  line  is  a 
least  squares  fit,  with  a  slope  of  -(1.03+0.10). 
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There  is  a  maximum  in  the  difference  spectrum  at  431  nm, 
and  a  minimum  at  410  nm,  an  observation  made  previously 
by  Morrison  et_  al .  (1970)  .  It  is  not  clear  if  there  is 

binding  of  iodide  ion  to  the  iron  of  lactoperoxidase . 

Perhaps  the  perturbation  of  the  Soret  spectrum  is  due  to 
a  structural  alteration  induced  by  the  high  concentration 
of  salt  in  the  enzyme  solution;  we  observed  that  nitrate 
ion  at  the  same  concentration  produced  the  same  sort  of 
difference  spectrum.  In  addition,  we  were  unable  to 
remove  cyanide  ion  bound  to  LP  (by  observing  a  shift  of  the 
Soret  maximum  of  LP-CN  at  432  nm  back  to  412-413  nm)  by 
adding  iodide  ion  up  to  very  high  concentrations  (ca.  0.25  M) . 
In  general  hemoprotein  halide  complexes  form  more  readily  at 
low  pH;  however,  our  spectral  results  give  no  indication  of 
such  a  trend  for  iodide  and  indeed  provide  little  evidence 
for  any  binding  of  iodide.  It  can  at  least  be  said  that 
if  iodide  ion  binds  to  the  iron  of  lactoperoxidase,  the 
extent  of  binding  is  much  less  than  that  of  fluoride  (Segal 
et  al . ,  1968)  and  cyanide  (Dolman  et  ad. ,  1968)  . 

Discussion 

The  plots  of  log  vs.  pH  and  log  k vs.  pH  are  linear 
with  a  slope  of  -1  within  experimental  error  over  the  pH 
range  of  the  study.  This  behavior  can  be  explained  by  a 
kinetically  important  ionization  on  either  the  enzyme  or 
substrate  (Dixon  and  Webb,  1964)  which  is  outside  the  pH 
range  of  this  study  to  low  pH.  The  pK^  of  hydriodic  acid 


! 
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has  been  estimated  to  be  about  -9  (Bell,  1959).  Using  this 
value,  it  can  be  shown  that  the  rate  constant  for  the 
reaction  of  LP-II  with  HI  would  exceed  the  maximum  value 
of  the  diffusion-controlled  limit  (10^ J  M  ^sec  by  about 
eight  orders  of  magnitude .  For  example,  if  the  system  were 


LP-II 


LP-II 


+  HI 
/ 

kl 

+  I 


4 


Products 


Then  d [Products] 
dt 


ka[LP-II] [HI]  +kb[LP-II][I  ] 


k  [LP-II] [I  ] total  k,  [LP-II] [I  ] total 
a  -  +  -2 - 


K, 


1  + 


[H+] 


1  + 


isii 

K„ 


.  (3) 


which  is  to  be  compared  with  the  experimentally  determined 
equation 


d[Products]  .  fTT^  T_,r_-,.  ,  , 

-d - — - L  =  kn  [LP-II]  [I  ] total 


(4) 


Therefore 


kl  = 


K, 


1  + 


[H+] 


1  + 

kl 


(5) 


If  we  suppose  that  HI  is  the  reactive  species,  then  kj 
and 

k 

k,  =  - ^ 


=  0, 


K 


(6) 


1  + 


[H+] 


6  “l  “1 

At  pH  3,  k^  is  about  10  M  sec  and  therefore 


k  =  106  (l  + 

a  V  10"3 


=  10 


18 


(7) 
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which  clearly  exceeds  the  diffusion-controlled  limit.  The 
same  is  true  at  any  pH  value.  Thus  the  kinetic  results  can 
be  explained  if  LP-II  exists  in  two  forms  with  the  protonated 
form  reacting  much  more  rapidly  than  the  unprotonated  form 
with  iodide  ion. 

A  simple  mechanism  which  incorporates  the  effects  of 
first  and  second  order  concentrations  of  iodide  ion  (but  not 
the  ionization  in  the  enzyme)  is: 


where  P  and  P  are  the  products. 

LP-II-I  is  a  complex  formed  between  LP-II  and  iodide 
ion,  and  K  is  its  dissociation  constant.  One  may  be  able 
to  make  some  choice  as  to  which  steps  are  more  important 
in  the  mechanism.  From  the  work  of  Roman  et  al. ,  (1971)  it 

is  known  that  HRP-II  does  not  form  a  long-lived  complex  with 
iodide  ion.  However,  the  second  order  rate  constant  k^  for 
the  reaction  of  LP-II  with  iodide  is  larger  by  a  factor  of 
12  over  the  comparable  constant  for  HRP-II.  This  indicates 
greater  efficiency  in  the  LP-II  oxidation  and  therefore 
favors  the  formation  of  the  LP-II-iodide  complex.  The 
direct  reaction,  with  rate  constant  k  ,  cannot  be  excluded 
on  the  basis  of  a  purely  kinetic  argument  however.  The 
reaction  which  is  second  order  in  iodide  points  clearly  to 
the  importance  of  a  complex  between  LP-II  and  iodide  ion. 
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The  probability  of  termolecular  collisions  would  appear  too 
small  to  account  for  a  significant  proportion  of  the  total 
rate;  this  possibility  has  not  been  shown  in  the  mechanism. 
The  alternative  involves  the  reaction  between  iodide  ion 
and  the  LP-II-iodide  complex,  an  example  of  substrate 
activation . 

It  is  possible  to  place  limits  on  the  value  of  the 

dissociation  constant,  K  ,  of  the  kinetically  important 

ionizing  group  on  the  enzyme.  From  the  lack  of  curvature 

at  low  pH  of  the  log  k^  vs.  pH  plot,  an  upper  limit  can  be 

assigned  of  pK  <2.8.  A  lower  limit  of  pK  is  obtained  by 

a  a 

extrapolation  of  the  log  k^  vs.  pH  plot  to  the  diffusion- 

g 

controlled  limit.  Depending  upon  whether  one  uses  10  or 

1010  M  ^sec  for  this  diffusion-controlled  rate  (Alberty 

and  Hammes,  1958)  one  obtains  pK  >1.2  or  pK  ^-0.8.  This 

a  a 

pK^  value  lies  outside  the  range  of  ionization  constants 
for  the  common  amino  acids  found  in  proteins,  and  may  be 
due  to  the  influence  of  the  heme  group  on  the  ionization  of 
an  acid,  or  the  ionizing  group  may  be  a  component  of  the 
heme  group.  An  acid  group  bound  in  the  fifth  or  sixth 
coordination  position  of  the  heme  iron  appears  an  attractive 
possibility,  but  an  acid  group  obtained  from  protonation 
of  the  porphyrin  ring  cannot  be  excluded. 

The  reaction  of  LP-II  with  iodide  is  similar  to  the 
HRP-II -iodide  reaction  in  that  the  catalytic  importance  of 
a  single  acid  group  on  the  enzyme  is  clearly  demonstrated 
in  both  reactions  (Roman  et  al . ,  1971).  However ,  the 


faster  rate  and  the  importance  of  a  term  second-order  in 
iodide  in  the  LP-II  reaction,  point  clearly  to  the  pitfall 
if  one  attempts  to  extrapolate  results  from  one  enzymatic 
reaction  to  the  analogous  reaction  for  a  related  enzyme. 
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Chapter  5 

The  Kinetics  of  the  Oxidation  of  p-Cresol  by 

Lactoperoxidase  Compound  II 

Introduction 

As  stated  in  previous  chapters,  lactoperoxidase 
catalyzes  the  oxidation  of  a  wide  variety  of  compounds 
by  hydrogen  peroxide  according  to  the  following 
generally  agreed  upon  reaction  scheme: 


LP 

+  H2°2 

- > 

LP 

-I 

LP-I 

+  ah2 

- > 

LP 

-II 

+ 

LP-II 

+  AH2 

\ 

/ 

LP 

+ 

AH 

2AH  •  - 

- > 

Products 

where  LP  represents  the  native  enzyme,  LP-I  and  LP-II  are 
the  oxidized  forms  of  the  enzyme  referred  to  as  compounds 

I  and  II,  respectively,  and  AH^  is  the  oxidizable  substrate. 

Previous  work  on  the  kinetics  of  reactions  of  compound 

II  has  involved  the  oxidation  of  iodide  ion  (Maguire  and 
Dunford,  1972),  described  in  Chapter  4.  The  present  study 
was  undertaken  with  a  view  to  determining  the  effect  of  a 
neutral  molecule  on  the  reduction  of  compound  II.  p-Cresol 
was  chosen  because  of  the  remarkable  efficiency  of  oxidation 
of  phenols  by  peroxidase  systems  (Saunders  ejt  al .  ,  1964)  , 
and  because  of  its  resemblance  to  the  amino  acid  tyrosine, 
while  having  fewer  acid-base  groups.  A  peroxidase  system 
has  been  suggested  to  be  active  in  the  thyroid  gland 
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(Hosoya  and  Morrison,  1967a),  which  produces  the  hormone 
thyroxine  by  coupling  iodinated  tyrosine  residues.  This 
chapter  describes  in  greater  detail  material  submitted  for 
publication . 

Experimental 

Lactoperoxidase  (LP)  v/as  obtained  from  Calbiochem  as 
a  lyophilized  powder,  and  purified  by  gel  filtration  at  4° 
on  a  Sephadex  G-200  Superfine  column  with  a  pH  7.0  phosphate 
buffer  as  eluant.  The  ratio  of  the  absorbances  at  412  nm 
and  280  nm  (P.N.)  was  never  less  than  0.8  for  any  enzyme 
sample  used  in  the  kinetic  experiments.  Solutions  of  the 

enzyme  were  stored  in  the  cold  at  concentrations  of  about 

-4 

10  M,  and  were  diluted  immediately  before  use.  Enzyme 

concentrations  were  determined  spectrophotometically  at 

5  -1  -1 

412  nm,  using  a  molar  absorptivity  of  1.14  x  10  M  cm 
(Morrison  et_  al.  ,  1957)  . 

Practical  grade  p-cresol  was  obtained  from  Eastman 
Organic  Chemicals,  and  was  purified  by  sublimation . ^ 
Solutions  of  p-cresol  were  prepared  by  weighing  and  the 
concentrations  checked  spectrophotometrically  at  277  nm. 

All  other  substances  were  of  reagent  grade  and  were  used 
without  further  purification.  Hydrogen  peroxide,  30%  by 
weight  obtained  from  the  Fisher  Scientific  Co.,  was  stored 


1 


We  are  grateful  to  Dr.  M. 


K. 


Evett  for  the  purification 


of  the  p-cresol. 


as  a  5  x 10  M  solution,  and  diluted  before  use;  its 
concentration  was  checked  periodically  by  the  method  of 
Ovenston  and  Rees  (1950) .  Water  was  distilled  from  alkaline 
potassium  permanganate  and  then  redistilled. 

Kinetic  investigations  were  carried  out  at  25°  using 
a  stopped-flow  apparatus  essentially  the  same  as  that 
described  previously  (Hasinoff,  1970)  except  that  a  poly¬ 
propylene  cell  was  used  instead  of  one  constructed  of 
plexiglass.  In  a  typical  experiment,  a  2  x 10  M  solution 
of  LP  in  pure  water  was  converted  to  compound  II  (LP-II) 
immediately  before  use  by  the  addition  of  one  equivalent 
of  hydrogen  peroxide.  The  compound  I  (LP-I)  initially 
formed  was  converted  by  reducing  impurities  in  the  enzyme 
solution  to  LP-II  almost  immediately.  This  solution  of 
LP-II  was  allowed  to  react  in  the  stopped-flow  apparatus 
with  a  solution  of  p-cresol  in  the  appropriate  buffer  of 
ionic  strength  0.1.  The  final  reaction  mixture  had  an  ionic 
strength  of  0.05.  At  high  pH,  the  p-cresol  made  a  contri¬ 
bution  towards  both  the  ionic  strength  and  the  buffering 
capacity,  and  the  former  was  allowed  for  by  suitable 
adjustment  of  the  ionic  strength  of  the  buffer. 

The  reduction  of  LP-II  to  LP  was  monitored  at  412  nm 
(the  Soret  band  maximum  of  native  lactoperoxidase)  by  follow¬ 
ing  the  amplified  photomultiplier  voltage  vs.  time  trace  on 
a  564B  Tektronix  storage  oscilloscope.  The  data  were 
recorded  in  the  form  of  a  four  figure  digital  print-out  at 
30  equally  spaced  intervals  of  time,  by  using  an 
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analog-to-digital  converter.  Pseudo  first  order  conditions 
were  achieved  by  keeping  the  concentration  of  p-cresol  in 
at  least  a  ten-fold  excess  over  that  of  LP-II,  which  was 
below  10  ^  M.  The  resulting  absorbance  changes  were  less 
than  0.04,  and  could  hence  be  regarded  as  proportional  to 
the  observed  voltage  changes.  An  average  rate  constant 
with  standard  deviation  was  obtained  from  5-10  traces  for 
each  set  of  conditions.  After  reaction  the  solutions  were 
collected  for  pH  measurement,  which  was  carried  out  with 
an  Orion  digital  pH  meter  in  conjunction  with  a  Fisher 
combination  electrode. 

The  degree  of  dissociation  of  p-cresol  at  the  same 
salt  concentrations  as  were  used  in  the  kinetic  experiments 
was  measured  spectrophotometrically ,  using  a  Cary  14 
spectrophotometer  with  the  cell  compartments  maintained  at 
25° .  The  pH  was  maintained  with  carbonate-bicarbonate  and 
glycinate  buffers,  which  were  also  added  to  the  reference 
solution.  Absorbance  measurements  were  made  at  295  nm , 
close  to  the  maximum  in  the  difference  spectrum  between 
the  acidic  and  basic  forms  of  p-cresol.  Fifteen  solutions 
were  studied  over  the  pH  range  9.39-11.32,  and  absorbance 
measurements  were  also  made  for  solutions  0.01  M  in  perchloric 
acid  and  sodium  hydroxide. 

Results 

The  aqueous  solutions  of  p-cresol  exhibited  maximum 

3  -1 

molar  absorptivities  of  1.70  x  10  M  cm 


at  277  nm  in 


/- 
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0.01  M  perchloric  acid  and  2.54  x  103  M“1cm"1  at  295  nm 

in  0.01  M  sodium  hydroxide.  These  values  compare  favourably 

with  the  maximum  values  of  1.71  x  103  M_1cm"1  at  277  nm  and 
3  -1  -1 

2.55  x  10  M  cm  at  295  nm,  respectively,  found  for  a 

sample  of  99.96%  purity  (Herington  and  Kynaston,  1957). 

The  effective  pK  of  p-cresol  under  the  conditions  of 

a 

the  present  work  may  be  defined  in  terms  of  concentrations 
and  the  operational  pH  scale  by  the  relation 

pK  =  log— — +  pH  (1) 

[ARCf] 

where  AROH  and  ARO  represent  the  protonated  and  unprotonated 
forms  of  p-cresol,  respectively.  Non-linear  least  squares 
analysis  of  the  spectrophotometric  data,  using  the  measured 
value  of  80  M  ^cm  3  for  the  molar  absorptivity  of  p-cresol 
at  295  nm  in  acid  solution,  yielded  a  pK^  value  of 
10.11  ±0.01.  The  largest  deviation  from  the  best  fit  curve 

was  of  4%  absorbance,  and  the  predicted  molar  absorptivity 

3-1-1  3 

of  the  anion  was  2.71  x 10  M  cm  compared  with  2.54  x 10 

M  3cm  ^  obtained  by  direct  experimental  measurement.  This 

value  of  the  pK^  is  in  agreement  with  the  value  10.10 

obtained  by  the  use  of  Davies'  equation  (1938)  together  with 

the  value  at  zero  ionic  strength  found  by  Chen  and  Laidler 

(1962)  . 

All  reaction  traces  were  simple  exponential  curves, 
showing  the  reaction  to  be  kinetically  first  order  in 
concentration  of  LP-II.  The  results  were  analyzed  by  a 


? 
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non-linear  least  squares  method  on  a  computer  with  equal 
weighting  of  all  voltage  readings. 

Over  the  entire  pH  region  of  this  study,  the  observed 
first  order  rate  constant  depended  on  the  concentration  of 
p-cresol  in  a  simple  linear  fashion : 


k  , 
obs 


k  [PC] 


(2) 


where  k^g  is  the  experimental  first  order  rate  constant, 
k  is  a  second  order  rate  constant,  and  [PC]  is  the  total 
concentration  of  p-cresol.  This  is  shown  in  Fig.  5-1,  a 
plot  to  check  the  order  of  the  reaction  with  respect  to 
p-cresol  at  pH  7.87,  from  which  k  is  (1.17  ±0.04)  x 10^ 

M  ^sec  ^  and  the  intercept  is  0.74  ±1.13  sec  ^ .  The  line¬ 
arity  of  the  plot  confirms  that  the  reaction  is  first  order 
with  respect  to  the  concentration  of  p-cresol.  The  rate  of 
spontaneous  decay  of  LP-II  was  negligible  compared  to  the 
fast  rates  observed  in  the  presence  of  p-cresol.  At  pH 
values  below  10,  the  range  of  concentration  of  p-cresol 

that  could  be  employed  was  restricted  to  one  order  of 

.  .  -5 

magnitude;  the  lower  limit  was  10  M  due  to  the  restriction 

of  observing  pseudo  first  order  conditions,  and  the  upper 

-4 

limit  was  10  M  because  the  rates  of  reaction  became  too 
fast  to  study  using  the  stopped-flow  apparatus.  Table  5-1 
contains  values  of  k  (with  standard  deviation)  as  a  function 
of  pH  over  the  entire  pH  range. 

The  plot  of  log  k  vs.  pH  in  Fig.  5-2  is  well  defined. 


and  can  be  seen  with  the  aid  of  Dixon's  rules  (Dixon  and 
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Table  5  1.  Values  of  k  (with  standard  deviation) 
function  of  pH. 


as  a 


V  vr-l  -1 

k,M  sec 

PH 

Buffer^" 

(9. 48  ±  0 . 78)  x  105 

2.08* 

GH 

(1. 67  ±  0.22)  x  106 

2.53 

GH 

(1.63  ±  0.21)  x  106 

2.71 

GH 

(1.94  ±  0.25)  x  106 

2.93 

F 

(2.32  ±  0.64)  x 106 

2 . 96* 

F 

(1. 83  ±  0.28)  x 106 

3.15 

F 

(1.85  ±  0.32)  x 106 

3 .35 

F 

(2 . 30  ±  0 .28)  x  106 

3.78 

F 

(2.17  ±  0.22)  x  106 

3.99 

F 

(2 . 66  ±  0 .14)  x  106 

4 . 00* 

F 

(2.27  ±  0.14)  x  106 

4.22 

F 

(2.40  ±  0.23)  x 106 

4.52 

F 

(2 . 13  ±  0. 40)  x  106 

4 . 75 

F 

g 

(2 . 09  ±  0.13)  x 10 

4 .96* 

A 

(2.23  ±0.24)  x 106 

4.96 

A 

(2 . 01  ±  0.15)  x 106 

5.19 

A 

(2.05  ±  0.36)  x 106 

5.20 

F 

(2.06  ±  0.09)  x 106 

5.37 

A 

(1.93  ±  0.13)  x 106 

5.53 

A 

(1.66  ±  0.10)  x 106 

5.83 

P 

(1.52  ±0.19)  x  106 

6 .00 

P 

(1. 61  ±  0. 03)  x 106 

6.01* 

P 

(1.56  ±  0.10)  x 106 

6.21 

P 

(1.31  ±  0.28)  x 106 

6.38 

P 

(1.47  ±0.16)  x  106 

6.62 

P 

(1.38  ±0.11)  x  1 0 6 

6.79 

P 

v  6 

* 

(1.28  ±  0.03)  xlO 

6.97 

P 

(1.30  ±0.08)  xlO6 

6.99 

P 

(1.27  ±0.09)  xlO6 

7.20 

P 

(1.09  ±0.17)  xlO6 

7.34 

T 

continued 


Table  5-1  continued 
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(1.09 

+ 

0.15) 

x  106 

7.61 

T 

(9 . 00 

+ 

1.89) 

x  105 

7.81 

T 

(1.17 

+ 

0.04) 

x  106 

7.87* 

T 

(1.02 

+ 

0.06) 

x  106 

8 . 03 

T 

(1.12 

+ 

0.08) 

x  106 

8.16 

P 

(1.10 

+ 

0.12) 

x  106 

8.27 

T 

(8 . 30 

+ 

1.57) 

x  105 

8.48 

T 

(9 . 50 

+ 

0 . 89) 

x  105 

8 .68 

T 

(8 .81 

+ 

0.79) 

x  10^ 

9.06 

C 

(7.82 

+ 

0.40) 

x  105 

9.24 

C 

(6.57 

+ 

1 . 31) 

x  105 

9.43 

C 

(6 . 35 

+ 

0.68) 

x  105 

9.46* 

c 

(5.40 

+ 

0.25) 

x  105 

9.59 

c 

(4.65 

+ 

0.25) 

x  105 

9.69 

G 

(3.37 

+ 

0.06) 

x  105 

9.87 

G 

(3.17 

+ 

0.16) 

x  105 

9 . 90* 

G 

(2.26 

+ 

0.19) 

x  105 

10.05 

G 

(1.71 

+ 

0.17) 

x  105 

10.16 

G 

(1.53 

+ 

0.23) 

x  105 

10.20 

G 

(1.05 

+ 

0 . 05) 

x  105 

10.33 

G 

(5.28 

+ 

0.46) 

x  104 

10.54 

C 

(5.09 

+ 

0.53) 

x  104 

10.54 

G 

(5.10 

+ 

0.37) 

4 

x  10 

10.55 

G 

(1.60 

+ 

0 . 15) 

x  104 

10.86 

G 

(4.85 

+ 

0.42) 

x  103 

11.15* 

C 

(4.26 

+ 

0.58) 

x  103 

11.18 

G 

^Buffer  key:  GH ,  glycine  hydrochloride-glycine;  F,  formic 
acid-sodium  formate;  A ,  acetic  acid-sodium  acetate;  P, 
potassium  dihydrogen  phosphate-disodium  hydrogen  phosphate; 

T ,  tris (hydroxymethyl ) aminomethane  hydrochloride-tris (hydroxy¬ 
methyl  )  aminomethane  ;  C,  sodium  bicarbonate-sodium  carbonate; 

G,  glycine-sodium  glycinate. 

pH  values  marked  by  an  asterisk  are  those  at  which  order 
studies  in  [PC]  were  performed.  All  other  rate  constants  are 
results  of  determinations  at  a  single  concentration  of  p-cresol. 


' 
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Fig.  52.  A  plot  of  log  k  vs.  pH  for  the  reaction  between  LP-II  and  p-cresol.  The 
solid  line  is  a  curve  calculated  by  computer  for  the  fitting  of  the  data  to  Eq.  3.  The 
error  bars  represent  standard  deviations  in  the  rate  constant,  k 
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Webb,  1964)  to  depend  on  the  acid  dissociation  of  the 

substrate  at  pH  10  and  the  ionization  of  three  kinetically 

important  enzyme  groups  having  pK  values  of  around  2,  6, 

a 

and  9.5.  The  minimum  reaction  scheme  consistent  with  the 
data  is  given  by  Scheme  I,  in  which  P-II ,  HP-II, 
and  H^P-II  are  the  four  kinetically  distinct  states  of 
protonation  of  compound  II,  the  K's  are  acid  dissociation 
constants,  and  the  k's  are  bimolecular  rate  constants: 


Scheme  I 


K. 


h3p-ii 

A 


k 


h2p-ii 


K. 


K 


V 

HP-II 
A 


/ 

P-II 


+ 


+ 


+ 


AROH 
A 


K 

^  a 


ARO 


AROH 
A 


K 

f  s 


ARO 


AROH 
A 


K 

^  a 


ARO 


AROH 
A 


K 

v  a 


ARO' 


-> 


'8 


Products 


The  vertical  arrows  denote  fast  non-rate  determining 
proton  transfers.  The  corresponding  kinetic  expression  is 
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k 


k1[H+]4 


KaK2K3K4 


+ 


[H+l3 

KaK3K4 


+ 


+ 


2 


4 


+ 


+ 


Va\ 

K4  / 


K 

a 


+ 


[H+l  [ H+ 1  2 
K4  K3K4 


isV-V] 

K2K3K4/V 


+ 


. (3) 

In  a  reaction  involving  two  or  more  ionizations  such  as 
that  depicted  by  Scheme  I ,  molecular  rather  than  group 
ionizations  are  necessarily  implied  (Dixon  and  Webb,  1964). 

Group  ionization  constants,  of  which  the  molecular  ioniza¬ 
tion  constants  consist,  cannot  be  determined  from  the  analysis 
of  rate  data  alone.  If  the  pK  values  of  the  group  ionization 
constants  are  widely  separated,  the  ionization  constants 
measured  indicate  the  true  value  of  the  group  ionization 
constants.  The  example  of  the  unsymmetrical  dibasic  acid 
given  by  Dixon  and  Webb  (1964)  illustrates  clearly  the 
distinction  between  group  and  molecular  ionization  constants : 


R-CHC0o 

K1  ^ 
lx^ 

CH2C02H 

R-CHC02H 

- - 

R-CHCCb 

1  ^ 

CH2C02H 

CH2C02 

R-CHCCbH 

•^K 

2x 

1  2 

CH„CO 
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Each  group  may  ionize  either  from  the  undissociated 
molecule  or  from  the  anion,  in  other  words  in  the  first  or 
second  stage;  we  therefore  need  four  ionization  constants 
to  represent  the  ionization  of  an  unsymmetrical  dibasic 
acid  completely .  The  relationship  between  the  molecular 
ionization  constants  and  ,  and  the  group  ionization 
constants  K^x ,  K^,  K^,  K ^  is  given  by  the  following 
expressions : 


K,  =  K,  +  K., 

1  lx  ly 


K. 


2x  2y 

K0  + 

2x  2y 


Two  pathways  of  Scheme  I  such  as 


H3P-II  +  ARO 


-> 


and 


H2P-II  +  AROH  — — v. 


(4) 

(5) 


(6) 

(7) 


are  kinetically  indistinguishable  and  are  combined  together 
in  one  term  of  the  expression  relating  k  to  the  specific 
rate  constants,  ionization  constants,  and  hydrogen  ion 
concentrations.  The  only  unambiguous  rate  constants  obtain¬ 
able  from  Eq.  3  are  k^  and  kg.  The  differential  rate  lav/ 
for  the  reaction  represented  by  k^  is 

d[H  P-II] 

-  - - -  =  k9  [H^P-II]  [ARO  ] 

dt 


and  for  the  reaction  represented  by  k^  is 


(8) 
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d [H2P-II] 

- -  =  k  [H0P-II] [AROH] 

dt  J  2 

(9) 

but  [H  P-II] 

[H  P-II] 

=  K  - - 

[H  ] 

(10) 

and  [AROH]  = 

[H+l  [ARO'1 

(ID 

K 

a 

Equations  (10)  and  (11)  substituted  in  equation  (9)  give 

d[H  p-li]  k  K 

- -  - -[H-P-II]  [ARO  ]  (12) 

dt  K  J 

a 

which  when  compared  with  Eq.  (8)  is  seen  to  be  of  identical 
form.  Hence  the  two  reactions  (6)  and  (7)  are  kinetically 
indistinguishable . 

The  k  vs.  pH  data  were  analyzed  using  a  non-linear 
least  squares  program  (described  in  Appendix  2) .  The  smooth 
curve  in  Fig.  5-2  shows  the  best  fit  to  Eq .  3,  and  is 
defined  by  the  first  column  of  parameters  in  Table  5-2. 

An  alternative  and  algebraically  more  simple  method  of 

analyzing  the  log  k  vs.  pH  plot  which  avoids  the  need  to 

consider  kinetically  indistinguishable  or  insignificant  rate 

processes  is  provided  by  the  use  of  transition  state  acid 

dissociation  constants.  As  is  shown  elsewhere  (Critchlow  and 

Dunford,  submitted  for  publication) ,  an  acid  dissociation 

in  the  transition  state  gives  rise  to  an  increase  in  slope 

in  the  plot  of  log  k  vs.  pH,  whereas  an  acid  dissociation  in 

the  ground  state  causes  a  decrease  in  slope.  The  number 

and  approximate  pK  values  of  such  ionizations  may  therefore 

a 


» 
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Table  5  2.  Parameter  values3  and  standard  deviations  for 
the  variation  of  k  with  pH,  obtained  by  non-linear  least 
squares  fit  of  data  to  Eqs.  (3)  and  (15). 


Equation  3 

Equation  15 

Equation  15 

pK  variable 
a 

pK  variable 

3 

pK  invariable 

3 

kl 

ob 

(k3 

+ 

k  K 

2  a. 

K2 

(2 . 30+0 . 05) xio6 

(k5 

+ 

k4Ka, 

K3 

(1.09+0. 03) xlO6 

(k7 

+ 

k^K 

6  a^ 

K4 

b 

0 

00 

ob 

pk2 

2. 3+0.1 

2.310.1 

2.210.1 

pk3 

5.810.1 

5.810.1 

5.910.1 

pK4 

9.710.1 

9.710.1 

9.810.1 

PKa 

10.3+0.1 

10.3+0.1 

10.11 

PK{ 

6.110.1 

6.210.1 

V 

k 

(2 .39+0 .05) xlO6 

(2 . 37  +0 . 05) xlO6 

aRate  constants  in  units  of  M  ^sec  ^ 


'zero  within  experimental  error. 
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be  read  off  the  plot  by  an  extension  of  Dixon's  rules 
(Dixon  and  Webb,  1964)  in  the  same  way  as  the  pK  values 
of  the  enzyme  or  enzyme-substrate  complex.  The  essential 
significance  of  this  treatment  may  be  arrived  at  more  simply 
by  considering  the  experimental  log  k — pH  profile  in  a  less 
complicated  plot,  as  in 


The  change  of  slope  on  passing  from  the  horizontal  portion 
of  the  plot  at  high  pH  to  the  portion  of  unit  negative  slope 
must  be  caused  by  a  change  in  the  reaction  route  from  one 
involving  a  given  set  of  reactants  to  one  involving  a 
different  set  containing  a  total  of  one  more  proton.  The 
transition  state  at  medium  (and  low)  pH  therefore  contains 
one  more  proton  than  that  for  reaction  at  high  pH,  and  the 
point  at  which  the  curve  begins  to  rise  in  going  to  lower 
pH  may  be  associated  formally  with  a  transition  state  acid 
dissociation.  The  subsequent  return  to  a  horizontal  relation¬ 
ship  at  lower  pH  is  of  course  due  to  the  protonation  of  the 
ground  state  and  is  covered  by  Dixon's  rules  (Dixon  and  Webb, 
1964)  . 
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The  use  of  this  treatment  in  a  description  of  the 

mechanism  of  a  reaction  requires  that  each  transition- 

state  acid  dissociation  constant  (pK^  )  be  paired  with  that 

for  the  ground-state  dissociation  (pK  )  considered  to  involve 

a 

the  same  ionizing  group,  although  both  members  of  each  pair 
may  not  be  found  in  the  pH  range  of  the  study.  The  pK 
difference  between  the  ground-  and  transition-states,  which 
may  be  regarded  as  expressing  the  sensitivity  of  the  rate 
to  the  ionization  in  question,  provides  a  qualitative 
measure  of  the  change  in  environment  accompanying  the 
activation  process.  For  example,  if  protonation  of  a 
particular  protein  residue  facilitates  a  reaction,  then 
PP^  >  P-^a '  as  "'"S  case  with  the  simple  reaction  considered 

above.  Conversely,  if  pK^  < pK  ,  then  protonation  of  the 
residue  will  retard  the  rate  of  the  reaction. 

In  the  actual  derivation  of  an  equation  describing 
the  observed  rate  constant  as  a  function  of  pH,  use  is  made 
of  dimensionless  Michaelis  functions  (Dixon  and  Webb,  1964)  , 
and  the  final  result  is  (Critchlow  and  Dunford,  submitted 
for  publication)  : 


k  , 
obs 


(13) 


where  k  is  a  pH-independent  rate  constant,  the  f's  are 
Michaelis  functions,  E  is  the  enzyme,  A,  B,  ...  z  are 
reactants,  and  the  subscripts  refer  to  particular  states 


of  ionization  of  the  enzyme,  substrates,  and  transition 
state.  For  example,  the  corresponding  equation  for  the 
simple  case  shown  above  is 


‘obs 


(i  +  KbrH+i) 

(1  +  K  /  [H+]  ) 


(14) 


For  the  oxidation  of  p-cresol  by  LP-II,  consideration 
of  Fig.  5-2  shows  a  single  kinetically  important  transition 
state  ionization  of  pK^  value  about  6,  so  that  the  reaction 

Cl 

may  be  represented  by  Scheme  II, 


Scheme  II 
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where  the  species  in  square  brackets  represent  transition 
state  complexes  in  different  states  of  protonation,  is 

the  associated  transition  state  acid  dissociation  constant, 
and  the  other  symbols  have  been  defined  for  Scheme  I.  The 
resulting  kinetic  expression  is 
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k  =  k 


1  +  K^/  [H+] 


/IH^I 


K. 


+  1  + 


+ 


[H+]  [H+]2/ 


(15) 


K 


[H  V 


where  k  is  the  pH -independent  rate  constant  associated  with 
the  pH-independent  part  of  the  curve  between  pH  3  and  pH  5 . 
The  significance  and  advantages  of  this  style  of  treatment 
of  pH-rate  effects  in  terms  of  transition  state  acid 
dissociation  constants  are  discussed  in  more  detail  else¬ 
where  (Critchlow  and  Dunford,  submitted  for  publication), 
kut  the  readily  apparent  advantage  of  simplicity  can  be  seen 
by  comparing  Eqs.  3  and  15. 

Non-linear  least  squares  treatment  of  the  data  using 
Eq.  15  resulted  in  the  same  smooth  curve  in  Fig.  5-2.  The 
best  fit  values  of  the  parameters,  with  pK  values  identical 
to  those  obtained  by  the  former  treatment,  are  given  by  the 
second  column  of  parameter  values  in  Table  5-2.  The  value 
of  pK^  for  p-cresol  obtained  kinetically  is  different  from 
the  spectrophotometric  result  mentioned  previously;  however, 
a  second  analysis  with  a  fixed  value  of  pK^  =  10.11  produced 
only  small  changes  in  the  other  parameters,  as  shown  in  the 
third  column  of  Table  5-2.  A  good  fit  to  the  experimental 
data  was  still  obtained,  so  that  the  difference  in  pK 

a 

values  may  not  be  significant. 


Discussion 


A  number  of  conclusions  may  be  drawn  from  the  nature 
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of  the  pH  rate  profile  for  the  oxidation  of  p-cresol  by 
compound  II  of  lactoperoxidase .  Critchlow  and  Dunford 
(submitted  for  publication)  have  performed  an  analogous 
study  of  the  oxidation  of  p-cresol  by  compound  II  of  horse¬ 
radish  peroxidase  (HRP-II) ,  and  they  obtained  a  pH-rate 
profile  very  similar  to  that  obtained  in  this  study.  They 
also  obtained  evidence  for  an  acid  group  on  the  enzyme  with 
a  ^Ka  va^ue  °f  2.3.  As  they  have  suggested,  the  decrease 
in  rate  associated  with  this  pK  may  arise  from  an  initial 
step  in  the  acid  catalyzed  splitting  of  the  porphyrin  ring 
from  the  apoprotein,  perhaps  by  protonation  of  the  acceptor 
atom  of  a  hydrogen  bond  involved  in  the  heme-protein 
stabilizing  linkage.  The  pK  of  LP-II  at  pH  5.8  is  notable 
for  its  closeness  to  the  transition  state  pK^  at  pH  6.1, 
which  may  mean  protonation  of  the  same  group;  if  such  is  the 
case,  protonation  of  this  group  would  then  exert  only  a 
slight  accelerating  influence  on  the  rate  of  reaction,  indicat¬ 
ing  that  the  group  is  positioned  close  to  the  active  center 
of  reaction,  but  plays  little  part  in  the  catalytic  process. 

It  is  also  possible  to  draw  some  conclusions  about  the 

nature  of  the  reactive  form  of  p-cresol.  Inspection  of  the 

pH-rate  profile  shows  that  if  p-cresol  were  reacting  in  the 

anionic  form,  the  second  order  rate  constant  would  exceed 

the  diffusion  controlled  limit  ( 1 0 0  M  ^sec  ^  or  less  for 

an  enzyme-substrate  system  in  which,  also,  the  direction  of 

approach  of  the  substrate  molecule  is  probably  important) 

4= 

below  pH  6.5.  Thus  the  pK^  at  pH  6.1  cannot  be  paired  to 
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the  pK^  of  p-cresol  at  pH  10.3  and  therefore  the  transition 

state  pK^  pertaining  to  the  protonation  of  p-cresol  must 

occur  outside  the  experimentally  covered  pH  range  to  high  pH. 

This  indicates  that  the  acid  form  of  the  substrate  is  the 

reactive  species  over  the  entire  pH  range  covered.  The  same 

conclusion  is  reached  when  Scheme  I  is  used.  In  order  to 

decide  which  of  the  kinetically  equivalent  pathways  included 

in  each  of  the  two  non-zero  rate  terms  of  Eq.  3  in  Table  5-2 

is  responsible  for  the  reaction  it  is  necessary  to  make  use 

of  the  maximum  value  which  can  be  placed  on  the  rate  of  a 

diffusion  controlled  process.  Since  K  /Kn  =  10~8  and 

a  2 

-  5 

^a/K^  =  3  x 10  ,  it  follows  that  even  for  the  maximum 

diffusion  controlled  limit  ( 10 1 0  M~1sec~1)  the  predominant 
pathways  are  those  described  by  k3  and  k5 .  Thus  the  proton- 
ated  forms  of  compound  II  would  appear  to  react  with  the 
neutral  form  of  the  substrate. 

It  is  possible  that  the  transition  state  pK^  at  pH  6.1 
could  be  paired  to  any  of  the  three  enzyme  ground  state 
pK ' s ;  however,  if  the  5.8  -  6.1  pairing  scheme  is  accepted, 
the  pK  at  9.7  can  only  be  paired  (again  by  using  the  concept 
of  a  second  order  rate  constant  exceeding  the  diffusion 
controlled  limit)  with  a  transition  state  pK^  lying  outside 
the  observed  pH  range  at  high  pH.  This  is  equivalent  to 
saying  that  the  group  of  pK  9.7  must  be  present  in  its  acid 
form  for  reaction  to  occur.  The  value  of  pK^  =  9.7  is 
notable  in  that  it  is  1  pH  unit  higher  than  the  value  of 
8.6  obtained  from  the  studies  of  the  reaction  of  HRP-II  with 
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ferrocyanide  (Hasinoff  and  Dunford,  1970)  and  p-cresol 
(Critchlow  and  Dunford,  submitted  for  publication) ,  which 
may  indicate  that  the  identity  of  the  ionizing  group  is 
different,  or  that  its  ionization  constant  is  shifted 
through  some  sort  of  interaction  such  as  hydrogen  bonding 
in  the  acid  form.  In  spite  of  this,  the  similarity  of  the 
pH-rate  profiles  for  the  oxidation  of  p-cresol  by  LP-II  and 
HRP-II  gives  some  indication  that  the  mechanisms  of  reaction 
in  both  cases  are  the  same,  except  for  the  unproductive 
binding  observed  in  the  latter  case.  This  work  can  also  be 
compared  to  the  reaction  between  LP-II  and  iodide  ion,  in 
which  the  second  order  rate  constant  varies  with  unit 
negative  slope  over  the  experimental  pH  range,  from 
4 . 2  x 10  M  sec  at  pH  2 . 9  to  2.6x10  ^  M  ^sec  ^  at  pH  10 . 1 
(Maguire  and  Dunford,  1972) .  Thus  the  second  order  rate 
constant  for  p-cresol  oxidation  by  LP-II  is  appreciably 
greater  than  the  rate  constant  for  iodide  oxidation  over 
most  of  the  pH  range,  except  for  the  low  pH  region.  The 
radical  difference  in  the  pH-rate  profiles  for  these  two 
substrates  (three  enzyme  pK  's  observed  in  the  LP-II — p-cresol 
reaction,  none  in  the  pH  range  2.9  -  10.1  observed  for  the 
LP-II — iodide  reaction)  may  be  a  result  of  the  lack  of 
steric  interactions  between  iodide  ion  and  LP-II,  or  it 
may  be  the  result  of  the  operation  of  a  different  mechanism. 
Since  p-cresol  reacts  in  its  unionized  form,  electrostatic 
interactions  cannot  be  a  factor  in  its  reaction. 
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Chapter  6 

The  Rate  of  Exchange  of  Water  with  Iron  of 

Horseradish  Peroxidase  Detected  by  017 
Nuclear  Magnetic  Resonance 

Introduction 

Horseradish  peroxidase  is  an  enzyme  which  catalyzes 
oxidations  in  the  presence  of  hydrogen  peroxide  and  has 
been  given  the  following  designation  by  the  Enzyme  Commission 
of  the  International  Union  of  Biochemistry  (1961):  1.11.1.7 

Donor:  ^2^2  ox^ore^uct3se.  The  first  isolation  of  horse¬ 

radish  peroxidase  (HRP) ,  in  partially  purified  form,  was 
carried  out  by  Bach  and  Chodat  (1903) .  It  was  crystallized 
by  Theorell  (1942)  using  electrophoretic  techniques.  HRP 
is  now  available  commercially  in  various  degrees  of  purity 
from  a  number  of  firms.  HRP  has  a  molecular  weight  of 
40,200  (Maehly,  1955) ,  and  it  has  been  found  that  about  18% 
by  weight  of  HRP  is  carbohydrate  material  (Shannon  et  al . , 

1966;  Theorell  and  Akeson,  1942),  although  the  nature  of  the 
linkage  to  the  protein  is  not  known.  HRP,  like  hemoglobin 
myoglobin,  and  catalase,  contains  as  a  prosthetic  group 
protoporphyrin  IX,  shown  in  Figure  1-la  of  Chapter  1.  The 
hematin  nature  of  HRP  was  established  by  Keilin  and  Mann  (1937) , 
and  the  importance  of  the  heme  group  was  demonstrated  by 
Theorell  (1940)  who  found  that  neither  the  protein  nor  the 
prosthetic  group  independently  has  any  appreciable  peroxidatic 
activity.  In  its  natural  form,  the  iron  of  the  heme  group  of 
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HRP  is  in  the  ferric  state.  To  date  there  has  been  no 
X-ray  crystallographic  work  reported  for  HRP ,  so  evidence 
regarding  the  nature  of  the  heme-protein  forces  is  mostly 
based  on  a  comparison  of  the  properties  of  HRP  with  those 
of  hemoglobin  and  myoglobin.  There  has  been  a  considerable 
amount  of  speculation  in  regard  to  the  identity  of  the 
ligands  occupying  the  fifth  and  sixth  coordination  positions 
on  the  heme  iron  of  HRP.  It  is  generally  agreed  that  a 
water  molecule  occupies  the  sixth  coordination  position 
(Brill  and  Williams,  1961;  Brill,  1966  ;  Nicholls,  1962). 

Nuclear  magnetic  resonance  (NMR)  methods  have  been 
extensively  applied  to  the  investigation  of  biological 
systems  (Kowalsky  and  Cohn,  1964) .  High  frequency  NMR 
(Kowalsky,  1962),  relaxation  time  measurements  (Jardetsky, 
1964),  and  halide  ion  probe  techniques  (Stengle  and 
Baldeschweiler ,  1966;  Marshall,  1968)  have  been  used  to 
determine  the  structure,  conformation,  and  motion  of 
enzymes  and  proteins  in  solution.  While  interpretation  of 
the  NMR  spectrum  of  a  macromolecule  is  limited  by  its 
inherent  complexity,  much  information  can  be  obtained  from 
the  NMR  spectrum  of  a  small  molecule  that  is  able  to  probe 
the  environment  of  the  macromolecule  by  rapid  exchange 
between  free  solution  and  attachment  to  the  macromolecule. 
This  study  was  undertaken  in  an  attempt  to  determine  the 
exchange  rate  and  activation  parameters  involved  in  the 
interaction  of  the  ferric  iron  of  the  heme  group  of  HRP  with 


water . 


Let  us 

consider 

the  Bloch  equations  in 

the  rotating 

frame  of  coordinates 

(Pople  et  al.,  1959)  : 

du 

u 

dt 
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If  a  complex  moment  of  magnetization,  G,  is  defined  as 


G  =  u  +  iv 


(3) 


then 


dG 

dt 


du  .  dv 

dt  1  dt 


(4) 


or 


dG 

dt 


+ 


T. 


i  (oJq-o) 


G  =  -iyH^Q 


(5) 


If  exchange  of  a  nucleus  between  two  sites  A  and  B  is 
allowed  for,  then 
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where  and  tb  are  the  lifetimes  of  a  nucleus  in  each 


environment 
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If  a  steady  state  is  assumed  in  which 


dG 


A 


dG 


B 


dt  dt  0  ' 

and  if  it  is  assumed  that  A  is  the  dominant  site  and  that 
the  contribution  is  negligible  (this  could  correspond 

to  water  exchange  between  bulk  solvent  (site  A)  and  a 
paramagnetic  ion  (site  B) ) ,  then 
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This  derivation  follows  that  of  Swift  and  Connick  (1962) 
for  two  sites. 

When  this  equation  is  compared  to  that  representing  no 
exchange ,  and  when  respective  real  and  imaginary  parts  are 
equated,  the  following  relations  result: 
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20BS 


=  X 
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tatb(Xb  +Aq3B  ) 
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TATB(AB  +Aa)B  ^ 


(9) 


(10) 
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An  experimental  situation  is  depicted  in  Figure  6-1  in  which, 
saY /  the  addition  of  a  paramagnetic  ion  to  water  results  in 
a  frequency  shift  of  the  solvent  peak.  The  half-width  at 


half-height  of  the  broadened  peak  is  equal  to  1/T  ,  and  1/T 

Z  Z 

is  equal  to  the  half— width  at  half-height  of  the  original 
peak. 


If  we  let 

1  =  1  _  1 

T  T  T 

2p  20BS  2A 


then 


(11) 


2p 


A 


(T2B)  ~2  +  <T2BTBrl  +  A(0B2 

((T2B)1  +  (TBr1)2  +  A“B2 


(12) 


The  temperature  effects  predicted  by  equation  (12)  were 
investigated  by  Swift  and  Connick  (1962)  ,  who  considered 
the  following  limiting  cases: 

Case  A 


>> 


1  1 
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then 
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(14) 


Relaxation  occurs  through  a  change  in  the  precessional 


frequency  and  is  rapid; 


is  controlled  by  the  rate  of 
^  n  [M] 


chemical  exchange.  is  given  by  ^ 


where  n  is  the 
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number  of  coordinated  water  molecules  per  ion  and  [M]  is 
the  molar  concentration  of  the  paramagnetic  ion. 

Case  B 


then 


>> 


>> 


T  T 
2B  B 


(15) 


Chemical  exchange  is 
of  relaxation  through 
Case  C 


tbAwb' 


(16) 


rapid;  — —  is  controlled  by  the  rate 
2p 

the  precessional  frequency. 


>> 


(17) 


then 
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Relaxation  by  T2b  is  fast;  - —  is  controlled  by  the  rate 

2p 

of  chemical  exchange. 

Case  D 
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Chemical  exchange  is  rapid; 
relaxation  process. 
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The  temperature  dependence  of  t  is  given  by 

B 


T 


B 


(21) 


T 


B 


(22) 


where  AH^  and  AS^  are  the  enthalpy  and  entropy  of  activation 
for  the  first  order  reaction  of  exchange  of  water  from  a 
cation,  and  k^  is  the  rate  constant  associated  with  such  a 
process. 

The  temperature  effects  predicted  by  the  preceding  four 

cases  are  illustrated  schematically  in  Figure  6-2  as  log 

— ^  as  a  function  of  .  From  consideration  of  experimental 

data  plotted  in  this  way,  one  may  be  able  to  determine  which 

particular  "kinetic  region"  describes  the  experimental 

system,  and  thus  it  may  be  possible  to  calculate  t  ,  hence 

B 

the  rate  constant  for  exchange  of  water  molecules  with  the 
ferric  iron  of  horseradish  peroxidase. 


Experimental  Methods 

17  .  .  17 

H^O  containing  4.56  atom  %  0  was  obtained  from 

Miles  Laboratories.  HRP  was  obtained  from  the  Boehringer- 

Mannheim  Corporation  with  a  P.17,  of  0.6  (for  HRP,  P.N.  is  the 

ratio  of  the  absorbance  at  403  nm  to  that  at  280  nm) .  The 

enzyme  was  purified  to  a  P.N.  of  2.8  by  gel  filtration  (on 

Sephadex  G-200  Superfine,  using  pH  7.0  phosphate  buffer  of 

ionic  strength  0.05  at  4°) ,  and  then  dialysis  against  doubly 

distilled  water,  followed  by  lyophilization .  For  HRP, 

samples  exhibiting  P.N.  of  3.0  are  commonly  considered  pure. 


Ill 
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Solutions  for  NMR  analysis  were  prepared  by  dissolving  the 

lyophilized  material  in  the  H  O'1'7. 

2 

HR?  solutions  (with  an  external  reference  of  I^O'*'7)  at 
six  pH  values  from  5.6  to  10.0  were  examined  by  NMR  at  temp¬ 
eratures  from  4°  to  55°.  The  sample  volume  was  about  2.5  ml. 
A  typical  HRP  concentration  was  2.20  x  io"4  m,  as  determined 
spectrophotometrically  (Cary  14)  with  an  extinction  coeffi- 
ient  of  9.1  x 10  M  1cm  1  at  403  nm  (Keilin  and  Hartree,  1951). 

The  NMR  spectra  were  recorded  on  a  Varian  model  HA-601 
spectrometer  equipped  with  a  Hewlett-Packard  5100  B 


Frequency  Synthesizer  and  a  standard  variable  temperature 
probe  and  temperature  controller.  The  temperature  of  the 
sample  was  controlled  by  a  stream  of  nitrogen,  previously 
heated  or  cooled,  and  calibrated  with  a  thermocouple  using 
ice-water  as  a  reference.  Operating  conditions  were  14,000 
gauss  and  8.116  MHz.  Resonance  was  recorded  as  the  deriva¬ 
tive  of  the  absorption  mode  (the  separation  of  the  two  peaks 
of  the  derivative  curve  is  2//3T2),  and  linewidth  calibrations 

were  made  with  the  frequency  synthesizer. 

17 

An  H20  reference  solution  was  used  before  and  after 
each  sample  measurement  at  each  temperature.  There  was  no 
detectable  shift  in  the  resonance  line  due  to  the  paramagnetic 
ion  introduced.  Table  6-1  contains  values  of  1/T  '(1/T  '  = 

7  2p  K  /x2p 

1/T9  xi/ [metal  ion]  at  six  different  values  of  pH  and  eight 
z  p 

to  ten  temperatures.  Figures  6-3  to  6-8  are  linear  least 
squares  plots  of  log  (1/T  ' )  vs.  1/T°K  for  each  pH  value. 

The  estimated  error  limit  per  point  is  +10%. 
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Table  6-1 
(1/T2  )  vs.  T°C 


PH  10.15  pH  8.87  pH  7.92 


T°C 

(1/T 

2p,>Xl0~5 

T°C 

(1/T 

2p,)Xl0"5 

T°C 

(1/T2p')xl0-5 

4 

2.23 

4 

1.71 

4 

2.05 

15 

0.80 

15 

1.35 

15 

1.98 

20 

0.54 

20 

0.95 

20 

1.82 

25 

0.74 

25 

1.05 

25 

2 . 04 

30 

0.81 

30 

1.05 

30 

2.02 

32 

0.85 

32 

0.98 

32 

1.95 

34.5 

0.83 

34.5 

1.23 

34.5 

1.97 

39 

0.66 

39 

0.85 

39 

1.71 

44 . 5 

0.59 

44.5 

0.74 

44.5 

1.79 

50 

0.64 

50 

0.76 

50 

1.53 

pH  6 

.85 

pH  5 

.94 

pH  5 

.61 

T°C 

(1/T 

i  -  g 

O  )  x  1 0 

2p 

T°C 

(1/T 

9  ')xl0~6 

2p 

T°C 

(1/T 

2p,>Xl0~6 

4 

0.88 

4 

0.41 

4 

0.79 

15 

0.50 

15 

0.26 

15 

0.54 

20 

0.42 

20 

0.17 

20 

0.45 

25 

0.33 

25 

0.15 

25 

0.48 
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Results  and  Discussion 

A  schematic  der ivative-of -absorption  curve  is  shown 
in  Fig.  6-9.  As  mentioned  above,  only  one  resonance  line  was 
seen;  there  was  no  observed  shift  in  the  line  due  to  the 
paramagnetic  ion. 

The  peak-to-peak  separation  of  the  experimental  curves 
decreased,  and  the  peaks  became  sharper,  in  every  case,  as 
the  temperature  of  the  samples  was  increased.  The  experimen- 
tally  determined  value  of  l/T^  for  the  reference 

solution  was  291  sec  ^  at  25°  and  276  sec  ^  at  30°.  Glasel 
(1966)  determined  a  value  of  195  sec  ^  at  25° ,  while  Garrett 
et  al .  (1967)  determined  a  value  of  143  sec  ^  at  29°.  How¬ 
ever,  Hindman  et  al .  (1970)  have  performed  T^  measurements 

which  indicate  the  presence  of  paramagnetic  impurites  in 
the  solutions  of  Glasel  et  al .  and  Garrett.  A  correction  to 
their  data  would  decrease  their  values  of  l/T^ •  In  view  of 
this  finding,  the  results  of  this  work  would  seem  to  indicate 
the  presence  of  paramagnetic  impurities  or  instrumental 
broadening.  These  errors,  if  systematic,  should  then  cancel 
in  the  treatment  of  data  because  l/T^  for  the  blank  solution 
is  subtracted  from  1/T2  for  the  various  sample  solutions. 

There  is  a  fair  degree  of  scatter  in  the  points  of  the 
plots  of  log  (1/T2  ' )  vs.  1/T°K ;  however,  the  points  define 
lines  of  positive  slope.  This  would  indicate  that  Swift  and 
Connick's  (1962)  "Case  B"  or  "Case  D"  (in  the  case  of  the 
pH  7.92  solution)  is  being  observed.  For  Case  B, 


/ 


121 


Fig.  6-9.  A  schematic  derivative-of-absorption  curve. 
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T2p  =  PMTBA“B2  <16> 

If  Case  B  were  the  situation,  this  could  be  shown  by  a 

study  at  another  frequency.  Since  Aco„  varies  with  the 

13 

frequency,  a  doubling  of  the  frequency  would  produce  an 

increase  in  the  value  of  l/T^^  by  a  factor  of  four.  Since 

the  frequency  was  not  varied  in  this  study,  there  is  no 

kinetic  information  available  on  this  point  because  there 

is  no  observable  shift  caused  by  H2017,  which  would  allow  a 

2 

calculation  of  Au)n  ,  hence  .  In  addition,  P„„  would  have  to 

B  B  M 

be  known ,  which  involves  a  knowledge  of  the  coordination 
number  of  the  iron  atom  in  HRP  for  water.  This  is  presumed 
to  be  one,  possibly  two,  .the  other  four  positions  being 
occupied  by  the  nitrogen  atoms  of  the  porphyrin  ring. 

If  Case  D  were  the  situation,  then  no  kinetic  information 
concerning  tb  could  be  derived  from  the  associated  equation 


1 


2p 


?2B 


(20) 


A  calculation  could  be  made  of  1/T  in  this  case,  but  t  , 

2  C 

the  correlation  time,  is  not  known  for  HRP,  nor  is  the 
Fe-0  distance. 

Table  6-2  contains  values  for  the  activation  energy 
for  the  exchange  reaction,  determined  from  the  slopes  of  the 
six  plots.  There  is  no  trend  apparent;  the  average  of  all 
these  values  is  4.0  kcal  mole  This  fairly  small  value 

(compared  to  values  for  ferric  and  other  ions  found  by  Swift 
and  Connick  (1962)  and  Luz  and  Shulman  (1965)  of  8-12  kcal 
mole  ^)  may  be  due  to  outer  sphere  exchange  of  water. 
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Table  6-2 

Activation  Energies  as  a  Function  of  pH 


of  HRP  -  Solutions 

PH 

Activation  Energy 

(kcal  mole 

5.61 

6.8 

5.94 

4.1 

6.85 

5.2 

7.92 

0.8 

8.87 

3 . 0 

10. 15 

4.3 
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Table  6-3  contains  values  for  the  lower  limit  of  the 

rate  constant  for  the  exchange  process,  calculated  according 

to  the  method  of  Connick  and  Poulson  (1959) .  Again,  there 

is  no  trend  with  pH,  and  an  average  value  for  the  lower 

limit  of  k  is  about  3  x 106  sec”1  at  25°.  Connick  (1965) 

has  reported  a  rate  constant  for  exchange  of  water  molecules 

from  the  first  coordination  sphere  of  ferric  ion  at  25°  to  be 
3  -1 

3  x 10  sec  .  The  difference  in  values  may  be  due  to  an 
exchange  of  water  molecules  between  bulk  water  and  an  outer 
coordination  sphere  of  the  ferric  iron  of  HRP ,  an  easier 
process.  The  problem  remains  that  of  the  large  experimental 
broadening  of  resonance  curves  reported  in  this  work  which 
is  incompatible  with  previous  observations  (Glasel ,  1966; 
Garrett  et  al_.  ,  1967)  ,  and  which  may  be  reflected  in  serious 
errors  in  values  of  activation  energies  and  lower  limits 
for  rate  constants  for  the  exchange  process,  but  which  does 
not  affect  the  conclusion  that  the  true  values  for  the  rate 
of  exchange  are  not  obtainable  in  this  temperature  region 
by  the  NMR  line  broadening  technique. 
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Table  6-3 

Lower  Limits 

of  the  Rate  Constant  for 

Exchange 

as  a  Function  of  pH 

PH 

kn  ,  -lx 

lower  (sec  ) 

5.61 

2 . 8  x  io5 

5.94 

8 .2  x 106 

6.85 

1.9  x io5 

7.92 

1.2  x io5 

8.87 

5.8  x io6 

10. 15 

4.1  x io6 

••• 
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Summary 

The  work  on  lactoperoxidase  reported  in  this  thesis 
contains  the  first  detailed  studies  on  individual  rate 
constants  in  the  lactoperoxidase  oxidation-reduction  cycle . 

A  measure  of  the  secondary  structure  of  lactoperoxidase 
has  been  established;  it  was  shown  that  lactoperoxidase 
contains  approximately  17%  a-helix  at  pH  7 . 0 .  In  addition, 
it  was  shown  that  anions  such  as  fluoride,  cyanide,  and 
azide  affect  the  ORD  spectrum  of  the  native  enzyme,  indicating 
that  these  ions  may  cause  an  alteration  of  the  geometry  of 
the  heme  group  with  respect  to  the  enzyme. 

A  detailed  pH  study  of  the  kinetics  of  the  formation 
of  lactoperoxidase  compound  I  over  a  wide  pH  range  was 
made,  and  the  pH-independent  rate  constant  is  (9.2  ±0.9)  x 
10^  M  ^sec  .  In  addition,  it  was  demonstrated  that  formate 
buffer  binds  to  lactoperoxidase  at  an  ionic  strength  of 
0.05. 

The  reaction  of  LP-II  with  two  substrates,  p-cresol  and 
iodide  ion,  was  studied  over  a  wide  pH  range.  It  was  shown 
in  the  case  of  iodide  ion  oxidation  that  a  rate  controlling 
ionization  on  the  enzyme  occurs  outside  the  experimental  pH 
range  at  low  pH.  In  addition,  a  second  order  dependence 
on  the  concentration  of  iodide  ion  was  noted,  and  a  complex 
between  LP-II  and  iodide  ion,  which  reacts  further  with 
iodide  ion,  was  proposed  to  account  for  this  observation. 

In  the  case  of  the  oxidation  of  p-cresol  by  LP-II,  it  was 
shown  that  the  rate  is  dependent  on  three  ionizations  of 


' 
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groups  in  the  enzyme  with  pK  values  of  2.3,  5.8,  and  9.7, 

a 

and  the  ionization  of  the  substrate.  It  was  also  demon¬ 
strated  that  the  neutral  form  of  p-cresol  was  the  reactive 
species  over  the  whole  pH  range. 

An  investigation  of  the  rate  of  exchange  of  water  with 

the  iron  atom  of  the  heme  group  of  horseradish  peroxidase 

17 

was  attempted  using  H2O  NMR.  Experimental  difficulties 
precluded  obtaining  accurate  values  of  line  widths,  but  the 
conclusion  was  made  that  the  true  values  for  the  rate  of 
exchange  could  not  be  obtained  in  the  experimental  tempera¬ 
ture  range . 

The  work  on  lactoperoxidase  reported  in  this  thesis 
provides  the  basis  for  further  studies.  For  example,  the 
kinetics  of  the  transition  from  LP-I  to  LP-II  using  various 
reducing  agents  has  yet  to  be  studied.  Preliminary  experi¬ 
ments  have  shown  several  experimental  difficulties.  The 
spontaneous  rate  of  decomposition  of  LP-I  at  pH  7  is  associ¬ 
ated  with  a  rate  constant  of  4  sec  ^  which  apparently 
increases  as  the  pH  is  either  raised  or  lowered.  The 
spontaneous  rate  also  varies  considerably  with  the  state  of 
purity  of  the  enzyme.  A  study  of  the  oxidation  of  iodide 
ion  by  LP-I  would  be  important  in  view  of  the  results  of 
Bjorksten  (1970)  and  Roman  and  Dunford  (submitted  for 
publication) ,  which  show  that  in  the  oxidation  of  iodide 
ion,  HRP-I  goes  directly  to  HRP  in  a  two-electron  transfer, 
without  the  appearance  of  HRP-II  in  the  reaction  cycle. 


-  , 
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Another  problem  of  interest  would  be  the  simulation 
of  the  action  of  thyroid  peroxidase,  i.e.,  an  investigation 
of  a  system  in  which  LP  catalyzes  the  oxidation  of  iodide 
ion  and  the  iodination  of  tyrosine.  However,  this  would 
require  studies  of  the  elementary  processes,  that  is,  the 
oxidation  of  iodide  ion  and  tyrosine  by  the  individual 
compounds  I  and  II  before  a  steady  state  reaction  is  under¬ 
taken.  Prior  knowledge  of  all  of  the  individual  rate 
constants  would  simplify  evaluation  of  the  steady  state 
results . 

As  far  as  the  nature  of  LP-II  is  concerned,  the 
difference  in  the  log  k-pH  profiles  for  the  oxidation  of 
iodide  ion  and  p-cresol  has  been  shown  in  this  thesis.  This 
difference  was  suggested  to  be  due  to  either  the  absence  of 
steric  interactions  in  the  case  of  the  oxidation  of  iodide 
ion,  or  to  the  operation  of  a  different  mechanism.  Further 
work  with  other  substrates  could  shed  light  on  the  structure- 
reactivity  relationship  of  the  active  site  of  LP-II.  For 
example,  an  homologous  series  of  phenols  could  be  used  as 
substrates,  with  substituents  on  the  benzene  ring  that  vary 
in  size  and  electronegativity.  On  the  other  hand,  other 
inorganic  ions  could  be  used  as  substrates,  such  as  ferro- 
cyanide ,  nitrite,  and  sulfite  ions.  The  role  that  lacto- 
peroxidase  plays  in  vivo  is  not  known,  so  that  even  if  the 
substrates  used  for  lactoperoxidase-catalyzed  oxidations 
are  of  little  physiological  importance,  a  detailed  study  of 
their  reactions  may  provide  a  meaningful  picture  of  the 
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mechanism  of  oxidation  by  lactoperoxidase ,  and  to  a  certain 
extent  by  peroxidases  in  general. 
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Appendix  1 

Purification  of  Lactoperoxidase 

Lactoperoxidase  was  supplied  by  Dr.  Martin  Morrison 
of  St.  Jude  Children's  Research  Hospital,  Memphis,  Tennessee, 
and  by  Calbiochem,  Los  Angeles.  The  lactoperoxidase  (LP) 
samples  from  both  sources  were  obtained  according  to  the 
method  of  Morrison  and  Hultquist  (1963) ,  described  below. 

All  purification  procedures  were  performed  at  4°. 

Fifteen  grams  of  the  sodium  form  of  carboxymethyl  cellulose 
cation  exchange  resin  (Bio-Rad)  were  added  to  each  liter  of 
raw  skim  milk,  and  the  pH  was  adjusted  to  7.0.  The  milk-resin 
suspension  was  stirred  for  1  hour  and  then  allowed  to  stand 
until  the  resin  settled.  The  milk  was  then  removed  by 
decantation,  and  an  additional  7.5  g  of  resin  per  liter  was 
added.  The  pH  was  readjusted  to  7.0  and  the  milk-resin 
suspension  was  stirred  again  for  1  hour  and  then  allowed  to 
settle.  The  combined  resin  was  washed  with  distilled  water 
in  batches  by  alternately  suspending  the  resin  and  allowing 
it  to  settle. 

The  resin  was  then  washed  into  a  chromatographic  tube 
with  distilled  water.  The  LP  was  eluted  with  0.5  M  sodium 
acetate.  The  enzyme  moved  down  the  column  as  a  dark  green 
band,  increasing  in  size  as  it  descended.  The  green  material 
which  eluted  was  precipitated  with  ammonium  sulfate  and 
centrifuged,  dissolved  in  a  small  volume  of  water,  and 
dialyzed  repeatedly  against  large  volumes  of  distilled  water. 
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The  resulting  solution  was  then  applied  to  a  column  of 
Sephadex  G-200  gel.  The  LP  fractions  collected  exhibited 
P.N.  >0.8  (P.N.  is  the  ratio  of  the  absorbance  of  an  LP 
solution  at  412  nm  to  that  at  280  nm) .  The  use  of  P.N.  was 
the  criterion  of  purity  of  samples  of  lactoperoxidase  used 
throughout  the  course  of  this  research.  A  P.N.  of  0.8  or 
greater  signified  that  the  sample  was  relatively  pure. 

P.N.  values  exceeding  0.95  for  LP  are  rare,  and  the  suspicion 
arises  that  extensive  chromatographic  purification  procedures 
cause  the  enzyme  to  lose  activity.  An  attempt  to  isolate  LP 
in  this  laboratory  using  the  method  outlined  above  yielded 
63  mg  LP  (P.N.  0.7)  from  5  gallons  of  milk. 

One  of  the  more  common  assay  methods  for  the  activity 
of  lactoperoxidase  samples  involves  the  oxidation  of 
guaiacol  (Maehly  and  Chance,  1958) .  The  major  disadvantage 
of  this  procedure  is  that  the  color  of  the  product  fades 
rapidly,  which  may  be  responsible  for  the  lack  of  reprod¬ 
ucibility  of  this  procedure.  It  appears  that  despite  the 
proliferation  of  activity  tests  for  peroxidases,  the 
absorbance  ratio  P.N.  is  the  most  widely  used  criterion  of 
purity. 

Further  purification  of  LP  obtained  from  Dr.  Morrison 
or  Calbiochem  was  carried  out,  and  the  procedure  will  now 
be  described.  The  LP  from  either  source,  which  was  obtained 
in  a  lyophilized  (freeze-dried)  form,  possessed  a  P.N.  of 
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about  0.6.  It  was  purified  to  a  P.N.  of  0.8  or  greater  by 
Sephadex  gel  chromatography. 

Sephadex  is  a  modified  dextran .  The  dextran  macromolecules 
are  cross-linked  to  give  a  three-dimensional  network  of 
polysaccharide  chains.  A  gel  filtration  experiment  can  be 
described  schematically  in  the  following  way.  Molecules 
larger  than  the  largest  pores  of  the  swollen  Sephadex,  i.e., 
above  the  exclusion  limit,  cannot  penetrate  the  gel  particles 
and  therefore  they  pass  through  the  bed  in  the  liquid  phase 
outside  the  Sephadex  particles.  They  are  thus  eluted  first. 
Smaller  molecules,  however,  penetrate  the  gel  particles  to 
varying  degrees  depending  on  their  size  and  shape.  Molecules 
are  therefore  eluted  from  a  Sephadex  bed  in  the  order  of 
decreasing  molecular  size.  As  the  molecules  pass  through 
the  bed  at  different  rates,  they  emerge  at  the  outlet  end  of 
the  column  separated  from  each  other.  The  zones  eluted  are 
somewhat  wider  than  the  applied  sample,  but  by  the  use  of 
suitable  experimental  conditions,  this  zone  broadening  can 
be  limited.  Sephadex  gels  and  chromatographic  columns  are 
obtained  from  Pharmacia  Fine  Chemicals,  Inc. 

All  procedures  described  were  carried  out  in  a  cold  room 
at  4°,  unless  otherwise  specified.  The  gel,  Sephadex  G-200 
Superfine,  was  swollen  by  suspension  in  pH  7.0  phosphate 
solution  of  ionic  strength  0.05  for  3  days.  A  chromatographic 
column  was  prepared  by  pouring  the  gel  into  a  glass  column 
of  dimensions  45  cm  in  length  and  2.5  cm  in  diameter, 
equipped  with  a  net  at  the  bottom  which  allowed  the  passage 
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of  water  and  enzyme,  but  not  the  gel  particles.  The  quality 
of  the  gel  packing  was  improved  by  the  passage  of  1  litre  of 
pH  7.0  buffer  through  the  column.  This  procedure  took  about 
2  days,  and  the  column  was  then  ready  for  application  of  the 
enzyme.  About  150  mg  of  lyophilized  LP  was  dissolved  in 
2  ml  of  pH  7.0  phosphate  buffer  of  ionic  strength  0.05.  The 
solution  was  clarified  by  passage  through  glass  wool. 

The  enzyme  solution  was  applied  to  the  top  of  the 

Sephadex  column  with  a  disposable  pipet.  After  the  green 

band  had  descended  under  the  top  of  the  gel  column,  a  buffer 

from  a  gravity  bottle  was  allowed  to  drain  into  the  column. 

The  eluant  was  pH  7.0  phosphate  buffer  of  ionic  strength 

0.05.  The  output  of  such  a  column  was  about  10  ml  per  hour, 

and  fractions  were  collected  every  10  min  (about  1 . 5  ml  per 

fraction)  in  test  tubes  with  an  LKB  Ultrorac  7000  fraction 

collector.  After  the  enzyme  had  eluted  from  the  column, 

the  test  tubes  containing  the  green  enzyme  (or  brownish- 

green  color  at  high  LP  concentrations,  about  10~^M)  were 

collected,  and  their  contents  analyzed  spectrcphotometrically 

(using  a  Beckman  DU  at  25°)  at  412  nm  and  280  nm.  The 

tubes  were  then  stored  in  the  cold  prior  to  use.  The  results 

of  one  such  purification  are  shown  in  Table  A-l ,  where  the 

absorbance  values  are  of  samples  of  1  mm  thickness,  and  the 

LP  concentrations  are  obtained  using  the  molar  absorptivity 

5  -1  -1 

£^22  =  1*14  x  10  M  cm  (Morrison  et_  a^L.  ,  1957)  .  Figure 
A-l  is  a  plot  of  P.N.  (from  Table  A-l)  vs_.  fraction  number. 
Cnly  those  samples  of  LP  of  P.N.  greater  than  0.8  were  used 
in  kinetic  experiments. 


Table  A-l:  The  results  of  a  chromatographic  purification  of 


lactoperoxidase  using  Sephadex  G-200  Superfine  gel. 


Fraction  # 

a 

A412 

a 

A280 

P.N. 

[LP]  ,  M 

73 

0.06 

0.10 

0.65 

5.70  x 10“6 

74 

0.13 

0.17 

0.77 

1.14  x io-5 

75 

0.24 

0.31 

0 . 80 

2 .10  x io“5 

76 

0.46 

0.53 

0.87 

4.03  x io“5 

77 

0.74 

0.84 

0.88 

6.49  x io"5 

78 

1.06 

1.22 

0.87 

9.30  x io-5 

79 

1.31 

1.50 

0.87 

1.15  x io~4 

80 

1.60 

1.86 

0.86 

1.40  x io"4 

81 

1.90 

2 . 00 

0.95 

1.67  x io"4 

82 

1.95 

2.14 

0.91 

1.71  x io"4 

83 

1.54 

1.68 

0 . 92 

1.35  x io-4 

84 

1.25 

1.37 

0 . 91 

1.10  x io-4 

85 

0.85 

0.93 

0.91 

LO 

1 

o 
1- 1 

X 

VO 

• 

r- 

86 

0.54 

0.60 

0.91 

4 .74  x  io"5 

87 

0.34 

0.38 

0.89 

2 .98  x  io"5 

88 

0.21 

0.24 

0.88 

1.84  x  io"5 

89 

0.13 

0.16 

0.86 

1.14  x  io"5 

90 

0.09 

0.12 

0 . 75 

VO 

1 

o 

1 — 1 

X 

o\ 

00 

• 

91 

0.06 

0.08 

'0.73 

5.26  x  10  6 

92 

0 . 04 

0.05 

0.73 

3.51 x  io”6 

Absorbance 

values  are 

obtained 

using  1 

mm  cells. 

/ 
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Fig.  A-l.  A  plot  of  P.N.  vs.  fraction  number  for  the  results  of  Table 
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Appendix  2 

Description  of  the  Non-Linear  Least  Squares  Program 

This  computer  program,  called  "BRT2"  was  supplied  by 
Dr.  M.  K.  Evett  (N.R.C.C.  post-doctoral  fellow  at  the 
University  of  Alberta,  1968-1970).  The  program  uses 
Gauss's  method  to  find  the  best  least-squares  values  of 
parameters  in  non-linear  functions  of  one  dependent  variable. 
A  discussion  of  the  minimization  methods  used  is  given  by 
Wentworth  (1965)  . 

The  program  is  listed  at  the  end  of  this  appendix. 

It  contains  some  features  that  are  not  often  used  (e.g. 
CNTROl(ll)  and  (12)  discussed  below).  The  subprograms 
OVERID,  WRNGWY ,  and  PCHANGE  may  be  removed  with  their  calls, 
or  made  dummies  to  reduce  space  requirements.  In  addition, 
the  "restrained  change  mechanism"  in  subroutine  GAUS  can  be 
removed . 

There  are  three  user-supplied  subprograms  for  BRT2  that 
must  be  programmed  for  use  on  a  particular  equation  that  one 
may  wish  to  solve.  The  first  subprogram,  called  MANIP, 
allows  one  to  modify  the  input  variables.  The  second  sub¬ 
program,  YCOMP ,  calculates  a  value  of  the  dependent  variable 
as  a  function  of  the  parameters  and  independent  variables. 

The  third  subprogram,  OUTDAT,  is  called,  if  needed,  at  the 
completion  of  the  "fit"  to  allow  computations  with  the 
parameters  determined  by  YCOMP.  These  three  subprograms  are 
found  at  the  end  of  the  listing. 

Examples  of  the  use  of  these  three  subprograms  will  be 
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given  : 

MANIP  —  if  the  data  were  obtained  as,  for  example,  log  k 
vs .  pH,  and  the  actual  equation  related  log  k  with  [H+] , 
then  MANIP  could  be  used  to  convert  the  input  variable  pH 
to  [H+]  by  the  following  cards 

DO  1  1=1 , NPTS 

1  Z (2 , 1)  =  EXP (-2 . 303*Z  (2 ,1) ) 

which  is  equivalent  to  saying  that  [H+]  =  e  2,303  PH ,  and 
performing  this  operation  for  as  many  input  points  as  is 
required . 

YCOMP  —  if  one  wished  to  program  a  first-order  kinetic 
expression  involving  absorbance,  A,  a  rate  constant,  k,  and 
time,  t,  which  are  related  by 


A 


oo 


A 

oo 


-kt 

e 


then  the  expression  is 

YCOMP  =  (B  (2) -B  (3)  )  *EXP  (-B  (1)  *Z  (2  ,N)  ) +B  (3) 
where 

YCOMP  =  A,  the  dependent  variable 

Z(2,N)  =  t,  the  independent  variable 

B  ( 1 )  =  k 

B  ( 2 )  =  A0 

B  ( 3 )  =  A 

00 

OUTDAT  —  if,  for  example,  one  obtained  from  YCOMP  a  value 
of  k  which  was  a  pseudo  first  order  rate  constant,  but 
wished  to  obtain  the  value  of  the  true  second  order  rate 
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constant  by  the  relation 

k*  =  k/[S] 

where  [S]  is  a  concentration  of  some  reactant,  then  the 
expression  would  merely  be 


K  =  B ( 1 ) /A 


with  an  accompanying  PRINT  statement. 

The  program  listing  may  be  consulted  for  specification 
cards.  The  subprograms  must  be  present  even  if  only  as 
dummies  (with  RETURN  as  the  only  executable  statement) . 

The  following  is  a  list  of  the  arrangement  of  data 
cards  for  the  BRT2  program: 


Card  No. 


Data 


Format 


1 

2 

3 


CNTROL ( I )  1=3,19  16,1614 

CNTROL ( 3 ) ,  TITLE  I6,18A4,2A1 

Number  of  parameters,  NC  15 

Number  of  parameters  varied,  NB  15 
Parameter  estimates, 

B(J)  J=1,NC  7E1 0 . 4/  ( 8E10 . 4 ) 


The  actual  data  cards  follow  card  3,  and  the  data  is  right 
justified  (E10.4)  with  the  dependent  variable,  Z(1,N),  in 
col.  1-10  and  the  independent  variable,  Z(2,N),  in  col.  11-20. 
At  the  end  of  all  the  data  cards  is  a  card  with  a  zero  in 
the  sixth  column.  Note  that  the  program  will  vary  the  first 
NB  of  the  NC  parameters. 

The  following  is  an  explanation  of  the  controls  used  by 
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the  BRT2  program: 


CNTROL ( 1 ) 

Number  of  parameters. 

CNTROL ( 2 ) 

Number  of  parameters  varied. 

CNTROL (3) 

Number  of  data  points.  On  CNTROL  card  make 

this  number  larger  than  the  largest  data  set. 

Use  the  exact  number  on  card  2;  this  number 

may  be  omitted  from  card  2  if  each  data  set  is 

terminated  by  a  card  with  999.  in  columns  7-10 

CNTROL (4) 

Number  of  parameters . 

CNTROL (5) 

Total  number  of  variables  (independent  +1) . 

CNTROL (6) 

Limit  on  the  number  of  iterations. 

CNTROL (7) 

=a+2b+4c+8d  where  a,b,c,  and  d  are  each  1  or  0 

according  to  whether  or  not  the  corresponding 

printout  is  desired. 

a.  data 

b.  sum  of  squares  and  parameters  for  each 

iteration . 

c.  normal  equation  matrices. 

d.  restrained  change  mechanism. 

CNTROL ( 8 ) 

Weighting  factor.  1)1.0,  2)l/obs,  3) 1/calc., 

4)obs,  5)calc,  6)1/Z(M,N),  7)Z(M,N),  where  in 

the  last  two  the  weight  is  read  in  as  the  last 

independent  variable  in  columns  21-30. 

CNTROL (9) 

Formerly,  OUTDAT  was  activated  if  the  value  of 

this  control  was  non-zero.  OUTDAT  is  now 

controlled  through  MANIP • 
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CNTROL (10) 

Corrections  to  parameters  are  scaled  by 

0 . 1 *CNTROL (10) .  If  CNTROL (10)  is  equal  to 

zero,  the  scaling  factor  is  1. 

CNTROL (11) 

Parameter  change  option.  Reduces  the  number 

of  parameters  by  one  this  many  times  when  no 

reduction  in  the  sum  of  squares  is  obtained. 

CNTROL (12) 

If  >0,  upon  termination  -  this  many  "over¬ 
rides"  (i.e.,  the  program  ignores  the  fact  that 

there  is  no  reduction  in  the  sum  of  squares) . 

If  <0,  the  program  performs  this  many  "wrong- 

ways"  (i.e.,  the  parameters  are  scaled  and 

changed  in  the  opposite  direction) . 

CNTROL (16) 

Allows  this  many  data  sets  to  be  run  under  one 

control  card.  May  be  -1  for  automatic  counting 

if  sets  are  terminated  by  a  card  with  STOP 

punched  in  columns  7-10. 

CNTROL (17) 

This  many  sets  of  starting  estimates  are  used. 

Controls  13,14,15,18,19,  and  higher  are  internal  program 
controls . 

This  program  has  been  tested  against  a  program  used 
previously  in  this  group  (Ellis,  1968)  on  a  number  of 
problems,  and  has  been  found  to  give  identical  results;  how¬ 
ever,  the  BRT2  program  has  the  advantages  of  several 
different  forms  of  weighting  of  data  points,  as  can  be  seen 
from  "CNTROL ( 8 ) " ,  and  also  the  advantage  that  each  data 
point  can  be  individually  weighted.  The  BRT2  program  is 
listed  on  the  following  pages. 


' 


iSIGNON  HAL  T= 
ON  AT  15:13.37 
SLIST  4SOURCE* 

5  C 

6 
7 

a 

9 

10 

it 

12 

13 

14 

15 

16 

1  7  1 

IS 

19 

20 
21 
22 

23 

24 

25 

26  2 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37  CALL 

38 

39 

40 

41 

42 

43 

44  3 

45 

46 

4  7  4 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60  5 

61 


1* , '/OBS* . • 
•  (M+l  •  .  •  ,N> 


CNTROL (161=1 
CNTROL! 171=1 


1=1 ,201 


1M  P=1 00  PRIO=H  'JIM  MAGUIRE' 

ON  01-14-72  LAST  ON  AT  15:11.05  ON  01-09-72 

PROGRAM  BRT2 
DIMENSION  BTRIAL(201 
DIMENSION  COR  E ( 20  00 1 
COMMON  CNTROL 
DIMENSION  NPR I  NT (41 
INTEGER  I  PR  1  NT ( 2 1 / • NO •  .  • YES • / 

INTEGER  CNTROL ( 25 1 . 00  I M, T 1 T( 20 1 
INTEGER  ST/'STOP'/ 

INTEGER  I  FORM (3,10  1/'  ','1.0','  ',' 

1'  '.'OBS'.'  ','  *  , ' CALC  *  , '  ',  *1/2', 

19*'  • / 

BIN(K,J)=0. 5—0 . 5* (  —  1 .  )**(K/2**(J-1 1  ) 

READ  12,  (CNTROLt  I  1  . 1=3.  191 
IF  (CNTROL( 3 1 .EO.O 1  GO  TO  It 
KNTRL=CNTROL (101 
MDI M=CNTROL( 51 
NDI M= I ABS ( CNTROL l 3 1 1 
IF  { CNTROL ( 161. EO.O) 

IF  (CNTROL! 1 7 1 .EQ .0 1 
I  I  I=CNTROL( 16) 

NPRO=0 
NPRO=NPRO+l 
CNTROL ( 23 1 =NPRO 
READ  13,  CNTR0L(3) ,  (TIT( I  1 
IF  (T IT( 1 1 .EO.ST 1  GO  TO  1 
PRINT  14,  TIT 
LAST=1 

NPROB=CNTROL( 17) 

PRINT  15 

DO  4  L=1 , NPROB 

READ  16.  IB,  IC,  (BTRIAL(  I  1  ,1  =  1  , IB) 

PRINT  17.  L. IB,  IC , ( BTR  I  AL (  I  1  ,  I  =  1  , IB  1 
TO  MANIP  ALLOWS  PARAMETERS  TO  BE  MANUPULATED. 
CALL  MANIP  (  1  ,BTR I AL, CORE (LOC2 1 , MDI M,ND I M,  BDI Ml 
1  =  0 

L  A  STP2=L  A  ST 4- 2 
NEXT=LASTP2+IB- 1 
DO  3  K=LASTP2 .NEXT 
1  =  1  +  1 

CORE(K)  =  STRIAL( I  1 
CORE (LAST  1  =  10 
CORE(LAST+l 1 =IC 
L AST=NEXT+ 1 
BD I M=CNTPOL ( 4 1 
L0CZ=4  *BDI M  +  L AST 
LOC  A=L  O  C  Z  +  MD  I  M*Nl)  IM  +  1 
NEXT  = 1 

DO  10  K= 1 , NPROB 
CNTROL ( l 0  >  =  KNTRL 
CNTROL! ?4)=K 
CNTROL  (2 1  =  CORE! NEXT  1 
CNTROL! 1  1 =COR  E { NEXT  +  1  1 
M3  =  C  N  T  ROL ( 2 1 
LOCB=L AST 
DO  5  I =1 . M3 

COREILAST+I-l )=CORE(NEXT+l+I 1 
NEXT=NLXT+M3+2 


BRT20000 

BRT20010 

BRT20020 

BRT20030 

BRT20040 

BRT20050 

BRT20060 

BRT20070 

,  '  1 /'  .  •  CALC •  , •  '.BRT20080 

•21  • .  • M+l ,  •  ,  ' N 1  • , BRT  20090 
BRT20100 
BRT20 1 1 0 
BRT  20120 
BRT  201 30 
BRT20140 
BRT201 50 
BRT20160 
BRT  20 170 
BRT20 1 80 
BRT201 90 
BRT  20  200 
BRT202 1 0 
BRT20220 
BRT20230 
BRT20240 
BRT 20  250 
BRT20260 
BRT20270 
BRT  202  30 
BRT20290 
BRT20300 
BRT20310 
BRT20320 
BRT 2033  0 
BRT  20340 
BRT20350 
BRT20360 
BRT20370 
BRT20380 
BPT20390 
BRT20400 
BRT204 10 
BRT20420 
BRT20430 
BRT20440 
BRT  20450 
BRT20460 
BRT20470 
BRT20480 
BRT  20490 
BRT20500 
BRT2051 0 
BR  T  20520 
BRT  20530 
BRT2054  0 
BRT20550 
BRT  20560 


62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

1  00 

101 

1  02 

1  03 

1  04 

105 

1  06 

1  07 

108 

1  09 

I  10 

1  1  1 

1  12 

1  1  3 

1  1  4 

1  15 

1  16 

1  1  7 

l  1  8 

1  I  9 

1  20 

1  21 
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CALL 


10 


1 1 


=1 .19), N8.NPTS.NVARS, LIMIT 
TO  6 


NB  =  CNT  ROL (  1  ) 

NPT S  = I ABS< CNT ROL ( 3 ) ( 

NVARS=CNTROL (5) 

L I M I T=CNTROL ( 6 ) 

PRINT  13.  (CNTROL(I).I 
IF  ( CNTROL ( 7 ) . NE . 0 )  CO 
PRINT  19 
GO  TO  8 
PRINT  20 
00  7  J=1 ,4 
L=  1 

IF  ( BI N< CNTROLI 7 ) , J ) 

NPRINT (J)=IPRINT (  L ) 

CONT I NUE 

PRINT  21,  (NPR INT ( J ) , J= 1 
LL  =  CNTROL<  8) 

PRINT  22.  CNTR0LI8) » ( IFORM( I ,LL) , 1=1 .3) 

IF  <  CNTROL (lO).EO.O)  PRINT  23 
ICHANG=1 O+CNTROLI 101 

IF  ( CNTROLI 1 0 ). NE. 0 )  PRINT  24,  ICHANG 

PRINT  25 

LASPM3=LAST+M3-l 

PRINT  26,  (CORE! I ) . I =LA ST, LASPM3 ) 

CNTROL ( 22) = 1 

PRINT  27,  MO  I M , NOI M, BOI M 

TO  MANIP  ALLOWS  DATA  TO  BE  READ  IN  AND  MODIFIED  IN  THAT 
IF  I CNTROLI 3 )  . LT . 0 . AND.K . EO. 1  )  CALL  MANIP  I  2 . CORE  I  LAST ) 
1  ,  MD  I  M, NDI M, BOI M ) 

IF  I NPTS. NE.  I A3S( CNTROLI 3  )  >  )  NPT S= I ABS I CNTROL l 3) ) 


>  G  T , 0 •  )  L  =  2 


>4) 


BRT20570 
BRT20580 
OR  T 2  059  0 
BRT20600 
BRT20610 
BRT20620 
BRT20630 
BRT2064  0 
BRT20650 
BRT20660 
BRT20670 
BRT20680 
BR T 2  069  0 
BR  T  2  0700 
BRT2O710 
8RT20720 
BRT20730 
BRT20740 
BRT20750 
BRT20760 
BRT20770 
BRT20780 
BRT20790 
BRT20800 
BR  T  20  8 1 0 
SUBROUT I NE.BRT20820 
, CORE! LOCZ ) BRT20830 
BRT20840 
BR  T2  085  0 


AND.K.EQ.l)  CALL  INDATA 


(CORE  I LOCB ) . CORE!  LOCZ ) . BRT2 0  86  0 

BRT20870 


IF  I CNTROLI 3 ) .GE. 0 
1 MDI M, NDI M, BDI M) 

IF  ( BI N( CNTROLI 7 ) . 2 ) .EO. 0 )  GO  TO  9 

CALL  PRINDA  (CNTROL , CORE! LOCZ ), MD IM,  NDI  M ) 

CALL  GAUSS  ( COR E{ LOC B ) . CORE! LOCZ » « M D I M, NDI M , BD 1 M , CORE ( LOC A ) ) 
PRINT  14,  TIT 

CALL  FINALE  I  CORE I LOCB )  , C ORE  I L OC Z ) , MD I M , NO  I M , BD I M . CORE  I LOC A)  ) 
IF  (CNTROLI21 ) .E0.0)  PRINT  33 

28 

29 

30 

31 

32 


( CNTROL (21 ) .EQ. 1 ) 

I CNTROL (21 ) • EQ . 2 ) 

( CNTROL (21  )  , 

(CNTROL ( 21  )  , 

( CNTROL (21  )  , 

CONT I NUE 

IF  ( NPRO.LT. I ABS(  I  I  I )  .OR. I  I  I .EQ.-l )  GO  TO  2 
GO  TO  1 
PRINT  34 


IF 
IF 
IF 
I  F 
IF 


,  EQ. 
.EQ. 
.  EQ  . 


3) 
4  ) 
5) 


PR  I  NT 
PRINT 
PRINT 
PR  I  NT 
PR  I  NT 


BRT20880 
BRT  20890 
BRT20900 
BRT20910 
BRT20920 
BRT  20930 
BRT20940 
BRT20950 
BRT20960 
BRT  20970 
BRT  20980 
BRT20990 
BRT21000 
BRT21010 
BRT21 020 


STOP 

BRT  2 1030 

c 

BRT21 04  0 

12 

FORMAT 

(16.1614) 

BRT21 050 

13 

FORMAT 

(16, 18A4 , 2 A  1  ) 

BRT2 1 060 

14 

FORMAT 

(  1H1 , 2  0  A  4 ) 

BRT2I070 

IS  FORMAT! 

lHO.'BTRIAL  READ-IN,  ECHO 

CHECK  FOLLOWS.' 

> 

BRT 2 1 060 

1  6 

FORMAT 

(  2  I  5 . 7E  1  0 . 4  /  (  BE  1  0 . 4  )  ) 

BRT2  l  090 

17 

FORMAT 

( 3 I 5 , 7E 1 5 . 5/ ( I5X.7E15.5)  ) 

BRT  2 1  100 

1  8 

format 

( 32H0  NONLINEAR  REGRESSION 

CONTROLS  19I5//30H 

NUMBER 

OF  PARBRT2 I  1  1 0 

1  AMETF.RS 

I12/30H  NUMBER  OF 

DATA  POINTS 

I  12/3  OH 

NUMBEBWT2 I 120 

2R  OF  VARIABLES  I12/30H  LIMIT  ON  NUM BER/ I  TCP AT  I ONS  I  I  2/ ) BRT2 1  1 30 

19  FOPMAT  (•  CNT  ROL ( 7  )  =  0.  SO  NO  INTERMEDIATE  OUTPUT  •)  BRT21140 

20  FORMAT  ('0  THE  FOLLOWING  INFORMATION  IS  TO  BE  PRINTED  OUT'//*  INPUBRT21I50 
IT  DATA  SUMSQ  AND  PARAMETER  VALUES  NORMAL  EQ  MATRICES  REST RA 1 BRT2 1 1 60 


1  22 

123 

124 

1  25 

1  26 

127 

1  28 

1  29 

1  30 

131 

1  32 

1  33 

134 

1  35 

1  36 

1  37 

1  38 

1  39 

140 

1  4  1 

1  42 

143 

144 

1  45 

1  46 

1  4  7 

148 

1  49 

1  50 

151 

152 

1  53 

1  54 

1  55 

1  56 

1  57 

1  58 

1  59 

1  60 

161 

1  62 

1  63 

1  64 

1  65 

1  66 

1  67 

1  68 

1  69 

1  70 

171 

172 

1  73 

1  74 

1  75 

1  76 

1  77 

1  78 

1  79 

180 

1  81 


151 


21 


25 

26 
27 


2NED  CHANGE  MECHANISM*  1 
FORMAT  (5X.AJ.18X.A3.21X.A3.25X.A3) 

22  FORMAT  (•  CNT  ROL ( 8 )  =*.13.’  SO  MINIMIZE  LEAST  SQUARES  WEIGHTED  BY 

1  ( • ,3A4 , •  ) • ) 

23  FORMAT  (•  CNTROL(IO)  =  O,  GO  UNC ONS T R A I  NED  PARAMETER  CHANGES  •) 

24  FORMAT  (  •  CNT  ROL (10)  NONZERO.  PARAMETERS  CAN  CHANGE  BY  AT  MOST 
1  *  .  I  4 .  •  PERCENT*  ) 

FORMAT  (42H0  STARTING  GUESSES  FOR  THE  PARAMETERS  D<J)> 

FORMAT  ( 1 X. 1 OE 1 3. 4 ) 

FORMAT  I  4  8H  0  NONLINEAR  REGRESSION  ARRAY  DECLARATIONS  ARE  /4H 
I  12. 1H.I  3.8H)  ANO  B(  12. 1H1/1H0) 

FORMAT  (*0  NORMAL  TERMINATION.  CONVERGENCE  OF  PARAMETERS  •) 

(*0  TERMINATION  DUE  TO  NfGLlGIBLE  CHANGE  IN  SUM/SQUARES 
(  '  0  TERMINATION  DUE  TO  LIMIT  ON  NJMBER  OF  ITERATIONS  •) 
(*0  TERMINATION  DUE  TO  SINGULARITY  OF  NORMAL  EQS  •) 


28 

29 

30 

31 

32 

33 

34 


FORMAT 
F  ORMAT 
FORMA  T 
FORMAT  ( ' 
FORMAT  ( 


BRT 2  I  1 70 
ORT21  I  80 
BRT21 190 
BRT21 200 
OR  12  12  1 0 
HRT21220 
BRT21230 
BRT2 1240 
BRT21250 
Z ( ORT2 1 260 
BRT21 270 
BRT2. 1280 
•  ) BRT  2 1 290 
8RT2 1300 
BRT21310 


OF  SQUARES  REDUCTION  IN  LEAST 
REDUNDANCY  -  -  -  EXPRESSION 
JOB  TERMINATION.') 


SQUARES  I NTER VAL* ) 8RT2 1 320 
CHANGE  INDICATED* (BRT21330 
BRT21340 


I  Q.Z.MDIM.NDI M.BDI  M) 


C 

c 

c 

c 

c 


0  NO  SUM 
0  OPTION 
FORMAT!' ISTANDARD 
END 

SUBROUTINE  INDATA 
COMMON  CNTHOL ( 25 ) 

DIMENSION  9(1),  Z(MOIM.NOIM) 

INTEGER  CNTROL.BDIM 
NVARS=CNTROL ( 5  I 
I BSET=CNTROL( 24 ) 

KTROL=CNTROL ( 14) +1 
M=CNTKOL ( 15) + 1 
IDSET=CNTROL(23) 

NPTS=CNTROL( 3 ) 

IF  ( CNTROL( 3 ) . EQ.O )  GO  TO  8 
GO  TO  (1.3.7),  KTROL 
NPTS  IS  THE  NUMBER  OF  DATA  POINTS 

NVARS  IS  THE  NUMBER  OF  VARIABLES  TO  BE  READ  IN,  INCLUDING  DEPENDENT 
IDSET  IS  THE  INDEX  OF  THE  PARTICULAR  DATA  SET  BEING  CONSIDERED 
13SET  IS  THE  INDEX  OF  THE  FITTING  FUNCTION  LOOP  OR  SIMPLY  THE  B  SET 
I  NDE X 


CNT  ROL  ( 14 )  =  BLANK,  IF  ALL  DATA  ARE  IN  STNADARO  ARRANGEMENT  (SEE  WRITEUP) 


C 

c 

c 

c 

c 

c 

1 


1,  if  ONE  DESIRES  TO  READ  INDEPENDENT  VBLS  ONLY  ONCE  FOR  EACH 
SET  OF  B  PARAMETERS. 

2.  IF  THE  INDEPENDENT  VBLES  ARE  TO  BE  READ  ONLY  ONCE  FOR  EACH 
THE  ENTIRE  PROGRAM  IS  RUN. 

M  (  CNTROL (15)  >  IS  BLANK  IF  ALL  DATA  ARE  ON  CARDS.  OTHERWISE.  (M=I ) 

DATA  ARE  EITHER-ALL  ON  TAPE,  OR  ONLY  THE  DEPENDENT  DATA  ON  TAPE. 

IF  ( I ABS ( CNTROL! 9 ) ) .EQ. I )  PRINT  12.  NPTS 
DO  2  1=1 .NPTS 

READ  (5.13)  ( Z I J .  I  )  . J=1  .NVARS ) 

IF  ( I ABS ( CNTROL ( 9 ) ) .EQ. 1 )  PRINT  14.  I . ( Z ( J , I ) . J= 1 . N V AR S ) 

CONTINUE 


8 

9 


GO  TO  6 

J=2, NVARS ) 

I ) . 1=1 , NPTS ) 


GO  TO  11 
CONTINUE 
IF  ( IBSET.NE. I ) 

DO  5  1=1 . NPTS 

RE  AD  13,  ( Z (  J .  I  ) 

READ  { M, 1 3 )  ( Z ( ) 

GO  TO  11 
CONTI NUE 

IF  ( IDSET. NE. 1 )  GO 
GO  TO  4 

IF  ( I ABS! CNTROL ( 9 ) 1 .EQ. 1 ) 
NPT  S=  NP TSF l 


TO  6 


PRINT  12,  NPTS 


BRT21350 
I NDA0000 
I NDAOO 1 0 
I NDA0020 
INDA0030 
I NDA0040 
I NDA0050 
INDA0060 
I NDA0070 
I NDA0080 
I  NO A0090 
I NDAO 1 00 
I NDAO 1 1 0 
I NDAO 1 20 
I NDAO I  30 
I NDAO 140 
l  NDAO 1 50 
I NDAO 1 60 
I NDAO 1 70 
INDA01 80 
I NDA01 90 
TIME INDA0200 
I NDA02 1 0 
I NDA0220 
I NDA0230 
1 NDA0240 
I NDA0250 
I NDA0260 
INDA0270 
I NDA0280 
I NDA0290 
I NDA0300 
I NDA03 I  0 
I NDA0320 
I NDA0330 
I NDA0340 
I NDAO  350 
I NDA0360 
I  NO A 0370 
I NDA0380 
I NDA0390 
I NDA0400 


1  82 

1  83 

1  84 

1  85 

1  86 

1  87 

1  88 

1  89 

190 

191 

1  92 

193 

1  94 

1  95 

196 

197 

198 

199 

200 

201 

202 

203 

204 

205 

206 

207 

208 

209 

210 

21  1 

212 

213 

214 

215 

216 

217 

218 

219 

220 

221 

222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

2  36 

237 

238 

2  39 

240 

24  I 


lb 


READ  (5.13)  ( Z( J .NPTS) . J=1 .NVARS) 

IF  <Z( 1 .NPTS) .EO. 999. )  GO  TO  10 

IF  ( T ABS ( CNTROL ( 9 ) ) .EO . 1 )  PRINT  14.  NPTS. (Z(J.NPTS) .J=l .NVARS) 

GO  TO  9 

10  NPTS=NPTS-t 
CNTROL(3)=NPTS 

CALL  TO  MANIP  ALLOWS  MODIFICATION  OF  INPOT  DATA. 

11  CALL  MANIP  ( 3 .B , Z .MD I M, NO  I M, BD1 M) 

RETURN 

C 

12  FORMAT  (  1  HO , 1  DATA  READ  CHECK.  NPTS  =  *,I3/) 

13  FORMAT  (8E10.4) 

14  FORMAT  (1H  , 12 . • . • .8E12.4 ) 

END 

SUBROUTINE  SETUP  ( B . Z . A . SUMSQ . MD I M, ND I M. BD  I  M ) 

COMMON  CNTROL 

INTEGER  CNTROL ( 25 ). 8DI M 

REAL  A(BOIM.BDIM) . B < 8D I M , 4 > . Z < MD I M, ND I M) 

CALCULATE  THE  TERMS  IN  THE  MATRIX  EQUATION  A  »  X  =  C 

COMPUTE  ONLY  THE  SUM  OF  SQUARES  IF  CNTROL! 10)  NEGATIVE,  SINCE  THE  GAUSS 

CYCLE  DID  NOT  REDUCE  THE  SUM  OF  SQUARES . 

IF  (CNTROL ( 1 0 ) .LE .-1 )  GO  TO  2 
JJ=I ABS(CNTROL(  1  )  ) 

LL=CNTR0L(8) 

DO  1  J=l , J J 
B( J .3 ) =0. 0 
DO  1  K=1,JJ 

1  A( J . K ) =0 . 0 

2  NUMBER=I ABS( CNTROLI 3 ) ) 

SUM  SQ=  0 . 0 

I  HO LD= CNTROL (21  ) 

CNTROL ( 2 1 )=68 

CNTROL (21)  =  68  TELLS  YCOMP  THIS  IS  FIRST  CALL  OF  AN  ITERATION  OR  SUMSQ 

CALCUL AT  I  ON 

ZERO= 1 . / ( NUMBER- J J ) 

IF  (CNTROL ( 1 0 ). EO .-2 )  ZERO=l. 

DO  19  N=1 .NUMBER 

YC= YCOMP (N.B.Z.MDIM, ND1 M, BO I M ) 

CNTROL (21  )  =  I  HOLD 

IF  ( NUMBER . NE . I ABS ( CNTROL (3 )) )  GO  TO  2 
ZN=  Z ( 1  . N  ) 

DELY=ZN-YC 

GO  TO  ( 9, 3 ,4 , 5 . 6 . 7 , 8 ) ,  LL 

3  DEL Y  =  DEL  Y/ZN 
GO  TO  9 

4  DEL Y= DEL  Y/YC 
GO  TO  9 

5  DELY=DELY*ZN 
GO  TO  9 

6  DELY=DELY*YC 
GO  TO  9 

7  DELY=DELY/Z< MDI M.N) 

GO  TO  9 

8  DELY=DELY*Z(MDIM,N) 

9  SUMSQ=SUMSO+ ( DEL Y**2 ) /( ( NUMBER- JJ )* ZERO ) 

IF  (CNTROLI 1 0) .LE 1 )  GO  TO  19 

DO  17  K= 1 , J J 

B( K ,4 ) =DE  R I  V ( K  .N.B.Z.BDIM. MO IM.NDIM) 

GO  TO  (16.10.11.12.13.14,1b).  LL 

10  0(K,4 )=B(K ,4 )/ZN 


1 NDA04 l 0 

1NDA0420 

INDA0430 

I ND A  04  4  0 

I NDA0450 

I NDA0460 

INDA04  70 

I ND AO  4  80 

I NDA0490 

I NDA0500 

I ND AOS  1 0 

I N DAO 520 

I NDA0530 

I  N  D  AO  54  0 

SETUOOOO 

SETU0010 

SETU0020 

SETU0030 

SETU0040 

SE  TUO  0  50 

SETU0060 

SETU007  0 

SETU0080 

SETU0090 

SETU0100 

SETU01 1 0 

SETU0120 

SETU0130 

SETU0140 

SETUOl 50 

SETU0160 

SETUO 170 

SETU0180 

SETUO 190 

SETU0200 

SETU021 0 

SETU0220 

SETU0230 

SETU0240 

SETU0250 

SETU0260 

SET  UO  27  0 

SETU0280 

SETU0290 

SETUO 300 

SE  TUO  3 1 0 

SETU0320 

SETU0330 

SETUO  34  0 

SETU0350 

SETU0360 

SETU037  0 

SETU0380 

SETU0390 

SETU0400 

SETU0410 

SETU0420 

SETUOA  30 

SETU0440 

SETU0450 


2*2 

2*3 

2** 

2*5 

2*6 

247 

2*8 

2*9 

250 

251 

2  52 

253 

25* 

255 

256 

257 

258 

259 

260 

261 

262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

27* 

275 

276 

277 

278 

279 

2  80 

281 

282 

2  83 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 

296 

297 

298 

299 

300 

301 
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GO  TO  16 

It  B(K.4)=fHK.4  l/YC 

GO  TO  16 

12  B ( K  *  4  )  =  8(K,4 1*ZN 

GO  TO  16 

13  8 { K , 4 >  =  B< K.4 ) *YC 
GO  TO  16 

14  B ( K , 4 )  =0( K.4  1 /Z ( MDI M, N  1 
GO  TO  16 

15  8<K,41=B<K.*)*Z<MOIM,N) 

16  B(K,3)=B(K,3)*B(K,*)*0£LY 

17  CONTINUE 

DO  IB  J - 1 . J J 
00  18  K=J,JJ 

A(J,K)=A(J#K)+a(J.4)*B(K,4) 

18  CONTINUE 

19  CONTINUE 

IF  ( CNTROL! 1 0 > .LE .-1 )  RETURN 
00  20  K=2,JJ 
L=K  — 1 

DO  20  J= 1 ,U 

20  A ( K , J )=A( J,K 1 
RETURN 

END 

SUBROUT  I NE  BTEST  ( CNTROL . B . BO  I M ) 

INTEGER  CNTROL ( 25 ). BDI M 
REAL  B(BDIM,4> 

JJ= I ABSI CNTROL (11) 

XL  I  M=  I  . 00E-05 

IF  ( CNTROL (9  )  . EQ . 3)  XL I M= 1 . 00E-02 
DO  1  J=1,JJ 
DENOM= ABS(B<  J ,2 )  ) 

IF  ( OENOM.LT . 1 . OE-8 )  DENOM=1.0 
TEST=ABS(B(J.3)/DEN0M> 

IF  { TEST. GT.XLI Ml  RETURN 

1  CONTINUE 

IF  (CNTR0LI91 .EQ.3)  GO  TO  2 

2  CNTROL (21 )=1 
RETURN 

END 

FUNCTION  DERIV  ( K . N . B , Z . 8D I M , MD I M , ND I M > 

INTEGER  8DIM, CNTROL 

REAL  B( BDI M . 4 ) , Z ( MDI M. ND I M ) 

COMMON  CNTROL I  25 ) 

COMPUTE  FINITE  DIFFERENCE  APPROXIMATIONS  TO  PARTIAL  DERIVATIVES  IN 
CASE  THIS  ROUTINE  IS  NOT  SUPPLIED  BY  USER. 

I HOLD=C  NT  ROL (211 
CNTROL 121 1=67 

CNTROL (21)  =  67  SUPPRESSES  PRINTING  OF  CALCULATED  QUANTITY  YCOMP 

J  =  K 

H=0.001*ABS(8(J.l  )) 

REM  ARK  =  B ( J .  1  ) 

81 J. 1  )  =  REMARK  +H 

Y2  =  YCOMP{ N  »8t  Z  • MDI  Mt  NOI  M. BDI  M ) 

8 ( J • 1 >=REMARK-H 

Y 1 = YCOMP ( N  ,B . Z . MDI M. NOI M. BDI Ml 
B( J , 1 ) =REMARK 
DERIV=(Y2— Y1 1/(2. 0*H) 

CNTROL  <21  1  =  1  HOLD 
RETURN 


SE  TUO *60 
SE 1 UO*  7  0 
SETU0480 
SETU0490 
SETU0S00 
SE  TU  05 1 0 
SETU0520 
SE TUO 530 
SETU0540 
SETU0550 
SETU0560 
SETU0570 
SETU0580 
SETU0590 
SETU0600 
SETU061 0 
SET  U0620 
SETU0630 
SETU064  0 
SET  U0650 
SETU0660 
SETU0670 
SETU0680 
SETU0690 
BTESOOOO 
BTES001 0 
BT  ESO  02  0 
BTES0030 
BTES0040 
BTES0050 
BT  E  SO  06  0 
BTES0970 
BTES0O80 
BTES0090 
BTESO 100 
BTES01 1 0 
BTES0120 
BTES01 30 
BTESO 140 
BTESOl 50 
DER l 0000 
DERI  0010 
DER I  0020 
DER I  0030 
DER I  004  0 
DER 10050 
DERI  0060 
OER 10070 
OERI0080 
DER I  0090 
DERI  0 100 
DER I  01 1 0 
DERI  0120 
0ERI01 30 
DER 10140 
DER 101 50 
DER 101 60 
DERI  01 70 
DER 101 80 
DERI  0190 


302 

303 

304 

305 

306 

307 

308 

309 

310 

31  1 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 

332 

333 

334 

335 

336 

337 

338 

339 

340 

34  1 

342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 
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GO  TO  1 


CNTROL ( 20  > 

»  JJ 

J.  C  (  J  I  .  (  A<  J  ,K  I  .  K  = 


l  .  J  J  ) 


AT  ITERATION  13.59H,  IN  THE  MATRIX 
C  AND  A  ARE  / 1 1 5H  J  C(J) 

A  (  J  ,  K  )  *************************** 


END 

SUBROUTINE  DMPMAT  ( CNTROL . M . A , C > 

INTEGER  CNTR0LI25) 

RFAL  A ( M , M )  ,C<  M ) 

J J= I ABSl CNTRUL (1)1 
IE  ( CNTROL ( 2 1  > • EQ.  70  > 

PRINT  4,  CNTROL (20) 

GO  TO  2 
PRINT  5, 

DO  3  J=1 
PRINT  6. 

PRINT  7 
RETURN 

FORMAT  ( 1 5H0 
1  THE  POWS  OF 

2*********** 

34****  ********  ) 

5  FORMAT  (•HO  AF TER  '  I N VER S I  ON  AT  ITERATION  •.14,* 

1  AT  I  ON  AX  =  C.  THE  ROWS  OF  C  AND  A  ARE’/ 

3115H  J  C(J)  ***************** 

4************************************************ 

FORMAT  (1H  I  2. El 5.5.E20 . 5, 4E 1 5.5/IE38.5,  4E1  5.5 )  ) 
FORMAT  (IX) 

END 

SUBROUTINE  GAUSS  (  B  .  Z  .  MD  I  W  .  NO  I  M.  BD  I  M  .  A  ) 

COMMON  CNTROL 

INTEGER  CNTROLI25) ,BOIM 

REAL  A(SDIM,BDIM),B(BDIM,4),Z(M0IM,NDIM),S{24,6) 
BIN(K,J)  =  0.5-0.5*(-l.  )**(K/2A*(J— 1 )  ) 

IF  (CNTR0H22)  .  EO  .  0  .  )  GO  TO  1 
DETERM= 1 . 

CNTROL(22)=0 

NN=0 

NOR=0 

NWW  =  0 

NWFLAG=0 

NVSTQR=0 

PCHGSS=0. 

IF  (CNTROL ( 1 0 ) .LE . 0 )  GO  TO  2 
FAC=0. 1 4CNTR0L ( 10) 

FACLI M=CNTROL (10) 

GO  TO  3 
FAC=1 . 0 
FACLIM=FAC 
OLFAC=0 . 

J J=CNTROL { 1 ) 

L  IM I T=CNTROL ( 6) 

L  =  CNT  ROL ( 7 ) 

CNTROL ( 2  0 ) =NN 
CNTROL (21  ) =6  9 
DO  4  J=1,JJ 
B( J . 2 )=B( J . 1 ) 

SCALE=0 . 5 

MAR  K 1 0=CNTR0L (10) 

CALL  SETUP  ( B. Z , A , SUMSQ, MOI M, ND I M.BDI M) 

IF  (DETERM. NE. 0. )  GO  TO  6 
PRINT  43 

CALL  DMPMAT  ( CNTROL, BDIM, A.0( 1 .3) ) 

CNTROL (21 ) =4 
GO  TO  19 


DER I  0200 
DMPMOOOO 
DMPMOO 1 O 
DMPM0020 
DMPM0030 
DMPM0040 
DMPM0050 
DMPM0060 
DMPM0070 
D  M  DM0 08  0 
DMPM0090 
DMPMO 100 
DMPMO 1 l 0 

EOUATION  A  X  =  C , DMPMO 120 
***  ****** QMPMO 1 30 

»*****************OMPM0140 

DMPMO 150 

IN  THE  MATRIX  EQUDMPMO 160 
DMPMO 1 70 

***  X(J.K)  *****DMPM01 30 


*  *  **  * 


) 


DMPMO 190 
DMPM0200 
DMPM02 1 0 
DMPM0220 
GAUSOOOO 
GAUS0010 
GAUS002  0 
GAUS0030 
GAUS0040 
GAUS0050 
GAUS0060 
GAUS0070 
G AUSOOQO 
GAUS0090 
GAUS0100 
GAUS01 1 0 
G AUSO 120 
GAUS0130 
GAUS0140 
GAUS01 50 
GAUSO 160 
G AUSO 170 
GAUS01 80 
GAUSO 190 
GAUS0200 
GAUS021 0 
GAUS0220 
GAUS0230 
GAUS0240 
GAUS0250 
GAUS0260 
GAUS0270 
G AUS0280 
G AUS0290 
G AUS0300 
GAUS03 1 0 
GAUS0320 
G AUS0330 
GAUS0340 
GAUS0350 


’ 


362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 

375 

376 

377 

378 

379 

380 

381 

382 

383 

384 

385 

386 

387 

388 

389 

390 

391 

392 

393 

394 

395 

396 

397 

39  8 

399 

400 

40  1 

402 

403 

4  04 

405 

406 

407 

408 

409 

410 

41  1 

4  1  2 

4  1  3 

414 

41  5 

416 

4  1  7 

4  1  8 

4  1  9 

420 

421 


o  o 
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6 

IF  ( CNTROL ( 1 0 ) .LE .- 1 . AND. B IN ( L. 4 ) .EQ . 1 . )  GO  TO  7 

GAUS0360 

IF  (CNT  ROL (10). LE .  —  1  )  GO  TO  2 6 

GAUS0370 

IF  (BIN (L , 3) . EQ . 1 . )  CALL  OMPMAT  ( CN THOL . BD I M , A , B ( 1 

.3  )  ) 

GAUS0380 

7 

IF  ( B I N ( L , 2 ) . EQ . 1 . )  PRINT  42 

GAUS0390 

IF  (NN.EQ.O)  SQMIN=SUMSO 

GAUS0400 

IF  ( B I M L , 2 ) . FQ . 1 . )  PRINI  44.  NN , SUM SQ . ( B ( J , 1 ) , J = 1 

.  JJ) 

GAUS04 1 0 

IF  ( CNTROL ( 1 0 ). LE.-l )  GO  TO  26 

GAUS04ZO 

IF  ( CNTROL ( 20 )  .LE  .L I M I T  )  GO  TO  8 

GAUS04 30 

CNTROL (21 1=3 

GAUS0440 

GO  TO  19 

GAUS0450 

CALCULATE  SOLUTION  FOR  NORMAL  EQUATIONS 

GAUS0460 

R 

IF  ( JJ.NE. 1 )  GO  TO  9 

GAUS0470 

B( 1 , 3 ) =0 ( 1 , 3 ) / A ( 1,1) 

GAUS0480 

A(1 , 1 )=1 . / A ( 1,1) 

G AUS0490 

GO  TO  10 

GAUS0500 

9 

CALL  MATINV  ( A , J J . B ( 1 , 3 ) . 1 , DE TERM , BD I M ) 

GAUS05 1 0 

10 

IF  (BI NIL , 3) . EQ. 0. )  GO  TO  11 

GAUS0520 

l  HOLD=CNTROL  (  2  1  ) 

CNTROL (2 1 )  =  70 

C  NTROL (21)  =  70  HAS  THE  INVERSE  MATRIX  PRINTED  OUT  OF  OMPMAT 

CALL  OMPMAT  ( CNTROL, BDI M, A ,8 ( 1 .3 ) > 

CNTROL (21 1=1 HOLO 
1  1  NN=NN+ 1 

CNTROL (20 )=NN 
RESCAL=0.5 

IF  ( DETERM. EQ .0. )  GO  TO  5 
DO  12  1=1 , JJ 

S  (  I  .  4 ) =- 1  . 

If"  (0(1, 31. LT.0..AND.B( I t 1).LT.0..0R.8(I,3).GT.0. .AND. B{ 1,1] .GT.O 
11  S( I ,41= 1 . 

C..  S ( I , 2 )  CONTAINS  THE  PERMITTED  CHANGE 
C..  S ( I , 3 1  CONTAINS  THE  PREDICTED  LSTSQ  C 
..  S(I,4)=1.  IF  THE  LSTSQ  CORRECTION  IS 

..  AND  =-l.  OTHERWISE. 

S ( I , 1 1=8 ( 1,1) 

S< I ,3)=B< 1,3) 

S< I , 6 )=Sl I .4 ) 

12  S(I.2)  =  (F  AC— OLFAC  )*S(I.1)*S(I,4) 

IF  (8IN(L,3) . EQ . 0 . )  GO  TO  13 
PRINT  45,  ( <S( J,  I  1 , 1  =  1 , 4 > , J=1 . JJ | 

CORRECT  THE  CURRENT  VALUES  OF  THE  PARAMETERS. 


AND  ITS  REASSESSMENTS 
:ORRECT I  ON 

OF  THE  SAME  SIGN  AS  S(J,1) 


1  3 
C.  , 


14 


CALL  BTEST  ( CNT ROL , B , BD l M ) 

IF  CONVERGENCE  DUE  TO  STEPPING  DOWN  B(J,3) 

IF  (CNTROLC21  ) . NE.  1  )  GO  TO  20 

IF  (( CNTROL ( 1 0 ) .EQ .- 1 ). AND .( FAC . NE. FACL I M ) )  GO 
IF  ( CNTROL ( 1 0 ). GE. 0 . AND. NVSTQR .NE . 0 )  GO  TO  16 
IF  { CNTROL( 1 0 ) .NE .-1  )  GO  TO  19 
CNTROL (21 )=5 

C..OVERRIOE  OPTION 

IF  (CNTROL!  12  )  .LE.O)  GO  TO  15 

CALL  OVER  I D  (L . NOR. SQM  I N, J J, A. S. BDI M ) 

IF  (NOR. GE. CNTROL! 12 ) )  GO  TO  19 
GO  TO  39 

C . . WRONGWAY  OPTION 

15  IF  ( CNTROL (  1 2  1.EQ.0 )  GO  TO  16 

CALL  WRNGWY  ( L , NW W , SO M I N , J J , A , S , B , N WFL AG . BD 1 M) 
IF  ( NWW.GT.  I  ABS( CNTROL(  12 )  )  )  GO  TO  16 
GO  TO  5 

C.. LABELS  940,960,1030  AND  1040  WERE  COMMENTS 

C  ..  PARAMETER  CHANGE  OPTION 


AND  NOT  FITTING.  OVERRIDE 


TO  27 


GAUS0530 

GAUS0540 

GAUS0550 

GAUS0560 

GAUS057  0 

GAUS0580 

GAUS0590 

GAUS0600 

GAUS061 0 

GAUS0620 

GAUS0630 

GAUS0640 

GAUS0650 

GAUS0660 

GAUS0670 

GAUS0680 

GAUS0690 

GAUS0700 

GAUS07 1 0 

GAUS0720 

GAUS0730 

GAUS0740 

GAUS0750 

GAUS0760 

GAUS0770 

G AUS0780 

GAUS0790 

GAUS0800 

GAUS081 0 

GAUS0820 

GAUS0830 

GAUS084 0 

GAUS0850 

GAUS0660 

GAUS0870 

GAUS0880 

GAUS0890 

GAUS0900 

GAUS091 0 

GAUS0920 

GAUS0930 

GAUS0950 


42? 

4?3 

424 

425 

426 

427 

4  28 

429 

430 

4  31 

4  32 

433 

434 

435 

4  36 

437 

4  38 

4  39 

440 

44  1 

442 

443 

444 

445 

446 

447 

446 

449 

450 

451 

4  52 

453 

4  54 

455 

456 

457 

458 

459 

460 

461 

462 

463 

464 

4  65 

466 

467 

468 

469 

4  70 

47  1 

472 

473 

474 

475 

476 

477 

4  78 

479 

480 

491 
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16  IF  (CNTROL  (  1  1  >  .EO.O  .  AND. CNTROL  (  1  2  )  .  EQ. 0)  GO  TO  17 
IF  (CNTROL ( 1 1 ) .EO.O )  GO  TO  19 

CALL  PCHANG  (  NVSTOR  .  SQM  I  N  ,  PCHGSS  •  J  J  .  S  .  A  .(J  .  HD  I  M  ,NN,  64  ,  G28  ) 

17  DO  1 8  J=1 . J J 

18  B ( J . 1  )  =  S  <  J  .  1  ) 

19  RETURN 

C  .  . RESTR A  1  NED  CHANGE  MECHANISM 

20  IF  (( MARK1 0 .EQ. 0 .) .UR .( CNTROL ( 10 > . EO 1 ) )  GO  TO  2A 
DO  22  J= 1 . J J 

IF  (  ABS<  B<  J  ,  3  )  )  .LE.  At?S(  FAC*S<  J  .  1  )  )  )  GO  TO  21 
B(J.3)=S(J.4)*FAC*S(J.l ) 

IF  ( ABS ( B ( J , 3 ) ) . EQ. AHS( S ( J . 1 ) ) . ANO. S ( J . 4 ) . EQ .- 1 • )  B(J.3)= 
1  GN  (  1  •  OE-  10.SU.1  )  ) 

GO  TO  22 
S(J  »  2 )  =  B( J ,3) 

S( J.4 )  =  0  • 

CONT I NUE 
DO  23  J=1,JJ 

IF  ( S( J.4 ) .NE.O. )  GO  TO  24 
CONTINUE 
F  AC=F  ACL  I M 
OO  25  J=1,JJ 
B(J.1)=(B(J.2)*8<J.3)> 

IF  ( CNTROL (10).NE.— 1)  CNTROL! 10 )=— 3 
GO  TO  5 

TEST= ABS(  ( SUMSQ— SQM  IN) /SOMI N ) 

IF  ( TEST. GT. 1 . OE- 1 0 )  GO  TO  29 
CNTROL (21 ) =2 

CNTR0L(21 )  =  2  SIGNALS  THAT  SUM  OF  SQUARES  CHANGES  NEGLIGIBLY. 

IF  (( CNTROL < 1 0 )  .EQ.- 1  )  .AND . (FAC . NE .FACLI M )  )  GO  TO  27 
PRINT  46 

IF  ( NVSTOR.EQ.O)  RETURN 
GO  TO  16 
RESCAL=0.5 
CNTROL (21 ) =69 

SET  RESTRAINED  CHANGE  FACTOR 
OLF  AC  =  F  AC 
FAC=2 . »FAC 

IF  (FAC.GT.FACLIM )  FAC=FACLIM 
CNTROL( 10 ) =MARK1 0 
NN=  NN+ 1 
DO  28  J=1 . J J 
8(J»3)=S(J»3) 

S( J .2 )= (FAC— OLF AC )*S(J.1)*S(J.4) 

GO  TO  20 

RESTRAINED  CHANGE  PARTITION  CONVERGENCE  TEST 


2  1 


22 


23 


24 

25 


26 


27 


28 


29 


30 


31 


IF  ( SUMSQ. LT .SQM  I N)  GO  TO  38 

IF  ( NWFLAG.LE. 0 )  GO  TO  30 

NWFL  AG=  0 

GO  TO  39 

CNTROL ( 10 )=-l 

D IFF AC=FAC-OLFAC 

SCAF  AC=RESCAL  »D IFF AC 

rescal=scale»rescal 

IF  (  ( MARK1 0 . NE.O )  .OR . B I N ( L . 4 ) . NE • 1  )  GO  TO  33 
DO  31  J=1 . JJ 

B ( J , 3 ) =SCALE*B( J . 3 ) 

CALL  BTEST  (CNTROL.H.flDIM) 

IF  ( CNTROL ( 2 1 ). NE . 1 )  GO  TO  32 
IF  ( B I N( L . 4 1 . EQ . 1 . )  PRINT  47 


GAUS0970 
GAUS0980 
GAUS0990 
G  AUS 1000 
GAUSI010 
GAUS1020 
G AUS1 050 
GAUS 1 060 
C.AUS1  070 
G AUS1 080 
GAUS 1 090 
B(  J  .  3)  +SI GAUSl  100 
GAUS 1 1 10 
GAUS 1 120 
GAUS 1 1 30 
G  AUS 1 140 
GAUS 1 150 
G  AUS 1  160 
GAUS 1 1 70 
GAUS l 180 
GAUS1 1 90 
GAUS 1200 
GAUS 1210 
GAUS 1 220 
GAUS 1 230 
GAUS1240 
GAUS 1 250 
GAUS1260 
GAUS 1 27  0 
GAUSl 230 
G  AUS 1 29  0 
GAUS 1 300 
GAUS 1310 
GAUSl 320 
GAUSl 330 
GAUSl 340 
GAUSl 350 
GAUSl 360 
GAUSl 370 
GAUSl 380 
GAUS 1 390 
GAUSl 400 
GAUS 1 4 10 
GAUSl 420 
GAUSl 430 
GAUSl 440 
GAUSl 450 
GAUS 1 460 
G  AUS 14  70 
GAUSl 480 
GAUS 1490 
GAUS 1 500 
GAUSl 51 0 
GAUSl 52  0 
GAUSl 530 
GAUSl 540 
GAUS 1 550 
GAUSl 560 
GAUS 1570 
GAUSl 580 


/ 


482 

483 

484 

485 

486 

487 

488 

469 

490 

491 

492 

493 

494 

495 

496 

497 

498 

499 

500 

50  1 

502 

503 

50  4 

505 

506 

50  7 

508 

509 

51  0 

51  1 

512 

513 

514 

515 

516 

517 

518 

519 

520 

521 

522 

523 

524 

525 

526 

527 

528 

529 

53  0 

531 

532 

533 

534 

5  35 

5  36 

537 

538 

539 

540 

54  1 


157 


GO 

TO 

1  4 

32 

IF 

(BINIL.  .4) 

GO 

TO 

24 

33 

00 

34 

EQ.l.)  PRINT  48,  SCAFAC.QLFAC 


S(J,2)=SCALE*S(J,2) 

8( J .3  I =S( J ,2 ) 

34  CONTINUE 

I NTER=CNTROL( 9 ) 

IF  (OLFAC.EO.O)  GO  TO  35 
CNTR0L(9)=3 

CNTR0L19)  =  3  INVOKES  AN  ALTERNATE  CONVERGENCE  LIMIT  IN  3TEST 

35  CALL  BTEST  (CNTROL.B.BOI «l 
CNTR0LI9 )= INTER 

00  36  J=1 , J J 

IF  ( A8S( S( J , 3 ) ) .GT.F AC*S( J . 1 > > 

36  CONTINUE 

IF  (CNTROL(21 l.NE.l )  GO  TO  37 
IF  (BIN1L.4) .EQ.l . )  PRINT  47 
GO  TO  14 

37  IF  ( BINIL , 4)  ,EO. 1  . )  PRINT  49, 

GO  TO  13 

38  SOM  I N  =  SUM  SQ 
NOR=0 
NWW=0 
NWFLAG=0 

39  IF  ( MARK  10, GE, 1  )  F AC  =  0 . 1  * MARK  1 0 
NN=CNTROL ( 20 ) 

OLFAC=0. 

DO  40  J=1 . JJ 

40  B( J .2  >  =  B ( J , 1  ) 

IF  I NVSTOR. NE. 0 . OR.CNTROL ( 1 1 ) .EQ. 1 )  GO  TO  41 
CNTROLI 1 0)=MARK10 
GO  TO  5 

41  CNTROL(10)=0 
GO  TO  5 


B(J,3)=S(J,2) +0LFAC4S (J.l)*S(J.4) 


SCAFAC.QLFAC 


GAUS1S90 
GAUSI600 
GAUSl 61 0 
GAUS 1620 
GAUS1630 
GAUSl 640 
GAUS 1650 
GAUS 1 660 
GAUS1670 
GAUS 1680 
GAUS1690 
GAUSl 700 
G  AUS 17  10 
GAUSl 720 
GAUSl 730 
GAUSl 740 
GAUS 1 750 
GAUSl 760 
GAUSl 770 
GAUSl 780 
GAUSl 790 
GAUS! 800 
GAUS 1810 
GAUS1820 
GAUS 1 830 
GAUSl 840 
GAUS 1 850 
GAUSl 860 
GAUSl 870 
GAUSl 880 
GAUS 1 890 
GAUS1900 
GAUS1910 
GAUSl 920 
GAUSl 930 


C  GAUS 1940 

42  FORMAT  (11 4H0CYCLE  SUM  OF  SQUARES  ******** *********»******GAUS1950 

43  FORMAT  (*0  THE  FOLLOWING  NORMAL  EQUATIONS  ARE  SINGULAR  •)  GAUS1970 

44  FORMAT  ( 1H0I 2.2E21 . 7 , 4E 1 5 . 7/ ( E4S . 7.4E15.7) )  GAUS1980 

45  FORMAT  ('0  S  PARAMETERS  ' , // ( 4E 1 2 . 4 ) )  GAUS1990 

46  FORMAT  (*0  NEGLIGIBLE  CHANGE  IN  SUM  OF  SQUARES*)  GAUS2000 

47  FORMAT  ('0  PARTITION  LESS  THAN  MESH  LIMIT  -  -  -  RECYCLE  RESTRAINT  GAUS2010 

1  MECHANISM*)  GAUS2020 

48  FORMAT  (*0  NO  REDUCTION  IN  SUM  OF  SQUARES.  RESCALE  PARAMETERS  »//*GAUS2030 

1  RESTRAINED  CHANGE  MECHANISM  SET  AT  ( *,Ell.4,'+* ,F7.4, • )  TIMES  THEGAUS2040 

2  PREDICTED  LSTSQ  CORRECTION.')  GAUS2050 

49  FORMAT  (*0  NO  REDUCTION  IN  SUM  OF  SQUARES.  RESCALE  PARAMETERS  *//'GAUS2060 

1  RESTRAINED  CHANGE  MECHANISM  SET  AT  < • , E l 1 . 4 , • + * . F7 . 4 , • )  TIMES  THEGAUS2070 

2  VALUE  OF  THE  ORIGINAL  PARAMETER')  GAUS2080 

END  GAUS2090 

SUBROUTINE  OVERID  ( L , NOR . SOM  I N , J J , A , S . B , BD I M )  OVEROOOO 

COMMON  CNTROL  OVEROOIO 

INTEGER  CNTROL( 25 )  , BO  I M  OVER0020 

DIMENSION  A( 8DI M ,BDI M ) ,  S(BDIM.6).  P(24),  BIBOIM.4)  OVER0030 

NOR=NOR+ 1  OVER0040 

IF  (NOR.GT.CNTROLl 12 ) )  GO  TO  2  OVEROOSO 

IF  (BINIL .2 ) .EO . 1  . )  PRINT  4.  NOR . SQM  I N  OVER0060 

IF  (NOP.NE.l)  RETURN  OVER0070 

DO  1  J=1.JJ  OVER0080 


•* 


' 


542 

543 

544 

545 

54 6 

54  7 

548 

549 

550 

551 

552 

553 

554 

555 

556 

557 

558 

559 

560 

561 

562 

563 

564 

565 

566 

567 

568 

569 

570 

571 

572 

573 

574 

575 

576 

577 

578 

579 

580 

581 

582 

583 

584 

585 

586 

587 

588 

589 

590 

59  1 

592 

593 

594 

595 

596 

597 

598 

599 

600 

601 


1  J  C 


SI  J.5)  =  AI J, J) 

OVER  009  0 

1 

PI J )  =  S<  J .  1  ) 

OVERO 100 

I  CYC=CNTROH  20  )-l 

OVERO I 1 0 

RETURN 

OVERO 1 20 

2 

DO  3  J=1 ,JJ 

OVERO 1 30 

AI J , J ) =SI J .5) 

OVERO 1 40 

3 

BI J . 1 )=P( J) 

OVERO 1 50 

NOR=NOR- 1 

OVEROlBO 

IF  IBINIL.2) .EQ. 1 . )  PRINT 

5*  NOR. ICYC 

OVER  0 1 70 

RETURN 

OVERO 1 80 

C 

OVERO 1 90 

4  FORMAT  1*0  NO  CONVERGENCE 

IN  SCALING  INTERVAL. 

OVERRIDE  MECHANISM  OVER0200 

1  NO*  OPERABLE  -  OVERRIDE  * ,13.*  SOM  OF  SQUARES  MINIMUM  IS'.E14.6) 
5  FORMAT  (*0  AFTER  ' •  1 3 • •  OVERRIDES  REFERENCE  ITERATION  ‘,14) 

END 

FUNCTION  BIN  (M.N) 

BIN=0.5-0.5*(-l. >  »«t  M/2** IN-1  )  ) 

RETURN 

END 

SUBROUTINE  WRNGWY  ( L . NW W , SQM  I N , J J . A , S , B . NWFL AG . BD I M) 

COMMON  CNTROL 

INTEGER  CNTROLI 25 > . BO  I M 

DIMENSION  A( BDI M, BD I  Ml .  S124.6),  P<24),  B(BDIM,4) 

IF  ( BI N( L . 2 ) .EQ . 1 . )  PRINT  7,  SQM I N 
IF  INWW.NE.O)  GO  TO  2 
DO  1  J=1,JJ 
S<J,5)=A< J. J) 

1  PI J)=S( J, 1 ) 

I C YC=CNT  ROLI20I  —  1 

2  SM AX=ABSI SI  1 . 3 )/S(  1  , 1  )  ) 

DO  3  J=2,JJ 

IF  I ABS I  St J . 3 )/Sl J ,  1  >  ) .GT .SMAX)  S M A X  =  AB S l S I J . 3 ) /S I J .  1 )  ) 

3  CONTINUE 

DO  4  J=1 . J J 

SFR  AC  =  ABS I S I J,3)/ISMAX*SIJ,1))) 

4  BI J,  1 )  =  S( J ,  1 1-0 . 0S*S< J .6) *St  J, 1  ) * SFR AC 
CNTROLI 1 0 )=-l 

IF  INWW. GE. I ABSl CNTROLI 1 2 ) ) )  GO  TO  5 
NWW=NWW+  1 

C..NWFLAG  LIMITS  WRONGW AY  PARTITION  TO  ONE  RESTRAINED  CHANGE  CYCLE 
NWFLAG=l 
RETURN 

5  DO  6  J=1,JJ 
AI J, J)=SI J.S) 

6  BI J, 1 )=Pt J) 

IF  I  BIN! L.2 ) • EQ. 1  . >  PRINT  8,  ICYC 
RETURN 
C 

7  FORMAT  t'O  NO  CONVERGENCE  IN  SCALING  INTERVAL.  WRONGWAY  MECHANISM 
1  NOW  OPERABLE  SUM  OF  SQUARES  MINIMUM  IS  '.£14.6) 

8  FORMAT  t'O  WRONGWAY  LIMIT  EXCEEDED  REFERENCE  I TE R AT  I  ON'  .  I  4 ) 

END 

SUBROUTINE  PCHANG  l NV5T OR . SO M I N , PCHGSS , J J , S . A , 8 . BO I M . NN , * . A ) 

COMMON  CNTROL 

INTEGER  CNTROL I  25)  . 3DI M 

DIMENSION  SI24.6).  P I  24 ) •  BIBOIM.l),  AI BUI BO  I  Ml 

IF  ICNTROLI21  )  .EQ. I  . OR . CNTROL I  2 1  )  .EQ.2)  CN TROL l 2 1 ) =69 

IF  I  NVSTOR . NE. 0 . AND. SQM  I N .GE . PCHGSS )  GO  TO  2 

DO  1  J=1 . J J 

CONI 3-CNT ROL 113) 


OVER0210 

CVER0220 

OVER0230 

B I N00000 

BI N000 1 0 

BI N00020 

BIN00030 

WRNG0000 

WRNGOO 1 0 

WRNG0020 

WRNGOO 30 

WRNG0040 

WRNG0050 

WRNG0060 

WRNG007  0 

WRNG0080 

WRNG0090 

WRNGO 1 00 

WRNGO 1 1 0 

WRNGOl 20 

WRNGO 1 30 

WRNGO 140 

WRNGOl 50 

WRNGOl 60 

WRNGOl 70 

WRNGO 1 80 

WRNGOl 90 

WRNG0200 

WRNG0210 

WRNG0220 

WRNG0230 

WRNG0240 

WRNG0250 

WRNG0260 

WRNG0270 

WRNG0280 

WRNG0290 

WRNG030  0 

WRNG031 0 

WRNG0320 

PCHAOOOO 

PCHAOO 1 0 

PCHA0020 

PCHAOOOO 

PCHA004  0 

PCHA0050 

PCHA0060 

PCHA007  0 


602 

603 

604 

605 

606 

607 

600 

609  ■ 

61  0 

611 

612 

613 

614 

615 

616 

617 

61  8 

619 

620 

621 

622 

623 

624 

625 

626 

627 

628 

629 

630 

631 

632 

633 

634 

635 

636 

637 

638 

639 

640 

64  1 

642 

643 

644 

645 

646 

647 

648 

649 

650 

651 

652 

653 

654 

655 

656 

657 

65Q 

659 

660 

661 


159 


t 

2 


3 

4 


5 


6 

7 


C 


C 

c 

c 


c 

c 

c 

1 

2 

C 

c 

c 


3 


NUMBF.E=  J  J 
S(J.5)=A( J,J) 

IF  {  CNTROL  (  1  2  )  .  C.E  .  0  .  OR. CNTROL  (  2  1  )  .  EQ  «  691  P(J»  =  S(J,1) 

CONTINUE 

IF  ( SOM  I N. NE.PCHGSS 1  GO  TO  5 
PRINT  8 

CN1 ROL (I )=IABS(CNTR0L(2>) 

NVSTOR=I STORE 

IF  ( NVSTOR.NE.O 1  GO  TO  4 

DO  3  J=1  .NUMflEE 

CNTROL (13) =C0N1 3 

CNTROL ( 1 ) = NUMBEE 

J  J=  NU  MBEE 

A( J . J) =S( J .5) 

B(J. 1 ) =P ( J ) 

RETURN  l 
NVST  OR  =  0 
J J=CNTROL ( 1 ) 

CNTROL ( 1 1 ) =0 
CNTROL (  1 2 )  =  0 
RETURN  2 
PCHGSS=SOMIN 
I STORE=CNTROL( 1 ) 

IF  (CNTROL ( 1 ) . NE. I ABS( CNTROL! 2) ) )  GO  TO  6 
CNTROL ( 1  )  =  IA3S(CNT  ROL ( 2 ) )-l 
GO  TO  7 

CNTROL  (D-IABSI  CNTROL  (  2  )  ) 

NVSTOR=I STORE 
NV  A 8= I A8S( NVSTOR ) 

NE  WNV= I ABS(CNTROL(  1  )  )' 

PRINT  9,  NVA8.NEWNV 
NN=CNTROL ( 20 ) 

J  J=C  NTROL ( 1 ) 

RETURN  2 

8  FORMAT  (*0  PARAMETER  CHANGE  OPTION  UNAVAILING*) 

9  FORMAT  ( *  0 THE  NUMBER  OF  LSTSO  PARAMETERS  HAS  BEEN  CHANGED  FROM* 

1 , •  TO* , 13) 

END 

SUBROUTINE  M  AT  I NV  (  A , N , 8 , M , DE TERM , I D I M ) 


PCHA0080 
PCHA0090 
PCHAO 100 
PCHAO 1 1 0 
PCHAO 1 20 
PCHAO 1 30 
PCHA01 40 
PCHAO 1 50 
PCHAO 1 60 
PCHAO 170 
PCHAO 1 80 
PCHAO 190 
PCHA0200 
PCHA021 0 
PCHA0220 
PCHA0230 
PCHA0240 
PCHA0250 
PCHA0260 
PCHA0270 
PCHA0280 
PCHA0290 
PCHA0300 
PCHA031 0 
PCHA0320 
PCHA0330 
PCHA0340 
PCHA0350 
PCHA0360 
PCHA0370 
PCHAO30O 
PCHA0390 
PCHA0400 
PCHA04 1 0 
PCHA0420 
PCHA0430 
I  3PCHA044  0 
PCHA0450 
PCHA0460 
MAT  I  0000 
MAT  10010 


..  THE  DIMENSIONS  OF  IPIVOT, PIVOT  AND  INDEX  MUbT 

DIMENSION  A(  I D I M .  ID  I M ) .  8(IDIM,IDIM> 

DIMENSION  I P I VOT ( 24 ) ■  INDEXI24.2),  PIVOT(24) 
EQUIVALENCE  ( I ROW , J ROW ) .  ( I C OLUM . JC OLUM ) ,  (AMAX.T, 

INITIALIZATION 

DETERM=1  .0 
DO  2  J=1 . N 
IPI VOT ( J  >  =  0 
DO  22  1  =  1  .N 

SEARCH  FOR  PIVOT  ELEMENT 

AMAX=0.0 
DO  7  J=1 • N 

IF  (IPI VOT ( J )  —  1  )  3,7,3 

00  6  K=1.N 


BE  EQUAL  TO  BDIMMAT10020 
MAT  I  0030 
MAT  10040 
MATI 0050 

SWAP)  MAT  10060 

MAT  I  0070 
MAT  10080 
MAT  I  0090 
MAT  I  0 1 00 
MAT  101  10 
MAT  10  120 
MATI01 30 
MAT  I  0 1 40 
MAT  I  01 50 
MATI0160 
MATI 0 1 70 
MAT! 0100 
MATI 0190 
MAT  10200 


. 


662 

663 

664 

665 

666 

667 

668 

669 

670 

671 

672 

673 

674 

675 

676 

677 

678 

679 

660 

681 

682 

683 

684 

685 

686 

687 

688 

689 

690 

691 

692 

693 

694 

695 

696 

697 

698 

699 

700 

701 

702 

703 

704 

705 

706 

707 

708 

709 

710 

711 

712 

7  1  3 

714 

715 

716 

717 

7  1  8 

719 

720 

721 


1  u  0 


IF  I  I  PI VOTIK )- 1  )  4,6,26 

4  IF  ( ABS! AM  AX l-ARS ( A ( J ,K )  )  )  5,6,6 

5  IROW=J 

I CC!LUM=K 
AMA  X  =  A  <  J  ,  K  > 

6  CONTINUE 

7  CONTINUE 

IF  (A»AX)  8,27,8 

8  I  PIVOT!  I COLUM )  =  I P I VOT!  ICOLUM >+l 

C 

C  INTERCHANGE  ROWS  TO  PUT  PIVOT  ELEMENT  ON  DIAGONAL 

C 

IF  (  I  ROW- I COLUM )  9,13,9 

9  DETERM=-OETERM 
DO  10  L= 1 , N 
SWAP= A ( 1 ROW.L ) 

A!  I  ROW ,  L  )  =  A  !  ICOLUM.L) 

10  A! ICOLUM , L >=SWAP 

IF  { M )  13,13.11 

11  DO  1 2  L=1 , M 
SWAP=B( IROW.L) 

B  !  I  POW,  L  )  =  E3! ICOLUM.L) 

12  B< ICOLUM, L)=SWAP 

13  INDEX! I , 1 )=IROW 
INDEX!  I , 2  )  =  I  C  OLUM 

PI VOT(I)=A! ICOLUM, ICOLUM) 

DETERM=DETERM*PI VOT! I ) 

C 

C  DIVIDE  PIVOT  ROW  BY  PIVOT  ELEMENT 

C 

A! ICOLUM, ICOLUM >=1 .0 
DO  14  L=1.N 

14  A! ICOLUM, L)=A( I COLUM , L ) /P I VOT ( I ) 

IF  I M)  17,17,15 

15  DO  16  L= 1 , M 

16  B! ICOLUM ,L )=B! I COLUM , L ) /P I VOT t I ) 

C 

C  REDUCE  NON-PIVOT  ROWS 

C 

17  DO  22  LI =1 , N 

IF  ! L 1— I COLUM )  18,22.18 

13  T=A(L1 . ICOLUM) 

A ! L 1 , ICOLUM)=0.0 
DO  19  L=1,N 

19  A ( L 1 ,L >=A IL 1 ,L >-A ! I COLUM, L ) *T 

IF  I M )  22,22,20 

20  DO  21  L=1 , M 

21  BIL1  ,L )  =  B!L 1  ,L >  —  B (  I COLUM, L > *T 

22  CONTINUE 

C 

C  INTERCHANGE  COLUMNS 

C 

DO  25  1=1, N 
L=N+1-I 

IF  ! INDFXIL, 1 )— INDEX1L.2) )  23,25,23 

23  J  RO W= I NDE  X ( L .  1  ) 

JCOLUM= INDEX (L . 2 ) 

DO  24  K=1,N 
SWAP= A! K . J  RO  W ) 

A(K,JR0W)=A(K, J COLUM ) 


MAT l 0? 1 0 
MAT  I  0220 
MAT  I  0230 
MATI 0240 
MAT  I  0250 
MATI 0260 
MA  T I  0270 
MAT  I  0280 
MATI 0290 
MAT  10300 
MAT  10310 
MAT  I  0320 
MAT  10330 
MATI 0340 
MATI 0350 
MAT  I  0360 
MATI 0370 
MAT  I  0380 
MATI 0390 
MATI 0400 
MAT  104 1 0 
MATI 0420 
MATI 0430 
MAT  I  044  0 
MATI 0450 
MAT  I  0460 
MATI 0470 
MAT  I  0480 
MAT  I  0490 
MAT  I  0500 
MATI 05 1 0 
MAT  I  0520 
MAT  I  0530 
MAT  I  054  0 
MATI 0550 
MATI 0560 
MAT  I  0570 
MATI 0580 
MATI 0590 
MAT  I  0600 
MATI 06 l 0 
MAT  I  0620 
MATI 0630 
MAT  I  064  0 
MATI 0650 
MATI 0660 
MATI 0670 
MATI 0680 
MATI 0690 
MAT  I  0700 
MATI 07 1 0 
MATI 0720 
MAT  10730 
MAT  I  074  0 
MATI 0750 
MATI 0760 
MATI 0770 
MAT  I  0780 
MATI 0790 
MATI 0800 


722 

723 

724 

72b 

726 

727 

728 

729 

730 

731 

7  32 

733 

7  34 

735 

736 

737 

738 

7  39 

740 

74  1 

742 

743 

744 

745 

746 

747 

748 

749 

750 

751 

752 

753 

754 

755 

756 

757 

758 

759 

760 

761 

762 

763 

764 

765 

766 

767 

768 

769 

770 

771 

772 

773 

774 

775 

776 

777 

778 

779 

780 

781 


161 


A (K, JCOLUM)=SWAP 

24  CONTINUE 

25  CONTINUE 

26  RETURN 

27  DETCRM=0. 

RETURN 

END 

SUBROUTINE  PRINDA  (IC.Z.M.N) 

DIMENSION  I C ( 25 ) ,  Z(M.N) 

PRINT  2 

NPTS  = I A9S  <  1C ( 3  I) 

NVAR= I C ( 5 ) 

DO  1  1=1 ,NPTS 

1  PRINT  3.  I . <Z< J. I ) . J=1 .NVAR) 

RETURN 

c 

2  FORMAT!  *0  INPUT  DATA'//) 

3  FORMAT  <I4,8E15.5> 

END 

SUBROUTINE  FINALE  (  8  .  Z  ,  MD  I  M.  ND  I  M  ,  ID  I  M  ,  A  ) 

COMMON  CNTROL 

DIMENSION  Z( MDI M . ND 1 M ) ,  B(IDIM),  A ( I D I M , I D I M 1 ,  VARI24) 

INTEGER  CNTROL ( 25 ) 

DIMENSION  FRCTOVI9) 

J J= I ABS( CNTROL l 1 ) ) 

NUMBER=I ABS( CNTROL! 3) ) 

AV-0.0 
AV 1=0*0 
AV2=0. 0 
YMAX=0 . 0 
ZMAX=0. 0 
ZZM AX  =  0  « 0 
SUM  SQ=  0 . 

CNTROL! 1 0 1 =— 2 

CNTROL! 1 0)=-2  DIRECTS  SETUP  TO  CALCULATE  VARIANCE  RATHER  THAN  SUMSQ 
CALL  SETUP  1 3 ,Z . A , SUMSQ .MDI M , NDI M , IDI M) 

DO  1  1=1 , JJ 

1  VAR! I )=SORT! A( I , I 14SUMSQ) 

SS= SUMSQ*! NUMBER- J J ) 

STDE V= SORT! SUMSQ ) 

PRINT  5,  CNTR0L120) , IBt J) . J=1 ,JJ) 

PRINT  6.  ! V AR! I ) , I =1 , J J > 

DO  2  1=1 .JJ 

2  FRCTDV! I  )  =  VAR(  I  >/B!  I  I 
PRINT  7,  (FRCTDV! I ) , 1=1 , JJ) 

PRINT  8,  SS, SUMSQ, STDEV 
PRINT  9 

IF  (CNTR0LI5  ).EQ.2)  PRINT  10 
IF  ( CNTROLI5 ) .NE. 2 )  PRINT  11 
•  DO  4  N= 1 , NUMBER 

•  YC=YCOMP ( N ,B , Z . MOI M , ND I M, IDIM) 

DEL  Y=  Z (  1  , NJ-YC 

IF  (Z(l.N).NE.O.O)  RAT  I 0=DELY/Z!  1  , N ) 

IF  I Z ( 1  , N ) .EQ .0 .0 )  RAT  I  0= 1  «  0E30 
ABSRAT=ABS(RATIO) 

AV= AV+DELY 
A  VI =A VI +RAT I  0 
A  V2  =  AV2 -f  ABSR  AT 

IF  (CNTR0L(5).EG.2)  PRINT  13,  N . Z ( 2 , N ) . Z ( 1  . N 1  . YC . DEL Y , R AT  I O 
IF  (CNTROL ( 5 > . NE. 2 )  PRINT  12,  N , Z ( 1 . N ) . YC . DEL Y , RA T I O 


MAT  I  08 10 
MAT  I  0820 
MAT  I  0830 
MAT  I  0840 
MAT l 0850 
MAT  10860 
MAT  I  0870 
PR  I NOOOO 
PR  I NOO 1 0 
PR! N0020 
PR  I NOO  30 
PR  I  NO 0 4  0 
PR  I N0050 
PR  I N00O0 
PR  I N0070 
PR  I  NO  0 80 
PR  I N0090 
PR  I  NO  1 00 
PR I  NO  1 10 
F I NAOOOO 
FI NA0010 
F I NA  0  02  0 
FI  NA0030 
F I NA0040 
F I NA  0  0  50 
F INA0060 
FI NA0070 
F I NA0080 
F I NA0090 
F I NAO 100 
F I NAO 1 1 0 
FI NAO 1 20 
F I NAO 1 30 
FINA01 40 
FI NAO 1 50 
F I NAO 160 
F I NAO 170 
F I NAOl 80 
FI NA0190 
F I NA0200 
FINA0210 
F 1 NA0220 
F I NA02 JO 
F  I  N  AO  24  0 
F  I  NA0250 
F I NA0260 
F1NA0270 
F I NA0280 
FI NA0290 
FI NA0300 
F I NA03 1 O 
F I NA0320 
F I NA0330 
F I NAO 340 
F I NAO  350 
FI NA0360 
F I NA0370 
F  I  N AO  380 
F  I  NAO 390 
F INA0400 


782 

783 

784 

785 

786 

737 

788 

789 

790 

791 

792 

793 

794 

795 

79  6 

797 

798 

799 

800 

801 

802 

803 

804 

805 

806 

807 

808 

809 

8  1  0 

81  1 

812 

8  1  3 

8  1  4 

815 

816 

817 

81  8 

81  9 

820 

821 

822 

823 

824 

825 

826 

82  7 

828 

629 

830 

831 

832 

833 

834 

835 

836 

8  37 

838 

839 

840 

84  1 


152 


c 

5 


1  i 

12 

13 

>4 


A0  S V  AL  =  ABS  (  DEL  Y 1 

IF  (  YMAX  .(it  •  A3SVAL  )  GO  TO  3 

YMAX= ABSVAL 

Y YMA  X  =  OELY 

MAHK=N 

IF  ( ZMA X . GT . ABSR AT )  GO  TO  4 

ZMAX=ABSRAT 

ZZM AX=RAT 10 

HAHK 1 =N 

CONTI  NUE 

D=NUMBER 

AV-AV/O 

A  V  l  =  A  V  1 /D 

AV2= AV2/D 

RTMNSO=SORT ( SUMSQ ) 

PPINT  14,  AV  ,  A  V  1  ,  AV2  ,  YYMAX  ,  MARK  ,  ZMA  X  ,  w  ARK.1  ,R  TMNSQ 
CALL  MANIP  ( 4 , 0, Z , MDI N, NO  I M,  ID  I M ) 

RETURN 


FORMAT  { 8H-  AFTER  I3.12H 


FI NA04 1 0 
FINA0420 
F I NA0430 
F I NA0440 
F I NA0450 
F I NA0460 
F I NA0470 
F I NA04  80 
F I NA0490 
F I NA0500 
F I NA051 0 
F I NA0520 
F I NA0530 
F I NA0540 
FI NA0550 
F I NA0560 
FI NA0570 
FI NA0580 
F I NA0590 

I TEPATI ONS/45HO  THE  MINIMIZING  VALUES  OFFINA0600 


'// 


1  THE  PARAMETERS  ARE// ( 1 X , 5E 1 5 • 7 ) ) 

6  FORMAT  ('0  WITH  THE  FOLLOWING  STO.  DEVIATIONS 
1 < 1X.SE15.7) ) 

7  FORMAT ( ■ 0  WITH  THE  FOLLOWING  FRACTIONAL  DEVIATIONS'// 

8  FORMAT  ('0  THE  SUM  OF  SOU ARES , V AR I  A NCE  AND  SAMPLE  STD 
1  FOR  THIS  FIT  ARE • 3E 16,5) 

FORMAT  ( 1 

1 

10  FORMAT  (1H0, 'NUMBER  X  OBSERVED 

1  ALCULATED  OBS  -  CALC  (OBS  -  CALC) 

FORMAT  (94H0NUMBER  Y  OBSERVED 

IBS  -  CALC  (OBS  -  CALC)  /  OBS  /) 

FORMAT  ( I5,E23.5,E17.5,2E19.5) 

FORMAT  ( I5,E23.5.2E17.5,2E19.5) 

FORMAT  ( 30H0  AVERAGE  DEVIATION  E14.5/30H 


FINA061 0 
F I NA0620 
FI NA0630 
1  X.5F15.3)FINA0t>40 
DEVIATION  FINA0650 
F I NA0660 

.**»»****t**PlNA0670 

NA0680 


Y  OBSERVED 
/  OBS  •/) 

Y  CALCULATED 


CFINA0690 
FI N A  0  7  00 
OF  I NA07 1 0 
FINA0720 
F I NA0730 
FI NA0740 
AVERAGE  RELFINA0750 


1  0EV  E14.5/30H  AVE  ABS  REL  DEV  E14.5/30H  MFINA0760 
2 AX  I  MUM  DEVIATION  E14.5.6X.10H  AT  POINT  I4/30H  MAXIMUM  RFINA0770 
-}EL  DEV  E14.5.6X.10H  AT  POINT  I4/30H  ROOT  MEAN  SQUARE  FINA0780 


4  DEV  I  AT  I  ON  E14.5) 

ENO 

SUBROUTINE  MANIP  ( M , B , Z , MQ I M , ND I M , I D I M ) 

CALLS  TO  THIS  SUBROUTINE  ALLOW  MODIFICATION  OF  THE  OATA 


C 

c 

c 

c 

c 


AND  DATA  ARE  READ  BY  THIS  SUB- 


M  =  1  AFTER  B'S  ARE  READ. 

M  =  2  IF  CNTROL ( 3 )  IS  LESS  THAN  0, 

ROUTINE. 

M  =  3  AFTER  DATA  ARE  READ  IN  BY  INDATA. 

M  =  4  AFTER  STATISTICS  HAVE  BEEN  COMPUTED  (CALL  OUTDAT ) 
COMMON  CNTROL 
INTEGER  CNTROL ( 25 ) 

DIMENSION  B(IDIM),  Z(MDIM.NOIM) 

GO  TO  (1.2.3,12).  M 

RETUPN 

RETURN 

NPTS=I ABS (CNTROL! 3) ) 

NVARS= I ABS( CNTROL! 5 ) ) 

DO  4  1=1,  NP T S 

STOR=Z ( 1 . I ) 

Z(1.I)=Z(2.I> 

Z  <  2 .  I  )=STOR 
Ml 0=NPTS/2 


F  I  NA 07  90 
F I NA0800 
MANI 0000 
MANI 0010 
MANI 0020 
MANI 0030 
MAN  10040 
MANI 0050 
MANI 0060 
MANI0070 
MAN  I  0080 
MANI 0090 
MANI 0 1 00 
MANI 01 1 0 
MANI0120 
MANI 01 30 
MANI 0140 
M  AN  10150 
MANI0160 
MANI 0 1 70 
MANI 0180 
MANI 0190 


842 

OO  5  1=1. NPTS 

M ANl 0200 

84  3 

5 

Z  (  2  ,  I  )  =  Z  (  2  ,  I  )  /  (  B( 4 ) /O ( 5  1  > 

MANI 02 1 0 

844 

IF  (6(2) )  7.6.7 

MAN  I 0220 

845 

6 

8<2>=Z(1.1H-<Z(1.1>-Z{1.2))*(Z(2,1>/(Z(2.2)-Z<2.11)) 

MANI 0230 

846 

7 

IF  (0(3) )  9,8,9 

MANI 0240 

847 

8 

B ( 3 ) =7 ( 1 . NPTS ) 

MAN  I 0250 

848 

9 

IF  (6(1))  11.10.11 

MANI 0260 

849 

1  0 

B(l>=-(AL0G((Z(l.M10)-B(3>)/{Z(l,l)-B(3>)>)/(Z<2,MI0)-Z(2.1) 

)  MANI 0270 

eso 

1  1 

PRINT  13,  (0(1),  1  =  1, 3) 

MAN l 0280 

851 

RETURN 

MANI 0290 

852 

1  2 

CALL  OUTDAT  ( 6 , Z . MD I M , ND I  M .  ID  I M ) 

MANI 0  300 

853 

RE  TURN 

MANI 03 I 0 

854 

C 

MAN  I 0320 

855 

13  FORMAT  (lHO.’THE  ESTIMATED  FIRST  ORDER  RATE  PARAMETERS 

are: 

•,  MANI 033  0 

856 

1 3G1 2.4) 

MANI 0340 

857 

END 

MANI 0350 

858 

FUNCTION  YCOMP  ( N , B . Z . MD I M , ND I M , BDI M ) 

YCUM0000 

859 

INTEGER  CNTROL ( 25 ) 

YCOMOO 1 0 

860 

COMMON  CNTROL 

YCOM0O20 

861 

INTEGER  BDIM 

YCOMOO  30 

862 

DIMENSION  8(  BDI M )  ,  Z ( MD I M . ND I M ) 

YCQM0040 

863 

1 

YCOMP= (B(2)-B(3))*EXP(-B( 1 )*Z(2,N))FB(3) 

YCOMOO 50 

864 

RETURN 

YCOM0060 

865 

END 

YCOM0070 

866 

SUBROUTINE  OUTDAT  ( 8 • Z , MD I  M  , NO  I  M. IDl  M) 

OUTDOOOO 

867 

DIMENSION  B(IDIM).  Z ( MD I M . ND I M ) 

OUTDOOl 0 

868 

COMMON  CNTROL 

OUTD0020 

869 

REAL  K2LQG 

OUTD0030 

870 

INTEGER  CNTROL ( 25 ) 

OUTD004  0 

87  1 

IF  ( CNTROL ( 2 )  . NL . 6. OR. DC  6 )  ,EQ .0 . 0 )  RETURN 

OUTD0050 

872 

K2LOG=ALOG10(S( 1 )/8(6)) 

OUTD0060 

873 

PRINT  1 ,  K2L0G 

OUTD007  0 

874 

RETURN 

OUTD0080 

87  5 

c 

OUTD0090 

876 

1  FORMAT ( 1  HO , •  THE  BASE  10  LOGARITHM  OF  THE  SECOND  ORDER 

RATE 

CONSTAOUTDOIOO 

877 

1  NT  IS:  •  .F6.31 

OUTDO  1 10 

878 

END 

OUTDO  120 

END  OF 

FILE 

